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I.  BACKGROUND 


Numerous  evaluations  of  the  acoustic  field  radiating  from  a 
baffled  transducer  have  appeared  in  the  published  literature.  An 
important  feature  is  that  these  theories  are  applicable  for  a 
wide  range  of  parameters.  Approximations/  such  as  those  describ¬ 
ing  an  axisymmetric  sound  beam  in  the  far  field  (Fraunhofer  zone) 
can  substantially  reduce  computational  cost#  but  they  are  not 
necessary.  Linear  theory  is  valid  when  the  source  level  is 
sufficiently  low.  Even  then,  diffraction  effects  in  the  near 
field,  which  lead  to  localized  cancellations  and  reinforcements, 
complicate  the  task  of  correlating  near  field  measurements  to  far 
field  propagation  properties.^ 

The  situation  becomes  more  complicated  when  one  tries  to 
increase  the  propagation  range  by  raising  the  source  level.  It 
is  logical  to  try  to  overcome  effects  such  as  dissipation  and 
scattering  by  generating  higher  level  signals.  Such  attempts 
inevitably  lead  to  a  greater  role  for  nonlinear  effects.  One  of 
the  effects  of  nonlinearity  is  to  divert  energy  from  the  fun¬ 
damental  signal  to  higher  harmonics,  which  is  equivalent  to 
lowering  the  efficiency  of  the  transducer.  In  the  face  of  these 
concurrent  effects  it  is  apparent  that  developing  a  unified 
theory  for  nonlinear  effects  in  sound  beams  is  a  challenging 
matter.  However,  such  a  theory  is  necessary  if  understanding  of 
the  distortion  phenomena  is  to  be  enhanced.  A  prime  example  of 
the  prior  lack  of  insight  is  the  observed  differences  between  the 
distortion  of  the  compression  and  rarefaction  phases  of  a  signal, 
which  had  no  analog  in  simpler  types  of  acoustic  waves. 

A  variety  of  approaches  have  been  employed  to  study  the 
effects  of  nonlinearity  in  this  system.  One  approach  has  relied 
on  a  conventional  perturbation  solution  of  an  approximate  non¬ 
linear  wave  equation.  Such  an  analysis  seems  to  give  very  good 
results  near  the  transducer  face.  However,  it  quickly  breaks 
down  with  increasing  range  due  to  assumptions  that  are  made  in 
the  perturbation  steps. 

An  investigation  of  properties  in  the  far  field  was 
developed  based  on  an  approximation  as  a  quasi-spherical  wave. 
Such  a  formulation  assumes  that  the  wave  arrives  at  the  transi¬ 
tion  to  the  far  field  (e.g.  the  Rayleigh  distance)  without 
substantial  prior  distortion.  Hence,  the  spherical  wave  descrip¬ 
tion  is  inherently  limited  to  cases  where  the  transducer 
excitation  is  comparitively  low  level.  This  type  of  analysis 
also  leads  to  certain  anomalies,  such  as  the  fact  that  the  level 
of  distortion  is  dependent  on  the  choice  for  the  spherical  tran¬ 
sition  distance,  which  may  be  arbitrarily  chosen  beyond  the 
Rayleigh  distance. 

Another  approach  that  has  been  widely  employed  in  the 
Russian  literature  is  founded  on  a  version  of  Burgers'  equation 
that  has  been  modified  to  account  for  spreading  and  diffraction. 
This  nonlinear  partial  differential  equation  has  been  solved 
numerically  for  several  types  of  boundary  motion.  The  primary 
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limitation  of  this  approach  are  the  approximations  on  which  the 
model  equation  are  based,  and  computational  difficulties  stemming 
from  rapid  transition  due  to  diffraction.  An  alternative  treat¬ 
ment  of  this  nonlinear  parabolic  equation  based  on  Fourier  series 
expansions  has  reduced  the  computational  problems,  but  quaestions 
regarding  the  adequacy  of  the  model  equation  still  remain. 


II.  RESEARCH  TECHNIQUE 

The  primary  goal  of  this  project  was  to  develop  an  analyti¬ 
cal  description  of  transducer  radiation  in  which  finite  amplitude 
effects,  diffraction,  and  spherical  spreading  are  treated  consis¬ 
tently,  without  limitation  to  a  speciific  spatial  domain.  The 
technique  employed  singular  perturbation  theory  in  conjunction 
with  asymptotic  analysis. 

The  general  approach  uses  the  King  integral  in  linear 
theory,  which  is  a  Fourier-Bessel  integral  transform,  to  develop 
the  second  order  source  terms  that  generate  nonlinearities  in  the 
response.  There  are  two  kinds  of  nonlinear  effects  that  arise  at 
the  second  order.  Some  produce  terms  that  remain  bounded  as  the 
signal  propagates.  (One  such  effect  is  associated  with  the  fact 
that  the  input  from  the  transducer  originates  from  a  moving 
boundary,  rather  than  the  much  simpler  description,  2=0.)  The 
smallness  of  the  acoustic  Mach  number  leads  to  the  conclusion 
that  these  fixed  magnitude  effects  cannot  account  for  measured 
levels  of  distortion.  The  other  group  of  nonlinear  effects  arise 
from  resonance-like  phenomena.  These  terms  lead  to  distortion 
that  grows  with  increasing  distance.  Shocks  ultimately  form  from 
this  effect,  unless  dissipation  is  adequate  to  overcome  the 
nonlinear  distortion  process.  It  is  this  cumulative  growth 
effect  that  needs  to  be  evaluated. 

The  growth  effects  in  the  second  order  terms  are  evaluated 
by  using  asymptotic  integration  techniques  to  identify  the  por¬ 
tion  of  the  second  order  terms  that  grow  most  rapidly  with 
increasing  range.  The  aforementioned  breakdown  of  conventional 
(i.e.  regular)  perturbation  solutions  is  avoided  by  introducing 
coordinate  transformations  that  essentially  are  based  on  recogni¬ 
tion  that  cumulative  growth  is  a  singularity. 

This  approach  was  the  basis  for  a  variety  of  studies.  Some 
were  devoted  to  developing  efficient  numerical  algorithms  for 
quantitative  evaluation.  Others  increased  the  generality  of  the 
transducer  vibration,  including  linkage  with  the  effects  of 
elastic  behavior.  Another  group  of  studies  endeavored  to  obtain 
insight  into  the  nature  of  physical  processes.  These  efforts  are 
surveyed  below. 


III.  RESEARCH  ACHIEVEMENTS 


A.  Axisymmetr ic  Monochromatic  Excitation 

The  analysis  described  in  the  previous  section  involved  a 
large  effort  to  place  some  of  the  steps  on  a  firm  mathematical 
foundation.  This  effort  has  also  clarified  the  physical  under¬ 
standing  of  the  manner  in  which  a  continuous  spectrum  of  modes 
interact  to  create  distortion.  The  initial  presentation  [2]  was 
improved  substantially  in  the  published  version  [6,  7].  Those 
works  compared  the  analytical  results  to  series  of  independent 
measurements.  The  predictions  for  amplitude  levels  of  the  har¬ 
monics  were  well  within  the  identifiable  experimental  error. 
Furthermore,  a  comparison  of  waveforms  showed  that  the  analysis 
does  describe  the  different  shape  of  the  waveform  in  the  compres¬ 
sion  and  arefaction  phases.  The  comparison  also  confirmed  a 
change  in  this  asymmetry  associated  with  the  transition  from  the 
near  field  to  the  far  field.  No  prior  analysis  had  anticipated 
this  phenomenon. 

The  original  version  of  the  computer  program  for  this  model 
was  quite  inefficient.  It  required  a  relatively  large  amount  of 
computer  memory  and  long  execution  times.  Of  even  greater  con¬ 
cern  was  the  dependence  of  the  analysis  on  a  hypothesis  regarding 
the  nature  of  the  distortion  process  far  from  the  axis  in  com¬ 
parison  to  very  close  to  the  axis.  Both  matters  were  effectively 
treated  by  a  recently  completed  Ph.  D.  thesis  [12,  14].  That  work 
developed  independent  solutions  in  each  region  by  using  different 
asymptotic  approximations.  The  signal  in  the  paraxial  region  was 
found  to  behave  like  a  spectrum  of  quasi-planar  waves,  whereas 
the  off-axis  signal  was  found  to  consist  of  spectra  of  inward  and 
outward  propagating  conical  waves.  The  individual  responses  were 
then  matched  to  obtain  uniformly  accurate  expressions. 
Significantly,  the  results  agreed  with  the  earlier  mathematical 
forms.  A  key  benefit  of  the  analysis  was  that  the  new  perspec¬ 
tive  led  to  a  Fourier  series  representation  that  decreased 
computational  time  by  as  much  as  a  factor  of  one  hundred. 

This  improved  computational  power  was  exploited  to  test  the 
theoretical  predictions  against  several  prior  series  of 
measurements.  The  agreement  was  generally  very  close,  with  one 
exception.  Experiments  by  Gould  in  the  late  1960 's  investigated 
the  signal  very  close  to  the  transducer  face  at  a  high  reduced 
frequency,  ka  =  114.  The  measured  second  harmonic  distribution 
was  quite  different  from  the  theory.  This  discrepancy  is  now 
under  investigation.  Current  indications  are  that  the  dis¬ 
crepancies  between  theory  and  measurement  correspond  to  the 
second  order  terms  that  are  discarded  in  the  analysis  because 
they  have  minor  influence  far  from  the  transducer. 

The  ease  with  which  computations  could  be  performed  made  it 
possible  to  generalize  the  nature  of  the  excitation.  A  study  of 
particular  interest  considered  the  case  where  the  transducer  is 
actually  a  membrane  that  is  subjected  to  steady-state  harmonic 
excitation  [11]  .  The  response  of  the  membrane  and  the  surface 


pressure  are  fully  coupled  in  this  situation.  The  analysis  of 
the  finite  amplitude  sound  beam  was  combined  with  a  vibration 
analysis  of  the  membrane  to  determine  all  aspects  of  the 
response.  It  was  shown  that  the  vibration  analysis  could  be 
safely  performed  without  recourse  to  nonlinear  theories  for  the 
acoustic  signal,  but  that  the  acoustic  signal  required  knowledge 
of  the  effects  of  nonlinear  elasticity  when  the  excitation  was  at 
a  system  resonance. 


B.  Nonsymmetric  Monochromatic  Excitation 

The  recent  analysis  of  axisymmetric  waves  [12]  also  con¬ 
sidered  a  situation  in  which  the  transducer  vibration  consists  of 
the  superposition  of  an  axisymmetric  component  and  an  azimuthally 
travelling  nonsymmetric  part.  The  latter  corresponds  to  a  spa¬ 
tially  phased  rocking  motion  that  resembles  the  wobbling  of  a 
rolling  coin  as  it  falls  to  the  ground.  The  analysis  was  much 
like  that  for  the  axisymmetric  case.  One  benefit  of  the  greater 
generality  was  a  clarification  in  the  task  of  matching  the  off- 
axis  and  paraxial  signals.  This  study  is  apparently  the  first  to 
address  nonsymmetric  finite  amplitude  sound  beams.  Indeed, 
nonsymmetric  situations  have  not  been  extensively  explored  in  the 
linear  case.  It  was  shown  that  because  azimuth  dependent  signals 
must  vanish  along  the  beam  axis,  the  interaction  with  axisym¬ 
metric  effects  cannot  affect  the  signal  on-axis. 


C.  Dual  Frequency  Axisymmetric  Excitation 

A  major  generalization  was  achieved  in  analyses  of  the 
propagation  of  a  signal  generated  by  axisymmetric  excitation  at 
two  arbitrary  frequencies.  In  the  limit  as  the  difference  be¬ 
tween  these  primary  frequencies  decreases,  one  obtains  a 
parametric  array.  The  analysis  followed  the  line  of  investiga¬ 
tion  originally  developed  for  the  single  frequency  case.  The 
initial  study  of  harmonic  formation  [5]  identified  the  mechanism 
by  which  the  primary  signals  interact,  but  it  was  limited  to  very 
short  ranges.  The  innovative  aspects  of  the  subsequent  research 
[13,  15]  was  in  the  identification  of  the  approrpriate  set  of 
coordinate  transformations.  It  is  such  properties  that  represent 
the  actual  interaction  effects.  It  was  shown  that  the  distortion 
of  each  primary  signal  is  dependent  in  equal  part  on  the  signal 
in  both  primaries. 

The  analytical  results  were  prohibitive  for  extensive 
numerical  evaluation,  particularly  at  long  ranges.  For  this 
reason,  an  interface  with  a  spherical  propagation  model  was 
developed.  The  same  type  of  interface  had  been  explored  earlier 
for  monochromatic  waves  [4] ,  but  had  been  abandoned  when  the 
Fourier  series  form  of  that  signal  was  identified.  The  remark¬ 
able  aspect  of  the  dual  frequency  study  was  that  it  showed  far 
better  agreement  than  earlier  theories  with  a  variety  of  prior 
experiments  on  parametric  arrays.  It  even  reproduced  features  in 
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the  difference  frequency  signal  that  had  earlier  seemed  to  be 
anomolous . 


D.  Validation 


An  initial  series  of  experiments  was  performed  by  Dr.  Mark 
Moffett  in  July  1982  at  the  NUSC  facility  at  Newport,  Rhode 
Island  in  order  to  obtain  data  for  a  comparison  with  the 
theoretical  model  [1].  Agreement  between  theory  and  experiment 
for  amplitudes  of  the  harmonics  was  reasonably  good.  The 
hydrophones  employed  for  measurement  showed  substantial  ir¬ 
regularity  in  their  frequency  response  and  no  phase  calibration 
was  performed.  The  discrepancies  between  theory  and  experiment 
were  shown  to  be  less  than  the  uncertainty  in  the  response  of  the 
transducer . 

In  September  1983,  Moffett  endeavored  to  improve  his  measure 
ments  by  resurrecting  the  special  purpose  transducers  that  were 
utilized  by  Browning  and  Mellen.  The  data  was  sent  to  Georgia 
Tech  in  its  original  digital  form  in  order  to  provide  a  complete 
data  base  for  comparison.  There  was  a  great  deal  of  difficulty 
in  reading  the  tapes  due  to  limitations  of  the  available 
equipment.  Eventually,  the  data  were  displayed  as  waveforms  and 
analyzed  for  frequency  content. 

As  the  data  was  being  analyzed,  it  became  apparent  that  the 
receiving  hydrophone  exhibited  extremely  anamolous  behavior  at 
high  sound  pressure  levels,  such  as  third  harmonic  levels  that 
execeeded  the  second,  and  fourth  harmonics  that  exceeded  the 
third.  The  causes  of  this  behavior  that  have  not  yet  been 
identified.  However,  an  extremely  important  verification  was 
obtained  at  a  far  field  location.  The  signal  there  had  decayed 
to  a  level  that  seemed  to  be  within  the  tolerance  of  the 
receiver,  even  though  the  source  was  being  driven  to  a  high 
amplitude  —  it  was  found  that  the  measured  and  predicted  signals 
coincided  at  that  location  [8]. 


E.  Other  Viewpoints 

A  key  aspect  of  the  King  integral  for  linear  radiation  is 
that  it  treats  the  signal  as  a  superposition  of  a  continuous 
spectrum  of  modes.  The  physical  significance  of  the  nonlinear 
interaction  of  these  modes  is  obscured  by  the  lack  of  a  discrete 
mode  that  could  be  traced  through  the  anaa lysis.  The  continuous 
spectrum  of  modes  in  the  King  integral  arises  because  the  baffled 
system  is  essentially  a  circular  wave  guide  of  infinite  diameter. 

A  finite  diameter  waveguide  has  a  discrete  spectrum  of 
modes.  By  exciting  only  one  such  mode  in  a  linear  sense,  it 
would  be  possible  to  follow  closely  the  manner  in  which  resonant 
interactions  take  place  nonlinearly.  Analyses  of  this  problem  in 
circular  [3]  and  two-dimensional  waveguides  [9,  10]  followed 
steps  that  were  suggested  by  the  analysis  of  sound  beams.  In  the 


circular  geometry  the  modes  excited  by  the  nonlinear ities  were 
similar  to  the  directly  excited  one.  As  a  result,  all  of  the 
harmonics  were  in  phase.  This  is  manifested  by  identical  types 
of  distrortion  for  the  compression  and  rarefaction  phases  of  a 
signal . 

The  two-dimensional  study  disclosed  the  existence  of  an 
internal  resonance  that  had  not  been  identified  in  earlier 
investigations.  Specifically,  it  was  found  that  the  transverse 
variation  nonlinearly  excites  the  planar  mode,  as  well  as  second 
and  higher  spatial  harmonics.  At  high  frequencies,  the  phase 
speeds  of  the  various  modes  coalesce,  which  results  in  nonlinear 
dispersive  interaction  between  the  modes.  Significantly,  the 
signal  in  this  case  was  shown  to  be  display  waveform  distortion 
of  the  type  observed  in  sound  beams.  This  verifies  the  nonlinear 
King  integral  approach  to  sound  beams,  which  treats  the  signal  as 
the  dispersive  interaction  of  neighboring  modes  in  an  infinite 
waveguide . 

F.  Parameter  Studies 

The  main  computer  program  for  sound  beams  was  modified  to 
provide  predictive  capabilities  for  an  assortment  of  transducer 
vibration  patterns  f{R).  Propagation  curves  showing  the  depend¬ 
ence  of  the  amplitudes  and  phase  angles  in  the  region  from  the 
transducer  to  the  far  field  have  been  carried  out  [8,  11,  14]  for 
a  piston  transducer,  in  which  the  particle  velocity  across  the 
piston  face  is  uniform.  Results  for  a  hypothetical  transducer, 
in  which  the  pressure  is  uniform  across  the  face,  were  also 
obtained,  because  that  had  been  the  basis  for  studies  using  the 
modified  Burgers'  equation.  Another  configuration  receiving 
consideration  was  the  elastic  membrane,  in  which  case  the  pattern 
for  the  transducer  vibration  is  a  Bessel  function. 

The  effect  of  diffraction  decreases  with  the  progression 
from  the  piston  to  the  membrane  configuration  because  spatial 
transitions  are  less  severe.  It  was  found  that,  although  the 
fundamental  frequency  signal  inside  the  Rayleigh  length  is  sub¬ 
stantially  different  between  the  cases  of  uniform  particle 
velocity  and  pressure,  the  second  and  higher  harmonics  were  quite 
close.  In  the  membrane  case,  the  rapid  spatial  variations  in  the 
fundamental  were  substantially  reduced,  and  the  higher  harmonics 
varied  quite  smoothly.  This  could  prove  to  be  useful  for  ap¬ 
plications  requiring  near  field  measurements. 


A  noteworthy  aspect  of  this  study  is  the  fact 
been  no  prior  extensive  evaluations  of  phase  s 
higher  harmonics  in  the  near  field  and  the  trans 
The  evaluations  revealed  an  interesting  trend 
angles.  The  higher  harmonics  are  close  in  phase  to 
tal  near  the  source.  As  the  wave  propagates,  each 
to  lag  further  behind  its  predecessor,  until  there 
difference  in  the  far  field.  This  is  significant 
spherical  wave  also  undergoes  a  90*,  phase  shift  in 
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from  the  near  field  to  the  far  field.  The  implication  of  this 
observation  is  that  asymmetrical  distortion  results  from  spheri¬ 
cal  transitions  for  higher  harmonics  which  are  delayed  by  the 
higher  frequencies. 
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V4.  Finite-amplitude  waveforms  produced  by  a  circular  piston  projector. 
Mark  B.  Molfelt  INaval  Underwater  Systems  Center,  New  London,  CT 
063201  and  Jerry  H.  Ginsberg  (Georgia  Institute  of  Technology,  Atlanta, 
GA  303321 

Measurements  were  made  of  the  waveforms  produced  at  six  different 
locations  on  the  axis  of  a  0.5 1-m-diam  projector  driven  at  60  kHz  m  the 
NUSC/Newport  large  acoustic  tank  facility.  The  locations  were  at  the 
last  three  pressure  maxima  in  the  neaiiield,  a  quasifarfield  point  at  5  m, 
and  two  farheld  positions  at  10  m  and  15  m.  The  projector  was  driven  at 
several  levels,  and  yielded  waveforms  ranging  from  sinusoidal  at  the  low¬ 
est  levels  and  shortest  ranges  to  shock  formation  at  the  highest  levels  and 
longer  ranges.  Two  different  hydrophones  were  used,  but  neither  had  a 
flat  enough  response  to  avoid  ringing  when  shocks  were  present.  The 
waveforms  exhibit  the  asymmetry  (sharp  pressure  peaks  and  rounded 
pressure  troughs)  previously  observed  by  Browning  and  Mellen  [J. 
Acoust.  Soc.  44,  644-646  (1968)]  and  predicted  by  recent  work  of  Gins¬ 
berg  which  accounts  for  diffraction  as  well  as  nonlinear  propagation  ef¬ 
fects.  [Work  supported  by  ONR  Code  425  UA,] 
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I  would  like  to  describe  some  work  Mark  Moffett  and  I  have 
been  doing  with  finite- amplitude  waveforms  in  cases  where 
diffraction  occurs.  No  doubt  you  are  all  aware  of  what  happens 
to  one- dimensional  finite  amplitude  waves,  such  as  plane  waves  or 
spherically-spreading  waves.  In  such  cases,  the  pressure  peaks 
travel  slightly  faster  than  the  troughs,  and  it  is  an  easy  matter 
to  predict  the  waveform  via  weak-shock  theory.  For  example,  a 
waveform  which  is  initially  sinusoidal  distorts  in  an  antisymmetric 
way  and  eventually  can  become  a  sawtooth  shape  because  the  peaks 
move  as  far  ahead  as  the  troughs  lag  behind.  A  much  more  difficult 
problem  is  the  prediction  of  the  waveforms  resulting  when  diffraction 
is  present,  as  in  the  nearfield  of  a  piston  projector.  The 
propagation  is  not  one- dimensional.  It  is  no  longer  a  simple  matter 
to  follow  a  pressure  signal  as  it  propagates  from  one  point  in  space 
to  another,  because  the  diffracted  field  at  any  point  results  from 
contributions  from  several  source  regions. 

Diffracted  fields  are  normally  calculated  via  linear  theory, 
and  so  nonlinear  effects  like  distortion  can't  be  handled  readily. 

We  have  been  working  at  Georgia  Tech  on  a  new  approach  which  can 
accoxint  for  diffraction  and  nonlinear  distortion  simultaneously, 
under  an  ONR  contract.  We  needed  some  experimental  data  for 
comparison  with  the  theory  for  the  case  of  a  circular  piston 
projector.  Mark  Moffett  was  asked  to  see  if  he  could  provide  some 
data.  I  will  show  you  some  of  those  results  shortly. 


<<  Viewgraph  1  -  Browning  and  Mellen  Waveforms  >> 

First,  I'd  like  to  show  the  kind  of  thing  that  happens  when 
diffraction  and  distortion  operate  together.  These  pictures  are 
from  a  1962  letter  to  the  editor  by  Browning  and  Mellen.  They  used 
a  tiny,  8/lOths-of-a-millimeter-diameter  microprobe  as  a  hydrophone 
and  looked  at  the  waveforms  on  the  axis  of  a  40-centimeter  square 
projector  driven  at  150  kHz.  The  top  left  picture  shows  the 
waveforms  measured  at  three  different  levels  at  the  longest  range, 
which  was  8  meters.  Then  the  next  two  photos  down  the  left  hand 
side  and  the  three  down  the  right  hand  side  show  the  progressive 
distortion  of  the  pressure  waveform  as  the  range  was  increased  from 
1  meter  to  5  meters. 

Positive  pressure  is  up  in  these  pictures,  and  so  you  can  see 
that  the  pressure  peaks  are  sharp  while  the  troughts  become  rounded. 
In  other  words,  we  are  not  evolving  toward  the  sawtooth  shape  which 
would  be  expected  for  plane  or  spherical  waves.  Diffraction  and 
geometric  dispersion  shift  the  phase  of  the  harmonics  relative  to 
the  fundamental. 
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Fig.  S.  Wavafoaa  at  1  m  (Ft  «O.S  bar). 


Fic.  9.  i^'av^wm  at  4  m  iFi>4).S  bar). 


Fic.  b  H'atWimi  al  2  m  (Ft^O.3  bar). 


Fig.  f.  Wat^orrn  mi  S  m  (P#*0.5  bar). 


•o  u  to  minimUc  wall  rellcctiona  For  (he  profKtor  dimensiont 
and  wavelength  of  this  experiment,  (he  Fresnel  tone  extends  to 
about  4  m  but  the  Fraunhofer  zone  spreading  is  quite  insigniheant 
over  the  remaining  tank  length.’  The  hydrophone  is  (hen  set  on 
the  projector  axis,  and  waveform  nhaervations  arc  made  at 
appropriate  disLancca. 


Figure  3  gives  a  oomparison  of  two  wavefotma;  the  more 
sinusoidal  one  xras  measured  at  2  m  and  the  more  distorted,  at 
8  m.  The  initial  peah  pressure  Pt  measured  at  I  m  was  ap|>roai- 
mately  I  bar.  The  oacllations  on  the  trailing  side  of  the  (race 
are  art  tets  caused  by  "shock"  e.xatau<in  of  h>xiraphone  reso¬ 
nances  that  (all  ahove  (be  useful  frequency  ranj[e. 
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<<  Viewgraph  2  -  Block  Diagram  >> 

Unfortxmately ,  we  couldn't  use  the  Browning  and  Mellen  data 
for  comparison  with  the  theory,  because  the  theory  hasn't  been 
worked  out  yet  for  a  square  projector,  which  is  a  three-dimensional 
problem.  We  therefore  decided  to  repeat  the  experiment  using  a 
circular  projector.  The  projector  was  a  20- inch- diameter  array  of 
60-kHz  tonepiltz  elements  made  by  Raytheon  Sub  Signal  Division. 

The  projector  was  driven  with  60-kHz  pulses  which  were  amplified 
with  one  of  Bill  Konrad's  20-kilowatt  drivers.  The  hydrophone 
output  was  captured  with  a  Biomation  transient  recorder,  which  is 
actually  a  digital  machine  containing  an  A-to-D  converter.  After 
capturing  and  storing  the  waveform,  it  was  plotted  on  an  X-Y  recorder. 
The  experiment  was  done  at  NUSC/s  large  tank  facility  at  Newport, 

Rhode  Island. 

It  would  have  been  nice  if  we  could  have  used  one  of  the 
Browning  and  Mellen  microprobe  hydrophones,  but  we  would  have  had 
to  operate  the  projector  just  below  the  surface.  Also,  we  weren't 
sure  any  of  the  probes  were  still  working.  We  tried  three  different 
hydrophones,  hoping  that  one  of  them  would  have  a  flat  enough  response 
to  accurately  reproduce  the  waveform.  The  first  was  an  LC-5  hydrophone 
made  by  Celesco,  or  what  used  to  be  Atlantic  Research.  The  LC-5  is  a 
l/16th-inch  cylinder.  It  is  just  about  the  smallest  commercially- 
available  hydrophone  there  is.  The  second  was  lent  to  us  by  Gerry 
Harris  of  the  Bureau  of  Radiological  Health.  It  consisted  of  a 
polymer  membrane  stretched  over  a  metal  hoop.  This  hydrophone 
turned  out  not  to  be  useful  below  a  megahertz,  apparently,  because 
of  resonances  involving  the  supporting  hoop. 


X-Y 

recordor 


<<  Viewgraph  3  -  Raytheon  Hydrophone  Sensitivity  >> 

The  third  hydrophone  we  had  was  also  a  polymer- type  used  in 
medical  ultrasonics,  but  in  this  case  the  polymer  is  backed  wi^sh 
a  silicone  rubber  absorber.  This  hydrophone  was  made  by  the 
Raytheon  Research  Division  and  lent  to  us  by  Roger  Tancrell  and 
Dave  Wilson.  This  is  a  plot  of  its  sensitivity  and  you  can  see 
that  it's  not  bad,  but  it's  not  really  flat  either.  Nevertheless 
the  Raytheon  hydrophone  turned  out  to  be  the  best  of  the  three. 
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<<  Viewgraph  4  -  Waveforms  at  2.6  3  m» 

Here  are  some  typical  results.  They  are  waveforms  on  the 
projector  axis  at  5  meters.  The  waveforms  correspond  to 
approximately  10- dB  increments  in  drive.  The  receiving  amplifier 
was  also  changed  in  10- dB  steps  to  make  the  signal  levels 
comparable.  Positive  pressure  is  up  in  these  plots. 

You  can  see  that  the  waveform  goes  from  nearly  sinusoidal 
at  the  lowest  level  (shown  in  black)  to  a  form  with  sharp  peaks 
and  rounded  troughs  (shown  in  red) . 

You  can  also  see  that  the  zero  crossings  on  the  ascending 
part  of  the  wave  shifts  back  as  the  level  increases.  In  contrast 
zero  crossings  are  unchanged  in  one- dimensional  waves.  The  rise 
in  the  peak  pressure  which  accompanies  the  narrowing  of  the 
compression  phase  is  consistent  with  conservation  of  momentum, 
as  is  the  broader  and  shallower  form  of  the  rarefaction  phase. 


<<  Viewgraph  5  -  Waveforms  At  All  Ranges  >> 

This  shows  how  the  waveform  evolves  with  increasing  range 
from  half  a  meter  to  a  little  over  15  meters  from  the  projector. 
All  measurements  were  made  on  the  projector  axis.  The  upper 
three  waveforms  were  taken  at  axial  maxima.  (Probing  was  done 
to  make  sure  that  they  were  maxima) .  The  5-meter  point  is 
beyond  the  last  axial  maximum  at  2.63  meters,  but  is  not  yet 
in  the  true  farfield.  The  10  and  15  meter  distances  do  qualify 
as  genuine  farfield  measurements. 

These  plots  are  for  a  source  level  of  approximately  243  dB 
with  respect  to  one  micropascal-meter.  You  can  again  see  the 
sharpening  of  the  peaks  and  the  rounding  of  the  troughs  as  the 
propagation  distance  is  increased.  Each  plot  begins  at  an 
instant  where  the  pressure  is  zero.  You  can  see  that  the  first 
zero  crossing  moves  back  in  retarded  time  as  the  range  increases. 

I  should  explain  that  the  zero  of  the  retarded  time  was  determined 
by  referring  to  the  signal  at  a  very  low  drive  level,  where  it  was 
sinusoidal. 

In  order  to  be  sure  that  the  nonlinear  distortion  we  measured 
was  not  due  to  hydrophone  nonlinearity,  we  determined  the  second 
order  sentitivity  in  a  subsequent  measurement.  It  was  low  enough 
that  we  didn't  have  to  worry  about  it.  Also,  we  can  see  that  the 
waveform  close  to  the  projector  is  nearly  sinusoidal.  Since  the 
highest  levels  were  measured  at  short  ranges,  the  short-range 
waveforms  would  be  more  distorted  than  those  farther  away  if  the 
hydrophone  were  behaving  in  a  nonlinear  way. 
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Viewgraph  7 -Computed  Waveform  >> 


Here  is  a  comparison  of  the  computed  waveform  (shown  in 
black)  and  the  measured  one  (shown  in  red)  at  2.63  meters.  The 
irregularities  in  the  analytical  result  arise  from  nonconvergence , 
but  the  results  are  reliable  at  most  time  instants, I  should 
mention  that  this  is  pretty  much  the  worst  case.  The  other 
computed  waveforms  were  much  more  regular. 

When  we  compare  this  result  to  the  measured  waveform  we  note 
some  significant  discrepancies.  However,  we  must  also  accotmt 
for  the  hydrophone  response.  We  did  not  have  a  phase  calibration 
for  the  Raytheon.  Instead,  we  wrote  a  microcomputer  program  that 
describes  the  analytical  result  in  terms  of  its  frequency  response. 
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<<  Viewgraph  8  -  Frequency  Response  >> 


A  Fourier  analysis  of  the  computed  waveform  was  used  to 
reconstruct  the  signal  from  the  lowest  five  harmonics.  This  is 
the  green  curve.  Using  the  known  amplitude  sensitivity  and 
modifying  the  relative  phases  of  the  harmonics  yielded  the 


result  shown  in  black.  For  comparison,  the  experimental  result 
is  shown  in  red.  The  agreement  with  the  measured  signal  is  remarkably 


good. 

Similar  analyses  for  the  other  cases  indicated  that  the 
theory  over-predicts  the  distortion  very  close  to  the  source, 
but  it  seems  to  be  quite  accurate  beyond  one  transducer 
diameter. 
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LLS.  An  improved  King  integral  algorithm  for  finite  amplitude  aound 
bcama.  Jerry  H.  Ginsberg  (School  of  Mechanical  Engineering.  Georgia 
tnsiiluie  of  Technology,  Atlanta,  GA  303321 

This  paper  supercedes  the  previous  analysis  (J  H.  Ginsberg,  J. 
Acoust  Soc  Am.  Suppl.  I  71.  $30(19821]  of  the  infinite  baffle  problem  for 
an  asisymmeinc  harmonic  excitation,  which  denved  a  nonlinear  King 
integral  describing  the  distoaion  associated  with  finite  acoustic  Mach 
numbers.  That  analysts  was  shown  [M.  B.  MolTctl  and  J.  H.  Ginsberg,  J. 
Acoust.  Soc.  Am.  Suppl.  I  72.  S40  (19821]  lo  exhibit  excessive  nearfield 
distortion  m  companson  lo  experiment.  Using  the  linear  King  integral  in 
ns  conventional  complex  function  form,  as  opposed  to  the  real  function 
analysis  employed  previously,  leads  to  formulation  of  the  second  order 
poleniial  in  terms  of  complex  funclions.  Asymptotic  inicgration  of  this 
potential  function  reveals  that  there  may  be  significant  contnbution  from 
the  evanescent  spectrum  (small  transverse  wavenumbers|,  as  well  as  from 
the  propagating  spectrum.  A  coordinate  straining  transformation  de¬ 
scribing  ihe  full  spc.'irum  is  deduced.  From  it,  the  previous  analysis  is 
shown  to  be  only  asymptotically  correct.  The  new  analysis  reveals  that  Ihe 
distortion  is  governed  by  a  transformation  that  involves  Fresnel  integrals 
for  the  propagating  spectrum  and  the  error  function  for  ihe  evanescent 

spectrum.  Some  comparisons  of  Ihe  analytical  prediction  and  the  results 
of  Mofleti  and  Ginsberg  are  presented  {Work  supponed  by  ONR,  Code 
420.) 


Stouffer’s  Cincinnati  Towers 
Cincinnati,  Ohio 
9-13  May  1983 


/Ix'Sy  KVI  'T  (C  5o  R  a  of  I CK  j- 1  Oy) 

« vSt^  -T^  ” 

J~~ I  1^  I  "I ^  Xl  ✓v?  pli-fuc/c  E'TRe.c'^s 


5/ 

COry)prf^S(\}ie  ^ ^ : 
IriK/iScid  ^  I 


syr<ni^<-i^y  y  I !  I  C‘\{  COJ^!^ 

-R-^ 

*1  2  ^ 

t 

^77>>77V  u  / 

ri  ^  r  C  /  c  ^  ><  c  /  T  ^ 

r€  -5  3(y(re  or  loci  i-y 


Ins/iScid 
RliA  I  d 


E  1 1 


'  S  S C/ 


'/I  f  -c  r7 


-e  X  «  ^ 


=  <£  -fCR;  .' 


w/n  ^  /^€ 


-<  -«c 


/  w  f  ^  I  r  I  i  e  , 


-z- 


I.  J,C.  Lack^ooA^  T.G,  Mu.r^  P.X  2>{^c\.sUci<  0^73;  g:; 
a  .  F<^  ^  fiel4  at.>n€x  I  ^  ^  1 5  \'v- 

b.  i  -e  6i.  i -e  ol  <s  s  ^  f  I*,., 

s  /?  /o  ^  (  c  ^  I  \Ayci  v/  €  ;  /.  ‘ 


F  .A.  F:a  lo  o  i  oi'ska  y<l  ^  “iR,  \/-  e  (' 

G  .  hi  p  ^<r  Fi  ^  (A  ci  I Y  sf  s  -  no'^  -  >'. 

ploi^Ci/^  yy^yet  \y  ^  S 

b.  Bu  iTp  )or\  nn  oMt  ft  lay  ^-. 

a<dJi^^  rfL  Lpy>lp.C(ai^  oy>-e -f- o  ^ 

for  ^  4" ^ ^  Airecifoy 

ff-t  e  ^c  /■n  sy?<!ic^ 

Q  r\d  ftioat,. 


.  %'  %*■ 


Po  I-  tG.f 


p O u  10  Aa  "Tf  0^  5 

/?  lO  «  ly  5 1  5 

Poc^cf 

-  cM  V'^  =  c,  ii 


-^oCo'  jf  -t  0(4^") 

>0^^  (Orf  Tecfoi^i^ui^b 

(  |)  Qi/^^c.'f  ryit-flooA  iASi>^^  per'  ioQp 
3<ri€S  -  pr  B'^f  OU^Iy  e  r>opi  oy^^  A  Por 
/)0^  '  of  t  fo  r  ry\  Pa,  r- 4 ^  ^  ^  eyh^Ar'f  C^l 

a^A  Sphe/^tcaf  sAytK'res 
M  J\lo^(ir\^ar  xy^avt  e.<^U(a  ’ftp^ 

v’-({,-^^_  (y.,;^  v'^<i> 

-K  (y<i>-'7'4>) 


’A  v:"  >7^  K' t 


ITr-».U.l»l  HllJ»^.yii.r^  ,..in  ., 


-4  - 


I  » 


I  ( 


C>  )  ^>  €  >^  I  ^  ^ 

(c)  of  r-€  n  o^x^^ohl  a-\  Ipy]  ~ 

Reji^/re  S  Q  ( ^') 

(2  j  yT'i  -  USC'^  4o 

cl?,  (  •^  80) 

(3  J  J-f<^f/Ofiary  plna^c  (y\i  e  p  roji  1 0^  -  uS 

-/o  d^4-e  '  ^  t  in  €  fA« 

^1  f icci  n4  pci^4  <f>^  (f9Sl) 


»  i 


r.  ..i( 


i . 


-  ■  ■  ■  1 

••:-■  \  .•> 

,  "v  ^  * 


i 


I 


e  5 


1-1  ^  (i  e  (  ira^sCo^^  { cz  y  / 1  ^  i  c  a  I  co  0 ai i  ) 

4^,  ^  $,  k  'Jq  C  k  R)  k 

~rt- Oi-^S  ^ 0  ^  ^  o(.e.  ^  ^  .  C. 

'^l  =  /  V,  +  ©,  )^  (AR)^A 


■}-  r-^  *1  5  Co  r 
Y^loci^/ 
ft  p/ (  f  K </< 


"tfe*  n  5  fo  f  (/ 
phase 
<3 le. 


Re  2^H 

-  / 

4>,  --  /  (<^05  %)jJkH)^l< 


e  e  c 


> 


[U  ^  i-  ~^.  Q  -  Ti.  T  CK1<  \<H  I  pUctSZ  \/cK  ,^(0  L  (e 


(kR  )  IS 

0  5C.I 

UgI-O  ^ 

fR=^  -iH  //^<a 5  i/f 

VA/a  ve 

/6 

1/  k  . 

1  In  €.  Spo  n  S  <. 

c  onsisH 

of 

Gc  (LO'^  1' /  •n  G  0  u S 

sp  eci K'tyf  KV1 

-4  rG  tns  e 

\yyc^  (/  €  /  ^  H  s  . 

L.  G.  ^  ^  s  V-C'^  5  ^  SA/^  /  ^s(^  /<■<  1  ) 

G  4  € 


p'^^p 


.A  / 


'  p  e  ^  I  K/  e  I  y' 


ly  -  1  ^  I  5 

k- 


He  p  i^G  s  e  sp-e  e  M ^ 

S  1 1  s  ve^se  \AyG  U^ip s  (k>l) 

—  /ip  I  laic  fo^ 


■e  \/<^>^  e  6  c  e 


•7- 


Scc<Oi^M.  T-e  ^  5 

1 .  pri»y\a  I  :  Ij  hfy  /  5  t  ^  ^ 

wk^cf^  c  ^  i/yid\ie^  c  ra  ib  /  C. 

o  f  ^  <t><  i  . 

2.  ^/5e  /o  ^Q/-^  50UrC‘C  'l-e>^  5  fo^  ^ 

Kj^^/^c-f'  ^  K/^  e  s  c  ^ 


•  ^ 

•■'  ^>V'‘ 

«  "  »  •  k 


(s>  e  ^  -p 


o  ^  ' 


^ ■  "aT-^  =  //t  Jo  >  Jo 


o  <t; 


4  -T^  ni^^m  A/c 


E-^r/lfr-  }y\/Qirk  /'I  C  Cf  f' \// h  e*i  ^  C  OOr- J  i  ^  S  ^ 

Q .  (p  ^  ■f’Ae  doo^<^ft^  \^l^e>^€, 

Sy  *^P  j-Q  k  fC  <Z/:p^^SfO'^s  gf  BtsSel 

\{yf^  c  I  5  «t  valid, 

b.  Ex^'Tlt-^C  -kke,  ca  iyt^p  /  c?  f  /  c 

s  o  I  U  'k  i  O  ^  <P  ^  ^  ^ 

I  -f"  -I^O  id  €  d  L4  c  e  -f^e.  f^O(r»^ 

P<f  <=t  ■/ Uf  3  3-g  g  5  €^  / 

•f  u  «^C  f  *  0  S. 


•-:--V:-rj 
u  -J 


i  « 


_  L _ I 


y.«'. 


5  r  A  c  ot- 

50u«ree 


(1>^  =  /''!/•', 
o 

■^[X n  j/(r)Ry«j  k  j/as) 


(y^s  e  A.  S  yt^  fi  fo  j- !  c  ■€.  )(  p  S  i  Qyy  i  k  1^  ^  R  >^  L 

V^  ^  fc(k^A)siy\(  %±  %]  cos  [(kR*-"/^) 

t  (k>  -  tr/^  )]  ^  ^ 

S  €.  Q  (A  (  O  !  '(‘Y  y\  -^  \<.  iscpf^ 

ca  /  v/e.  ^ 

0.  Lei  y\  =■  k'i-  ^  ^  ^  }  4c^l  1 


C.  P  (  ^  f  ^•^e.'nc  C  Siph5  '=^  / /q /e^  /5 

C  5  I /yj  ^  C'f  f0<^  op  ^  ly 


4  o  i  (  c 


Re  I  ut 1 4  : 


/-I  k  /^  >  >  I 

o  '^k 

-t  Co  5  2  Clo5  (z  ^  '  ^^2.  ) 


5(^2 


A I  />  p  ‘ 

^er\e'ro^\  Solcil'iO'n  fo^ 

H-[^  <4n] 

/ipo  ^^4  ill  1 5  -To^  n-t  k/?> 
Ded^cz  ^  (ijofk 

i-^clt  p)  ^ 'id.^ 'll  of  R)  Si-/c4  i' 4 

C\  %y  yyy  tO  \' O  !  I  C-  f  ^  ^ ^  CP  ^  ^  ^  C  i^Ci  { 

s  ol  iA  4  /  <p.n  iviatfcAeS  -/'^C 

o^^Y^ploii^  ^csul  4, 


'T1 


I  f  »  ^  P  Po  f<«/  t  lA  ^  c  f  n 


5“  \/<?t  /  FC^) 


.  C.O'^5€ ^  ^  '  O C>  ^  ^yi  ^  ^  5 


h.  /)  y  /  ^  I  v^«-  ioc  (!'  y  £  t 

c<.  pl<^y>€  2"  =  (^(9 5 ^  ■pc?'^ 

^  ^  t  o 

c.  t/^  /u  €-  yt^U'^i'  \ya  f  ^  L 


cl  .  Co>^  C  I  [A^ 


h(2) a  o 


c^  ^  lo  !  'Ct,  ritSuHs. 


^--  ^/  ^  {(^<=5  n  )  j;()t/^) 

^  cos  ifjLj-.fU)'-  j/()=0']}'^^ 
4  £  ( 4’,')«v  ■*■  0  ^g"')  <«  Z''  I?  ^  H 


Piiy‘5/c£il  ^e5^c7^5e 


sjCAocdy  \ 


N/^  -  e« 


-eJ> 

3  ?<’ 


=  c, 


9t/> 


p- 


p  - 


o  4  €.  ir&  O'^  f>S  v/'*  I  e^  lo  ( e  ^  Lxi 


G>  f  yr  e  ^  ^  I  f 


o  ^ 


0* 


,  -  f  [t  >A-Pk,  +) 

+£’-  R^i)  ]  ‘^  l‘ 

-f(uj  >  -/-  0^^^) 


5/^.0/  p/<"<l#lC€  op  SCC0*^<^ 


-/“e 


5 


^  S  2- 


/^/<?^  -  n  Acc  of  ^  ^Cy 


\AJ  h  !  "I  i/\  <Q>^1 1  ^  t  •!  y 

I .  \^  \  <n  e  ^  ^o  ^  cC  t> 

^  SSe^h  ^ If y  QO  ^ ^ 


2,  Loc^^{Oy7  Q  o ^  ^fO  o^  *t  t  ^  ^  /t? 

sp^ct  Cc  pUoS^  /5  sAz-P/d**/ 


a( oi  €.  {^o  ^  o^h  ^  f  fy 


ct  fp-^^  A  iP>n  /A€, 


SAk^  Ot  \A  ^  4^  t/  /VT  fo-^r 


e  ^  0/^yyict  1 1  2  )  on 


\n[Jlept  ^  \/Cin  t  Cc  IqI  ^  ^ 

L  m  r 

t^On  line 

^k 

kR 

A 

5  c  n'0\^  (X  t>  "B  .  7~*^yr 

•<^  T-  £  z'^*-  Fj  ,^5*,^) 

kR  --  ^  fi,Ck,>,/,^^^,t) 

C  A  £?  o  5  -e  f"?.  ^  ^  g  /VT  q  II  OCc-^) 

•fen  3  fl^^f  ^  n  e.  m  o  f  w*o>  fo  nynly 


c:»CC(/<*^<^fC 


Li  ^  of^J  <  /*?  g . 

.  I 

u,  i  J  {s 

o 

no  f 


I  l)£.  Cl  /^C-  no  ^  Afrd  Tf^  ut^^fcA 

ca^ce  {  cell  0(^C^‘S '''*')  =#> 

—  os  <L(llct-l’0''y  «/-e^*-kki  /5  iT  €  icf  e 

<x^//4  Jn^/'^.SLS/^^  2- 

■f^t5  '1^/*^  ,5 

R-e.oi^O'n:  (Lo  ^ t  ^  cf't*^  s~l  ^  O/  *7/  »T  ^ 

^  >»0i^_  2^*.^ 

V<^ltA^  y^Uick  15  C«=?  no  €> /g  g/  L/ 

4  <<y\  a  I  r. 


I,  d  a:?  I 4  f 


^^  2  -  cl  y-^i)( ^  )  ^kV^  s/n(-{  ^  e^ 
kR'yS^  ^  ^  <2os[{  ■*■  e^ 

-  '<,  -  f)  J'X  f/^,) 

ressc/'^  c 

O  '^'|i  ^ 

*iLr.,)c(f^rk V,  [J^  U j; (^,fj}  J k 


L 


\A/Ci  \/e  u  ^  k 
CX.  o^  -  C/f  »  Cl  p  /g» 


^  I  V  ^  s 

K^ya  V  e  le-l  r£-  S€  Icl  I 


o<y\  e. 


<2.  / /^dU  I  <5 CLy  1 1  y\  d  C^  I  uLya\/f^ 


/  y\ 

u  iJ  e 


'Z  .  S  c  l\  \A/<^  \/€  /  ■€  /  p  /^op<^i  p-  £x4  ^  5 

(d\  spe  i  \y  €.  ly  -fAe  <^xi‘=^( 

ol\  ^  c  {  i  o  yy 


3-  T  k  ^  g/  /  3  ■/•(::?  y''f'/ o^n  (i,e,  7"^^  eooy'^ 
i  i^ci  5-Po^^  ^4- f  o<y\)  o-f  ec(cC^ 

)  5  uy\^i'^’^^c'f-€ai  i^y 

'A/g?  Ic  i  S, 


A- ^  ~T  h  Q_0  O  r'ol )  y\^^C  ~l  y"<^  5 -f^o  yyyi  ^ 


R  e  f  C  •(■  ( O 


^  Ina  O  ^  /"  t s  e  1*^  ^  uu^<?  v  € 

^e.ipc>/\5i  AkiaI:  TV^*^sv<’^5e: 

[/aria^blc  -f  *<•  ^k.~ '^k  ^JL 


jDr^5$ure  I  p)^ 
a.  X  /Cil 

ve  loci  4y  :  (  J|, 

■f  r  5  V  ^i^s'e 
\/<l  I oc  i  4  y  :  (_\/^)^^ 


L  ) 

S  (  ) 


cos 


(  ) 


->0  1  J 


J<,(  ) 

j;(  ^ 


posi-l  (0^ 


Ayi^i  : 

k 


Axi^l:  Tr^  ^ 


•Sl-^  i  ) 


Tr  ^/>5  verse  T  Ic  R 


cos  C  ^ 


-  2.0  ' 


yf>i5 


N  L/<  1^  C  ^  I  C  q  (  \  r- a  i  !  o  ^ 

(<^)  D(  ^  C  r^i  !  1  ^  ^  '  ^o~  ^  j  ~  ^  ^  ^  ^  ^ 

TunTTrY  W-  ■tert'Or 

Q.  i  K 

(  c)  E  ^iA  c{  I  tp  a  c  I'^p.  >^ecesS£>^y) 

k,„  -~ki  =  h 

[  J)  !  e  r  po  (a  iiy(p  pofy  *no‘^  I  ‘q  I  ^  -  n  poi'^i 

{  Kj -e  iA/ 0 ^  Coin,  lae.) 

<2.rror  -  ^ 

0  (1^”  'J  I P  K»  '3  e^/f^ 

(  <2^  C  r-  /  i  €  ^ffeciiy^^  cl^O(C€ 

Noi  ^  -f  b  ^  /  r(A)  <o5cilUl^5  rapidly 

\A/lqe-A  2  Or  R  15  (<qrp  e . 

(l)  V\\yl^€r  a>i  ^  <r  ty/ <qh\/ i  i  ECi^) 

(z )  J)^  c  r  ^<si  e  lr\  a.  %  r  -  [  ?  '''-t  R  )  />qc  re-x^s 

(f)  C.-orc  I ci  5/0  r 

TrcpezOtdo^l  rf^le.  ^  <  iy^o 


^rror  - 


0{lo 


IA\C  re^Sf$ 


K- ^ 


2,  ,  hi  [A  I  C  ^  ^  ^ ^  ^  i  ^  '  N-  '/h 

P~y^oC„^e^  'vA/  I  ^I^[i -t 

'  t-o  ''  A|^ 


4  C„'  £  Z  W;  (  ^k 

L=^/ 

-f-  ^  ^p  (-  ^1, 

^rapiA  CO'^ 'Ae  t 

P"  k/O,  i  iA  £(  i- t  0  O  -T  CX>O^AtiACL\ C  i^Af^sCo^ 


p^  C  ( 

fy  h 

^  ( 

7,  0=^ 

1  ^ 

f  *vf 

o  yr  ~  Rg 

pf^  5  0  y^ 

tn/t  ti  od 

(a.) 

Gui^'S^ 

^Oo\r  e(  -  c<^' 

V 

’;  A 

(b) 

Tbu  < 

roO  f' 

15  a  4* 

k 

(c^ 

C  00^ 

U 

\ in^S^0r/y\  Q.  i 

/(?/1  ,' 

;>,  ^  = 

pf 

^v 

,fiy)  J  l<R  = 

^k) 

I  0<'^\Or  5 


(e)  Solve 


fa-i)"' 

(ij/" 

.(^  y 

1  A<'-^ 

5 l^oc  s 

uur 

’  (.7  ^"’ 

kR-  f(^'‘'A^'' 


ife  *"«  ^ 


(Hr  (If)' 


>T<»  c of"  -fr^<^sfc?/ 

aiiOn  a/- 


yj~a^  CO'^l  ^O  5A  0  c  it. 

(e.^.  i  ^O'l-  M-ii^we). 

0*^  -  COnveA'H^^c^  ^(o^s 


/5  ro^clytfff. 


-53- 


Co^A^4fO^ 


1.  P/*  £. 

S  5urc 

^y;  '' 

P(,-,. 

^  /^o  C 

£:  -S^ 

7.  ^ 

A  C  / 

4r  <kr\  5  P  < 

j^K»i  iri\/e.^  1  ( / 

Pl,-0 

o 

-1  yj^^Ak!^) 

3.  E^\/a  lu  a  f  < 

-e^  r  // 

%olu\tOr  0^  ?- 

^^-0  i  \c  R 

PU-0 

©^  )  \Jq( (<  R)  k  y- 

-t  /c<u'£  f°  ^  COS  (h  d^)jJkk)Jk 
'  \ 

4,Maicl^  \oOoi*^J(>y"y  uai  : 

|<  <  ,  ;  R&pUc^  ^  bv  4  ©1^  by 

k>,;  by  p,4 

5.  St  ^  ^  ^  I'/’y  '  -  vA/Af^  1^  ^  / 

ci^alyfic, 

V^/oc/Zy  b,c.  ;  j^ylc^rai^ds  arc, 

ir\-( ^ p  rA.\o  le , 


M 


]  = 

{^<R  <  ■^^ 

c^ 

( 

4' o 

(  0 

R  >^ 

or  C4.  pi^j-O'^  (^l,e.  velocify  la.c,) 

B  eci  ^  ■(o/'ioi  / toXco  >  I 

Par  ^  (  e  !  c(  :  r  ^  (  R^-f- i^)''''^  -  Q  ^  ^ 


^  ^  5  -P O/’-^  3  ; 

^  C-O  I  -r~  /  d  C«5 


p,-  "if  f 


Pk  ^  i  (^)"  -  ^  ° 

lr\pt4l  pa  f' d  f  I-  <rs  ; 

Specified  f^  {)  ^ 

^  ^  filepe^d  oriy  o^  ^  ^  e(rt<) 

W<  ^  :  T-/.4  <=“ 

lA/aif^  ;  _§.  s-  ^  S' 

■■  /t 

,  _  6'ti).,,  _  - 


k  ^  o 


^  i^a  /< 


-  a32  £•,,, 


Vg.  /o c  1  ly  u,  C. 

{J/\r  f  td  SiA^ui<^ir  t  I  y 


'  k  k 


kF  ^(y+'X^^)  '  kVj^  cor(i^6^-<^- 

K 


^  Wki€/n  k  ^1  O, 

(<^)  Bol^  0  ( t  )  -I  ^  U  r^  C  <si  (/  f, 

(t)  ^  ^Inock  <*»  ■/  5<^«?//  ^  =^ 

no-f  p^yst^^Hy  accepi^ble. 

5.  Noi  ^  '. 

(^)  '^/'^  pl^ik  S€.  ^ 

(y)  f  0(e  e''’-)^  i  ^)^"= 

4 .  l-ly  poJ  ^15 

'T~h€  c:^taO[/€  -f S  foryyiah  0^  »5  -f  A  ^ 

exp<!it^s/o^  fo^  {^r^e  ?  (?r  ^gs .s g / 

^  LAy\  ci  •0'^  ✓Vi 

71,^2  --  -f-  e^'Yi-OtV^ 

-  CO  5  (t  +  ©^)j^(-<,.)]  j;^A> 

kR  =  /3^  £(r<-i)k  v;  2[sM(^f  ©1^  J  ) 

■4  C  £35  (■^  ■'■  ©L  )  'J’o  V ✓k')]  vTo^  CA ) 


e-/o 


Revue  d'acoustique 


11®  CONGRES 
INTERNATIONAL 
D'ACOUSTIQUE 

11th  International 
Congress  on  Acoustics 

11  Internationaler 
Kongress  fiir  Akustik 


HORS  SERIE 


VOl.  1 


IV  ICA 


Paris  1983 
Lyon-Toulouse 


1  6  _ 


FINITE  AMPLITUDE  ACOUSTIC  WAVE  PROPAGATION  IN  A  CYLINDRICAL  WAVEGUIDE 


GINSBERG,  Jerry  H.  and  MIAO,  H.  C. 


Georgia  Institute  of  Technology 
School  of  Mechanical  Engineering 
Atlanta,  Georgia  30332 
USA 


Introduction 

A  signal  propagating  through  a  hard-walled  circular  duct  can  undergo 
significant  nonlinear  distortion,  even  if  the  acoustic  Mach  number  is  a 
small  fraction.  Keller  and  Mlllman  [1]  and  Nayfeh  [2]  investigated  this 
problem  under  the  assumption  that  other  duct  modes  are  not  resonantly 
excited  by  driven  modes.  The  present  work  discloses  that  there  is  a  non¬ 
linear  mechanism  in  which  cumulative  distortion  is  generated.  The  method 
for  analyzing  this  phenomenon  is  a  direct  perturbation  scheme  using 
coordinate  straining  transformations.  Much  of  the  development  draws  on 
techniques  used  by  the  first  author  to  analyze  cylindrically  propagating 
waves  [3,4].  For  simplicity,  only  axisymmetric  duct  modes  are  considered 
here. 


Analytical  Formulation 


Let  (z/k,  R/k)  denote  dimensional  axial  and  transverse  coordinates, 
respectively,  for  the  cylindrical  duct,  and  let  dimensional  time  be  t/Q, 
where  is  the  frequency  of  an  excitation  at  z  =  0  and  k  =  with  c^ 

being  the  speed  of  sound  in  the  reference  state. 


The  nonlinear  wave  equation  governing  the  potential  function  in 
this  case  is 


- 


=  2(B 


-  1)  +|-  (v<J)  •  V(}))  + 


where  6q  is  the  coefficient  of  nonlinearity  (”  1  +  B/2A  for  a  liquid), 
a  hard-walled  duct  the  potential  function  must  satisfy 


(1) 

Fo 


'R-ka 


9<l>| 

3r| 


R“ka 
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where  a  is  the  dimensional  radius  of  the  duct.  The  harmonic  excitation 
considered  here  gives  rise  to  only  one  duct  mode  in  a  linearized  analysis, 
which  is  the  case  when 
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where  denotes  the  Bessel  function  of  the  first  kind  of  order  zero.  Noti 
that  yka  >  0  can  be  any  of  the  zeroes  of  -Iq'* 

The  acoustic  Ilach  number  £  is  much  less  than  unity.  It  is  used  as 
the  perturbation  parameter  for  the  potential  function. 

4)  =  +  e^(j>2  +  ...  (4) 

The  leading  term  in  this  expansion  is 

4*2^  =  -  sin(t  -  az)  J^CyR)  (5) 

where  the  wave  number  a  is  found  by  satisfying  the  linear  wave  equation  to 

9  y 

a  =  (1  -  y^)"  (6) 


Second  Order  Analysis 


Using  the  first  order  solution  4>2^  to  form  the  source  terms  exciting 
the  second  order  potential  4>2  obviously  leads  to  products  of  Bessel  func¬ 
tions.  Such  terms  are  not  amenable  to  conventional  techniques  for  finding 
particular  solutions.  Consider  instead  the  region  off-axis  where  yR  is 
■sufficiently  large  to  replace  the  Bessel  function  by  its  asymptotic 
expansion.  Such  a  region  exists  if  ka  is  sufficiently  large.  In  cases 
where  ka  is  not  large,  the  analysis  may  be  conceptualized  as  temporarily 
removing  the  walls  of  the  duct,  thereby  converting  the  system  to  an 
infinite  half  plane  in  z  >  0.  Carrying  out  the  analysis  in  the  fictitious 
region  where  yR  is  large  nevertheless  leads  to  a  response  which  satisfies 
the  boundary  conditions  for  the  original  problem. 


When  yR  »  1,  eq,  (5)  may  be  represented  as 


-  —  ^  sin(t  -  az)  cos(yR  -  7)  +  0[(yR) 


-3/2 


1  a  '-yR'  ^  4' 

Substituting  this  expression  into  the  second  order  part  of  the  wave 


equation  yields 


V^4>2  - 


3^4>^ 


3t' 


~  [(Sq  ■  sin(2yR)]  sin(2t  -  2az) 


TO  yR 


+  0((yR)"^] 


(7) 


(8) 


Only  the  part  of  the  particular  solution  which  exhibits  growth  with 
Increasing  z  need  be  evaluated.  The  remainder  of  the  particular  solution, 
as  well  as  the  complementary  solution,  remain  very  small  in  the  entire 
domain.  The  result  is  that  in  the  region  where  yR  »  1,  for  all  z, 

a 

4>2  - 5 —  z  cos(2t  -  2az)  3in(2yR)  +  0[(yR)”^]  (9) 

ATO'^yR 

The  foregoing  is  recognizable  as  the  asymptotic  expansion  for  large  yR  of 

2  2 

<t>2  *  -  2  cos(2t  -  2az)  [J„(yR)^  -  J  '(yR)'^]  +  0(1) 


(10 


where  0(1)  represents  terms  that  are  bounded  for  all  z.  This  expression 
is  descriptive  of  the  most  significant  part  of  the  second  order  potential 
at  all  pR. 


Renormalization 


Expressions  for  the  particle  velocity  components  v^,  Vp  and  the 
pressure  p  may  be  obtained  by  taking  the  appropriate  derivatives  of  the 
potential  function  as  obtained  from  eqs.  (4),  (5),  and  (9).  The  general 
form  of  these  variables  is 

Uj^  =  Efj^(t  -  az,  yR)  +  E^z  f2(t  -  az,  yR)  +  O(e^)  (11) 

where  f^^  and  f^  denote  bounded  functions.  These  expressions  lose  validity 
when  z  =  Od/E).  where  the  second  order  term  is  no  longer  small.  This 
situation  is  corrected  by  replacing  the  independent  variables  by  a  set  of 
strained  coordinates.  Because  time  appears  only  in  combination  with  az,  no 
separate  transformation  of  t  is  required.  The  general  form  of  this 
transformation  is 


az  =  C  +  £zgj^(C,n),  yR  =  n  + 


(12) 


Equations  (12)  are  substituted  into  the  expressions  described  by  eqs.  (11). 
Then  the  functions  g.  and  §2  may  be  determined  by  requiring  that  the  0(e2) 
terms  in  a  Taylor  series  expansion  remain  bounded  for  all  z.  The  result  of 
this  procedure  is  that 


az  =  5  +  — 2  z  cos(t  -  C)JQ(nR) 
2a 

yR  =  ri - j  z  sin(t  -  C)dg'  (^) 

2a 


(13) 


The  particle  velocity  and  pressure  resulting  from  this  transformation  are 


P  = 


*^0^0 


p  c  ^ 

V  =  E  -2-^  cos(t  - 
z  y  0 


(14) 


=  -  ecg  sin(t  -  5)JQ'(n)  +  O(e^) 


Discussion  and  Example 


Equations  (13)  and  (14)  jointly  describe  the  signal.  For  the 
purpose  of  interpreting  the  results  it  is  Instructive  to  first  observe 
that  the  coordinate  transformation  may  be  rewritten  as 


az 


C  + 


^0  /z, 

—  2  (— ) 

2a  0 


(15) 


Lines  of  constant  C  and  n  represent  wavefronts  and  rays  of  constant  phase, 
respectively,  for  the  signal  described  by  eqs.  (14).  Thus  eq.  (15) 
describes  a  process  of  self-refraction,  in  which  rays  (constant  n)  are 
bent  by  the  transverse  velocity.  This  is  in  addition  to  the  amplitude 
dispersion  which  shifts  the  wavefronts  (constant  e)  in  one-dimensional  waves 
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The  waveforms  in  Figure  1  describe  a  signal  in  air  when  the  funda¬ 
mental  nonplanar  duct  mode  is  excited.  The  corresponding  sound  pressure 
level  at  the  source  is  156.5  dB  (re  20  yPa)  and  ka  =  9.5.  Each  waveform 
is  plotted  in  retarded  time,  such  that  the  linearized  signal  for  each  case 
would  appear  in  the  figure  as  one  period  of  a  sine  curve  having  essentially 
the  same  amplitude  as  that  for  the  nonlinear  signal.  The  signal  at  R  =  6.0 
shows  a  predominant  second  harmonic  because  the  position  is  near  a  nodal 
ray  for  the  linear  signal.  This  is  one  of  the  effects  of  self-refraction. 

Figures  2  and  3  describe  the  amplitudes  P  of  the  first  six  harmonics 
as  a  function  of  the  source  amplitude  Pq.  Figure  2  (on-axis)  is  reminiscent 
of  the  result  for  planar  waves.  Figure  3  is  for  an  off-axis  location  near 
a  nodal  line  for  the  linear  theory.  The  higher  harmonics  grow  much  more 
rapidly  than  the  fundamental  here,  as  a  result  of  self-refraction  of  the 
nodal  ray.  The  gap  in  the  plot  of  P^  corresponds  to  a  thus  far  unexplained 
null  when  Pq  =  156  dB. 
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Fig  1.  Waveform  at  position  (z,r),  as  indicated 
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Nonlinear  King  integral  for  arbitrary  axisymmetric  sound  beams  at 
finite  amplitudes.  I.  Asymptotic  evaluation  of  the  velocity  potential 
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This  paper  initiates  the  derivation  of  a  general  analytical  model  for  nonlinear  effects  in  sound 
beams  driven  at  high  source  pressure  levels.  The  excitation  is  generated  by  a  planar  transducer 
that  is  in  harmonic  motion  in  an  arbitrary  axisymmetric  pattern.  The  analysis  develops  a 
perturbation  solution  of  a  nonlinear  equation  for  the  velocity  potential.  The  first-order  term, 
which  is  derived  with  the  aid  of  a  Hankel  transform  to  represent  the  transverse  dependence,  is  the 
King  integral  for  a  linear  sound  beam.  Using  this  integral  to  form  the  source  terms  exciting  the 
second-order  potential  leads  to  a  dual  Hankel  transform.  Reduction  to  a  single  integral  is  achieved 
with  the  aid  of  an  asymptotic  integration  following  Laplace’s  method.  The  second-order  term  that 
is  derived  in  this  manner  describes  the  tendency  for  the  second  harmonic  to  grow  with  increasing 
distance  from  the  source.  This  result  is  an  intermediate  step  in  the  overall  development,  because 
the  integrand  loses  validity  in  the  spectrum  of  transverse  wavenumbers  near  the  transition 
between  evanescent  and  propagating  wavelets,  as  well  as  for  increasing  distance  from  the 
transducer. 

PACS  numbers;  43.25. Cb 
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INTRODUCTION 

Recent  surveys'  ‘  have  noted  that  the  “infinite  baffle” 
problem  has  been  described  by  several  alternative  formula¬ 
tions.  The  specific  concern  in  this  subject  is  the  signal  gener¬ 
ated  within  a  fluid  by  small  amplitude  oscillations  of  a  trans¬ 
ducer  which  is  contained  within  an  infinite  planar  boundary. 
Such  results  are  valid  for  very  weak  signals,  in  which  case 
material  and  convective  nonlinearities  are  negligible  effects. 
Two  basic  formulations  of  the  linear  problem  are  the  Ray¬ 
leigh  and  King  integrals. 

The  Rayleigh  integral^  treats  the  signal  as  a  superposi¬ 
tion  of  spherical  wavelets  which  are  generated  by  infinitesi¬ 
mal  sources  on  the  transducer  face.  In  contrast,  the  King 
integral'*  results  from  a  Hankel  (Fourier-Bessel)  integral 
transform  transverse  to  the  axis  of  symmetry.  The  acoustical 
medium  in  such  an  analysis  becomes  a  waveguide  of  infinite 
diameter.  The  transducer  then  seems  to  generate  a  spectrum 
of  guided  planar  mode  wavelets  whose  strength  varies  with 
the  transverse  position.  The  significant  aspect  of  both  qua¬ 
drature  solutions  is  that  they  provide  a  convenient  frame¬ 
work  for  quantitative  evaluations  of  the  signal  at  any  loca¬ 
tion.  They  also  lead  to  analytical  approximations  that  are 
valid  in  certain  ranges,  such  as  the  Fraunhofer  (farfield)  re¬ 
gion. 

The  same  is  not  true  for  treatments  of  nonlinear  effects 
which  arise  when  the  transducer  is  driven  to  large  ampli¬ 
tudes.  One  type  of  analysis  of  this  question  was  performed  by 
Lockwood  era/.’  They  considered  the  case  where  the  excita¬ 
tion  is  reasonably  small,  so  that  nonlinearity  is  not  signifi¬ 
cant  in  the  nearfield.  Such  a  restriction  leads  to  a  farfield 
description  based  on  Lockwood’s  analysis  of  nonuniform 
spherical  waves.**  Obviously,  such  an  analysis  provides  no 
information  regarding  nonlinear  effects  in  the  nearfield. 


The  nearfield  was  the  specific  concern  of  the  analysis 
performed  by  Ingenito  and  Williams.^  They  employed  a  per¬ 
turbation  series  for  the  potential  function,  in  which  the  lead¬ 
ing  term  was  described  by  the  Rayleigh  integral.  That  result 
was  then  used  to  evaluate  the  source  terms  exciting  the  sec¬ 
ond-order  potential.  Neglecting  backscattering  at  the  second 
order  and  introducing  some  additional  approximations  then 
led  to  a  description  of  second  harmonic  formation  that  had  a 
quadrature  form. 

One  limitation  of  that  analysis  is  that  it  is  valid  only  for 
very  high  frequencies:  ka  >  100  according  to  Ref.  7.  Another 
shortcoming  is  one  that  is  often  encountered  in  perturbation 
analyses.  Specifically,  if  a  dependent  variable  is  expanded  in 
a  perturbation  scries,  then  the  results  are  only  valid  when  the 
second-order  terms  are  very  small  compared  to  the  first  or¬ 
der  ones.  The  analysis  performed  by  Ingenito  and  Williams 
indicated  that  the  second  harmonic  grows  with  increasing 
axial  distance,  whereas  the  first  harmonic  (i.e.,  linear  result) 
shows  no  such  growth.  It  follows  that  these  results  shed  light 
on  how  harmonics  begin  to  form,  but  further  extrapolation 
might  lead  to  errors.  (This  seems  to  be  the  case  for  their  Fig. 
1) 

Another  formulation  of  finite  amplitude  sound  beams, 
which  has  been  prominent  in  the  Soviet  literature,  employs  a 
modified  Burgers’  equation.  The  basic  assumption  made  in 
the  derivation  of  this  equation'  *  is  that  there  are  three  spa¬ 
tial  scales  for  the  signal.  The  shortest  scale  is  the  axial  wave¬ 
length  and  the  longest  scale  describes  the  development  of 
nonlinear  effects  axially.  The  intermediate  scale  describes 
the  variation  transverse  to  the  axis  of  symmetry. 

These  approximations  seem  reasonable  for  the  high-fre¬ 
quency  limit.  Unfortunately,  solutions  for  monochromatic 
transducer  motion  have  only  been  obtained  by  finite-differ¬ 
ence  techniques.  Typical  of  such  investigations  are  Refs.  10- 
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1 2,  which  all  seem  to  have  employed  the  same  (incompletely 
described)  computer  code.  It  is  significant  that  this  group  of 
investigations  have  only  considered  situations  where  the 
boundary  excitation  is  a  prescribed  pressure.  Extending 
those  analyses  to  cases  where  the  particle  velocity  on  the 
boundary  is  known,  as  is  the  case  for  most  transducers,  re¬ 
quires  a  relation  between  pressure  and  velocity  that  is  not 
contained  in  the  basic  theory.  Also,  it  should  be  noted  that 
the  small  scale  of  some  diffraction  effects'^  introduces  some 
doubt  regarding  the  length  scales  assumed  to  derive  the 
modified  Burgers'  equation. 

The  present  investigation  is  a  perturbation  analysis,  as 
was  the  work  by  Ingenito  and  Williams,  although  the  King 
integral  is  used  here  to  describe  the  first-order  term.  The 
analysis  is  founded  on  the  recognition  that  nonlinearities 
arise  in  two  forms  when  the  signal  level  is  moderately  high. 
Some  nonlinear  effects  maintain  their  level  or  die  out  as  the 
signal  propagates  away  from  the  transducer.  Typical  of  such 
an  effect  is  the  fact  that  the  transducer  face  represents  a  mov¬ 
ing,  rather  than  a  fixed  boundary,  for  satisfying  continuity  of 
particle  velocity. 

Order  of  magnitude  considerations  indicate  that  such 
effects  are  too  small  to  describe  the  levels  of  higher  harmon¬ 
ics  that  have  been  measured.  As  is  the  case  for  planar  and 
other  one-dimensional  waves,'*  significant  distortion  phe¬ 
nomena  stem  from  cumulative  growth  of  higher  harmonics. 
Such  action  is  a  result  of  the  fact  that  the  acoustic  medium  is 
nondispersive.  so  that  higher  harmonics  propagate  with,  and 
resonantly  interact  with,  the  primary  signal. 

The  present  analysis  consistently  accounts  for  cumula¬ 
tive  growth  effects.  The  sole  assumption  introduced  in  the 
course  of  the  analysis  is  that  the  nonlinear  mechanism  caus¬ 
ing  harmonic  generation  has  the  same  behavior  at  all  loca¬ 
tions  in  the  acoustic  field.  The  first  part  of  the  investigation, 
described  in  this  paper,  obtains  an  expression  for  the  first 
two  orders  of  approximation  of  the  velocity  potential.  The 
second  part  of  the  investigation”  will  employ  coordinate 
straining  transformations  to  correct  irregularities  in  the  re¬ 
sponse  associated  with  the  derived  potential  function.  The 
acoustic  signal  will  be  described  in  a  quadrature  form  that  is 
reminiscent  of  the  linear  King  integral.  A  quantitative  exam¬ 
ple  will  compare  the  harmonic  content  of  the  waveform  to 
measurements  recently  reported  by  Gallego- Juarez  and 
Gaete-Garreton. 

The  overall  analytical  procedure  may  be  traced  back  to 
the  author’s  previous  investigation  of  two-dimensional  radi¬ 
ation  from  a  boundary.'*’”  However,  the  use  of  complex 
functional  forms  and  the  introduction  of  Hankel  transforms 
to  treat  the  axisymmetrical  geometry  require  substantial  al¬ 
terations  from  previous  work. 


I.  FUNDAMENTAL  EQUATIONS 

The  propagation  speed  of  infinitesimal  planar  waves  is 
denoted  as  Cg,  and  (i3„  —  1)  is  the  nonlinearity  parameter  in 
the  pressure-density  relation  at  fixed  entropy^ 


p=pA  + - 

l\  \  Po  ' 


(1) 


Thus/^o  =  (}'  -I-  l)/2  for  an  ideal  gas,  where  is  the  ratio  of 
specific  heats.  Let  \R,  z]  denote  nondimensional  cylindrical 
coordinates  (axisymmetric),  with  R  =  0  corresponding  to 
the  center  of  the  transducer  which  is  on  the  boundary  c  =  0. 
Also,  denote  the  nondimensional  time  variable  as  r.  The  cor¬ 
responding  dimensional  position  coordinates  are  (R  /k,  z/k  i 
and  the  dimensional  time  is  t  /o),  -where  a>  is  the  frequency  of 
the  monochromatic  excitation  and  k  =  ci/c„  is  the  wave- 
number  for  a  nominal  planar  wave. 

The  dimensionless  velocity  potential  <1>  is  related  to  the 
particle  velocity  components  by 

dd>  36  . 

dz  3R 

The  boundary  condition  corresponding  to  the  axisymmetric 
motion  of  an  arbitrary  transducer  may  then  be  written  in 
complex  form  as 


36 

3z 


=  0 


=  —f(R  )exp(;7  )  +  c.c., 
li 


(3) 


where  f(R  )  is  any  complex  function  whose  magnitude  ap¬ 
proaches  zero  with  increasing  R.  In  general,  c.c.  will  be  used 
to  denote  the  complex  conjugate  of  the  preceding  term. 
(Forming  products  of  complex  functions  necessitates  retain¬ 
ing  conjugate  parts,  rather  than  identifying  only  real  parts.  I 
For  weakly  nonlinear  waves,  the  acoustic  Mach  number  €  is 
a  finite  parameter  with  |e|  •<  1 .  The  nonlinear  wave  equation 
governing  6  is‘° 

=  2(^0  - 

ot  “  at 

+  4-IV6-'^6)  +  0{6^).  |4a| 

3t 

where 

3-6  ,  \  36  ,  3-6 

V  0  - - r  H - T - \ - ♦ 

3R-  R  3R  3z-  (45j 

In  addition  to  Eq  (3),  the  other  boundary  condition  for  6  is 
that  the  signal  should  be  either  an  outgoing  wave  or  an 
evanescent  wave  at  large  z,  and  that  it  show  suitable  decay 
with  increasing  R. 

The  velocity  potential  is  expanded  in  a  perturbation  se¬ 
ries 

6  =  €6\  -b  €-62  -b  (5) 

The  equations  governing  6\  and  62  are  found  by  collecting 
like  powers  of  e  in  Eqs.  (3)  and  (4a).  The  first-order  equations 
are 


(6) 

^  j  =  )exp(/> )  +  c.c. 

dz  2i 

Equations  (6)  are  the  statement  of  the  linearized  prob¬ 
lem.  The  nonlinear  effects  are  contained  in  62  succeed¬ 
ing  terms.  A  complete  solution  for  62  requires  satisfaction  of 
the  boundary  conditions,  which  involves  evaluation  of  the 
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complementary  solution,  as  well  as  of  the  particular  solution 
associated  with  the  source  terms  arising  from  However, 
the  complementary  solution  is  bounded  and  therefore  repre¬ 
sents  a  noncumulative  O  (e^)  contribution  to  the  signal  at  all 
locations.  As  noted  earlier  such  effects  are  usually  insignifi¬ 
cant  compared  to  the  observed  levels  of  nonlinear  distortion. 
Thus  it  is  only  necessary  to  find  a  particular  solution  of  the 
second-order  equation.  The  first  of  Eqs.  (6)  provides  a  simple 
identity  for  The  resulting  second-order  equation  aris¬ 
ing  from  Eq.  (4a)  is 


dt- 


4-  V(b*^d 


(7) 


II.  LINEARIZED  SOLUTION 


Two  approaches  that  have  been  employed  to  solve  the 
linearized  problem,  Eqs.  (6),  are  the  Rayleigh  integral  and 
the  King  integral.  The  latter,  which  is  essentially  an  inver¬ 
sion  of  a  Hankel  (Fourier-Bessel)  transform,  is  more  suitable 
for  the  task  of  evaluating  Hence  let 


d, 


■»  -x 

Ju 


n(p^\n.z,t  \J„\nR  ]dn  -i-  c.c. 


(8) 


Substituting  this  expression  into  Eqs.  (6)  and  using  the  fsct 
that  J„[nR  )  is  a  solution  of  Bessel’s  equation  leads  to 


dz' 


d-<P, 

3i- 


-  n-0,  =  0. 


d<p, : 

dz  ! 

where 


=  f,  expi/r  I, 


F,=~  RfRiJ„\nR\dR 

ll.n 


(9) 


(10) 


is  the  transform  of  the  icomple.xi  spatial  excitation  function 

f\R  I. 

The  solution  of  Eqs.  (9)  is  a  propagating  wave  when 
0  <  n  <  1 ,  or  an  evanescent  wave  when  n  >  1 .  This  solution  is 

<A,  =  -  (F,/^jexp(/r (11) 
where 

/z;  =  n-  -  1.  (12al 


Satisfying  the  radiation  condition  as  2— oo  leads  to  the  fol¬ 
lowing  choice  for  the  branch  cuts: 


/|1  —  n‘)''%  0<n  <  1, 

|/|--l)''^  n>l. 


(12b) 


The  result  of  substituting  Eq.  (11)  into  Eq.  (8)  is 


^-^p-^exp(/r )<fn  +  c.c.  (13) 

This  is  the  King  integral  representation  of  the  potential  func¬ 
tion  for  the  linearized  sound  beam. 


III.  SECOND-ORDER  POTENTIAL 

The  first-order  solution  d^  in  Eq.  (13)  is  used  to  form  the 
source  terms  driving  dz  in  Eq.  (7).  Forming  quadratic  pro¬ 
ducts  of  d>  requires  that  different  symbols  be  used  to  repre¬ 
sent  the  transverse  wavenumber  forming  each  term  in  the 
product.  Also,  care  must  be  taken  to  include  the  complex 
conjugate  parts  in  the  product.  A  quadratic  product  of  sinu¬ 
soidal  terms  generally  leads  to  a  term  having  a  nonzero  mean 
value,  but  the  time  derivative  appearing  on  the  right  side  of 
Eq.  (7)  removes  such  an  effect.  The  result  is 


XJoimR  )Jo{nR  )  —  mnJ,{mR  )J](nR  )] 

XeJip[2it  -  \ju„  + /n„\z]dm  dn  +  C.C.,  (14) 

where  the  symmetry  of  the  integrand  has  been  exploited  to 
reduce  the  integration  over  the  first  wavenumber  to  a  finite 
range.  (This  introduces  an  additional  factor  of  two.) 

It  is  consistent  with  the  form  of  Eq,  (14)  to  try  to  con¬ 
struct  the  particular  solution  for  dz  as  the  sum  of  two  dual 
Hankel  transforms.  The  kernel  of  one  transform  would  be 
mn  JolmR  ]Jo{nR  )  and  the  kernel  of  the  second  would  be  mn 
J^(mR  )Ji{nR  ).  The  following  equivalent  form,  which  utilizes 
linear  combinations  of  the  aforementioned  kernel  functions, 
leads  to  significant  analytical  simplifications. 


02  dn,  (15a| 

0 

where 


0z  =  n(  m02i(z,t.m.n\[Jo(mR  ]Jt)[nR  ) 
Jo 


—  J^[mR  )  J^[nR  )\dm  -f  n  m0zz(z,t,m,n] 
Jo 


X  [Jol'ri/? )  Jo[nR  )  -I-  Jx{mR  ]J\(nR  )]dm.  (15b) 


The  following  identities,  which  are  derived  from  the 
recursion  relations  for  Bessel  functions,*'  are  useful  for  eval¬ 
uating  the  transverse  derivatives  of  dz- 

(^41-)'"'"''' 

=  -  (m*  -I-  n']Jo{mR  )Jo{nR  )  +  2mnJ,{mR  ]Ji(nR  I, 


=  2mnJo(mR  )Jo(nR  )-{m^  +  n^  -  4/R  -)Jf{mR  ] 
XJ^inR  )  -  2{m/R  )J(,(mR  )Ji{nR  ) 

—  2{n/R  )y,(mR  )• 


(16) 


In  view  of  these  relations,  the  result  of  substituting  Eqs  1  ^ 
into  Eq.  (14)  and  matching  the  integrands  on  either  side 
the  equality  sign  is 
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^  -  (m  +  [JoimR  UoinR  I  -  )y,M  )1 

z-  df  i 

^  _  (m  -  «|'022]  [JoimR  UoinR  I  +  l-'iM  )) 

dz-  df 

-  [)4/A  'iJ>\mR  )J,inR  )  -  2{m/R  }Jo{mR  )J,{nR  )  -  2in/R  )J,{mR  )Jo{nR  )](0,,  -  -P,,) 

=  _  4i{F„FJfi^  MJ[1^o  f^MmR  UoinR  )  -  mnJ^{mR  )J,{nR  )]exp[2/f  -  {ti„+fi„]z]^  c.c. 


If  the  last  bracketed  term  on  the  left  side  of  this  equation 
was  not  present,  it  would  be  a  simple  matter  to  match  like 
functions  of  R  on  either  side.  In  order  to  address  this  matter, 
note  that  the  functions  0,,  and  0;;  are  independent  of  R. 
Conditions  governing  them  in  a  specific  region  of  R  should 
be  applicable  for  all  R.  This  is  significant  because  the  brack¬ 
eted  term  causing  difficulty  in  Eq.  (17)  decreases  at  a  rate  1/ 
R  faster  than  the  other  terms  as  R  increases.  The  foregoing 
argument  suggests  that  because  the  term  is  negligible  at  large 
R,  It  should  be  negligible  in  the  evaluation  of  02,  and  0,2 
any  R 

The  validity  of  this  hypothesis  might  be  questioned  for 
situations  where  n  and  m  are  small.  In  that  case,  the  asymp¬ 
totic  decay  of  Bessel  functions  having  arguments  mR  or  nR 
might  be  approached  at  unacceptably  large  values  o(R.  This 
question  may  be  examined  by  using  the  series  expansions  of 
the  Bessel  functions  for  small  arguments.  Specifically,  when 
mR  <  1  and  n/?  <  1 ,  it  may  be  shown  that 

-^J,{mR  )J,[nR  )  -  2^JoimR  )Ji[nR  )  -  2-^t{mR  ) 

R '  R  R 

xJ„{nR  )  =  —  mn{\  —  l(m^  -f  n')R  ^ 

-t-  i(5m=n=  -m*  -  n*]R  ''-•••].  (18) 

Thus  the  third  bracketed  term  in  Eq.  ( 1 7)  is  O  [mn]  when  mR 
and  nR  are  small,  whereas  other  terms  in  that  equation  are 
order  unity.  Thus  the  troublesome  bracketed  term  should 
have  negligible  influence  in  this  region  also. 

Neglecting  the  aforementioned  term  has  a  physical  jus¬ 
tification.  Recall  that  in  the  King  integral  formulation,  the 
acoustic  signal  is  viewed  as  a  spectrum  of  modes  in  an  infinite 
waveguide.  The  wavenumbers  m  and  n  are  merely  param¬ 
eters  characterizing  the  transverse  rate  of  variation  of  these 
modes.  It  is  reasonable  to  expect  that  the  nonlinear  mecha¬ 
nism  generating  the  second-order  contributions  to  these 
modes  are  described  by  the  same  differential  equations  at  all 
values  of  R,  and  for  all  values  of  m  and  n. 

When  the  third  bracketed  term  in  Eq.  (17)  is  ignored, 
matching  like  functions  of  R  in  that  equation  leads  to 

<9-0,,  5^0,, 

^  -{m  +  /i)*0:, 

dr  df 

=  -  -  I  -FnFm  -b  nm) 

xexp[2if  -  ->r n„)z]  +  C.C., 

d^<p22  ^<t>22 

=  -  -  1  -  nm) 

Xexp[2/r-(/i„  -t-/i„)r]  -t-c.c. 
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The  virtue  of  constructing  the  solution  for  6,  in  the  form  of 
Eqs.  (15)  is  now  evident;  any  other  form  would  not  have 
resulted  in  uncoupling  of  the  equations  for  the  transform 
variables  0,,  and  0;2- 

The  form  of  the  solution  of  Eqs.  (19)  is  suggested  by  the 
physical  implication  of  the  linear  King  integral,  which  con¬ 
sists  of  a  continuous  spectrum  of  modes  in  a  circular  wave¬ 
guide  of  infinite  extent.  The  axis  of  the  sound  beam  (i.e.,  the  z 
axis)  is  the  direction  in  which  these  modes  propagate  In 
general,  nonlinear  generation  of  harmonics  increases  with 
increasing  propagation  distance.  Hence  the  particular  solu¬ 
tions  may  be  written  as 

0,,  =  u(2,m,/i)exp[2if  -  iij.„  +  tj.„]z]  +  c.c., 

(20) 

0;2  =  h(z,m,rt)exp[2;t  -  (m,  +  c.c.. 

where  cumulative  growth  will  be  manifested  by  increasing 
values  of  the  amplitudes  a  and  b.  Substitution  of  Eqs.  (20) 
into  Eqs.  (19)  leads  to  a  set  of  uncoupled  ordinary  differential 
equations  for  these  amplitudes. 

■^-2(A4,  +  H„]‘ -  (m  ^  n]‘ +  *]a 

=  -  2t^'’^'"(^o-  1  +  nm),  (21a) 

-^-2(^4,  -f-/4„)^-(m-n)=-f  4]A 

=  -21^"^"' [00-  1  -  nm).  (21b) 

At  this  stage,  it  is  appropriate  to  recall  that  the  analysis 
of  (ft,  requires  evaluation  of  only  the  portion  that  exhibits 
growth.  If  the  values  of  m  and  n  are  such  that  the  coefficient 
of  u  or  h  in  Eqs.  (21)  does  not  vanish,  then  the  particular 
solution  is  independent  of  z.  In  contrast,  if  this  coefficient 
should  vanish,  then  the  corresponding  particular  solution 
for  a  or  h  is  proportional  to  z.  (Vanishing  of  the  coefficient  is 
equivalent  to  secularity  in  perturbation  analyses  of  nonlinear 
vibrations.)  It  is  found  with  the  aid  of  Eq.  (12b)  that  the 
condition  of  a  vanishing  coefficient  only  occurs  when  m  =  ri 
in  Eq.  (19a).  Therefore  only  the  contribution  of  a  needs  to  be 
evaluated.  The  magnitude  of  b  is  bounded  at  all  z,  which 
means  that  b  does  not  represent  a  cumulative  distortion  ef¬ 
fect. 

IV.  INTEGRATION  BY  LAPLACE’S  METHOD 

The  condition  where  the  solution  to  Eq.  (2 la)  grows 
with  increasing  z  has  been  shown  to  arise  as  m—n.  In  con¬ 
trast.  regions  far  from  m  =:  n  represent  contributions  that  do 
not  change  in  overall  magnitude  with  increasing  distance. 
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t,*  *  ^  aC.. 


The  contribution  of  the  region  around  m  =  n  may  be  deter¬ 
mined  by  following  Laplace’s  asymptotic  integration  meth¬ 
od"  based  on  an  expansion  using 

m  =  n  —  qA,  J<1,  9  =  0(ll,  (22a| 

where  J  and  q  are  positive  because  m  <  /t  for  the  integrals  in 
Eqs  il5i.  Note  that  J  is  a  fixed  parameter  indicating  the 
scale  of  the  difference  between  m  and  n.  The  Taylor  series  for 
the  coefficient  defined  in  Eq.  1 12b|  is  found  to  be 

-  9'^  +  •••)•  t22b) 

These  expressions  for  m  and  fi„  lead  to  the  following 
representation  of  Eq.  (21a)  in  the  region  where  m::zn: 


-  [^’d  V/i;  -I-  C»(d  -'ija,  =  -  2100^1/^1  ■  (22c) 

Now  observe  that  when  ^  =  0  (m  =  n),  the  particular  solu¬ 
tion  of  this  equation  is 

(23) 

In  contrast,  the  general  solution  for  a  when  q^O  has  the 
form 


Equation  (29)  gives  the  behavior  of  a  in  the  region  where 
m~n.  According  to  Laplace’s  integration  technique,  the  be¬ 
havior  at  large  z  may  be  found  by  using  that  relation  only.  In 
order  to  demonstrate  this  feature,  the  region  of  integration  is 
broken  into  two  intervals.  The  boundary  between  these  in¬ 
tervals  is  defined  to  be  m  =  n  —  <5,  where  <5  is  a  small  finite 
quantity.  Then 


J-n  -  A  ^1. 

d<P.  -(-  d02- 

D  Jn  -d 


(31) 


Because  the  first  domain  of  integration  does  not  contain  the 
secularity  condition,  the  oscillatory  nature  of  the  integrand 
results  in  boundedness  of  the  first  integral  for  all  z.'^  This 
term,  like  comparable  effects  that  occurred  earlier,  does  not 
contribute  to  the  cumulative  distortion  process,  so  it  is  dis¬ 
carded.  In  essence,  this  region  features  destructive  interfer¬ 
ence  between  the  m  and  n  waves. 

In  order  to  evaluate  the  second  integral,  Eqs.  (22)  and 
(29)  are  substituted  into  Eq.  (30).  The  difference  between  m 
and  n  is  less  than  <5  in  the  region  of  integration,  so  m  may  be 
replaced  by  n  in  the  Bessel  function.  Similarly,  may  be 
replaced  by  /i„  in  the  exponential  term  in  Eq.  (30).  The  inte¬ 
gral  arising  from  Eq.  (3 1 )  then  simplifies  to 


a  =  A,  exp(o,z)  +  exp(<7^)  +  2i0oP'l/q-A  (24) 

where  the  coefficients  cr,  and  <7,  are  the  roots  of  the  charac¬ 
teristic  equation  governing  the  complementary  solution: 


cr  -241  „[2  -  nq A  /fx\  -  q'A  -/2fi* )a  -  q'A  V/rJ  =  0. 

(25) 

Solving  this  quadratic  equation  yields 

=  -  q'A  -h  0{A  ^),  <72  =  4^,  +0\A  ). 

(26) 

As  ^—4),  Eq.  (24)  must  approach  Eq.  (23).  Because  <72  is 
Oil)  and  Eq.  (23)  has  no  term  that  varies  exponentially  in  z, 
set  .-1;  =0.  Also,  because  <7,  is  0[A  ’),  exp(cr,z)  may  be  re¬ 
placed  by  the  leading  terms  in  a  series.  The  '-ondition  that 
Eq.  (24)  approach  Eq.  (23)  then  leads  to 

lim[4,(l  +  <7,z|  +  2i0oF\/q-A =  [ipj^\/2iii]z.  (27) 

1  ■« 


The  first  of  Eqs.  (26)  shows  that  this  condition  is  satisfied  at 
all  values  of  z  by 

A,=  -2i0,^UrA-.  (28) 

Substituting  this  expression  into  Eq.  (24)  yields  the  general 
solution 

a  =  {2i0oFl/9'‘^  ’)[  J  -  eJ'P)  -  9^^  V4^’)].  (29) 

The  next  step  is  to  evaluate  the  total  contribution  to  02 
of  the  value  of  a  associated  with  ail  wavenumbers  m.  For  this 
determination  the  first  wavenumber  n  is  held  constant  at  an 
arbitrary  value.  The  combined  effect  is  defined  by  Eqs.  (IS) 
and  (20)  to  be 

02  =  n  r  w{a(z,m,n)exp[2/r  -  (/<„ -(-/u,„)z]-)- c.c.} 

Jo 

X  IJoimR  )J„{nR  )  -  )]dm  -f  0(1). 

(30) 


02  =  {2i0on'FlE (SAA//^i)exp{2it  -  2n„z]  +  c.c.} 


xU„(n/J  ?-J,{nR  )=], 


(32) 


where 


-  txp(  -  q^A'z/An),] 
q^A  ’ 


(33) 


Evaluation  of  the  function  E  introduces  a  square  root  of 
z/fi]„  for  which  it  is  important  to  account  for  the  fact  that^„ 
is  imaginary  when  n  <  1.  Specifically, 


(34) 


where  an  overbar  denotes  the  complex  conjugate  of  the 
marked  quantity.  Integration  by  parts  of  Eq.  (33)  then  gives 


\  J  \  24i„fi„  / 

-(1/<5)[1  -exp(  -5  V4/iJ)], 


(35) 


where  erf  denotes  the  error  function. 

One  noteworthy  feature  of  Eqs.  (32)  and  (35)  is  that  the 
only  remaining  parameter  associated  with  the  asymptotic 
integration  is  the  integration  limit  S,  which  is  finite  value. 
Consider  Eq.  (35)  as  z  increases  while  S  is  held  fixed.  An 
expansion  of  the  error  function  for  large  arguments  leads  to 


S(irfz„z)''^ 


(36) 


E{S4iA/nl)~[{irfi„z)''^24i,fi„]  -  l/S. 


The  growth  effect  comes  from  the  first  term  in  E  above. 
In  general,  the  behavior  as  z— »oo  is  said  to  be  the  “dominant 
term.”^^  The  dominant  term  in  02  originates  from  the  por¬ 
tion  of  the  particular  solution  associated  with  the  region  of 
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seculanty.  mzzn.  Subdominant  terms,  such  as  the  particular 
solution  associated  with  m^n,  have  already  been  neglected 
because  they  do  not  represent  a  growth  effect.  Thus  the  func¬ 
tion  E  in  Eq.  i32i  may  be  replaced  by  its  dominant  term,  as 
given  by  Eq.  i36|  when  \/8  is  neglected.  When  the  resulting 
e.xpression  for  0,  is  used  to  form  the  second-order  potential 
according  to  Eq.  |15a|,  and  then  combined  with  the  first- 
order  potential  in  Eq.  1 13),  the  result  is 


6  = 


r 

Jo 


0  dn. 


where 


(37a) 


0  =  _  - ^l„z'Jo[nR  )  -I-  ripQ—^-" 

Xe.\p|2ir  -  lfi^z){Jn\nR  1^  -  J^nR  )’] 

-I-  c.c.  -I-  subdominant  terms.  (37b) 


V,  CLOSURE 

.\n  expression  for  the  pressure  may  be  obtained  by  dif¬ 
ferentiating  Eqs.  (371  with  respect  to  time,  but  the  result  has 
some  problematical  aspects.  First,  it  is  clear  that  the  0{e^) 
term  grows  as  c'  ^  while  the  0[e]  term  remains  bounded. 
Thus  the  second-order  term  satisfies  the  smallness  assump¬ 
tions  inherent  to  a  perturbation  series  only  when  z  is  small 
compared  to  l/f’.-^  From  a  practical  viewpoint  this  limits 
the  validity  of  the  result  to  distances  that  are  very  small  com¬ 
pared  to  the  shock  formation  distance. 

Another  aspect  relating  to  the  validity  of  the  result  is 
less  obvious  Consider  the  situation  where  n— *1,  in  which 
case/i„— 0.  The  expression  for  pressure  derived  from  Eq. 
|37)has^„  in  the  denominator,  so  it  is  singular  at  n  =  l.The 
key  aspect  of  the  singularity  is  that  the  O  (e)  term  will  contain 
a  factor  l//r,,  while  the  O(e^)  term  will  contain  a  factor 
l//i,/i!,  r  Thus,  theO(f^)  term  grows  more  rapidly  than  the 
O  (el  term  as  n—  1 .  This  is  another  instance  where  the  magni¬ 
tude  of  the  second-order  term  grows  relative  to  the  first- 
order  term.  As  was  true  for  thez'^*  dependence,  nonuniform 
accuracy  limits  the  usefulness  of  the  pressure  derived  from 
Eqs.  137). 

The  lack  of  uniform  accuracy  in  z  is  not  surprising,  be¬ 
cause  It  IS  the  equivalent  of  secular  terms  in  nonlinear  vibra¬ 
tion  analyses.  The  nonuniform  accuracy  in  the  wavenumber 
n  IS  a  result  of  the  analytical  procedure  that  was  followed. 
The  integration  by  Laplace's  method  assumed  that  \n„  \  is 
not  very  small.  This  is  most  clearly  indicated  in  Eq.  (22b|, 
where  the  truncation  of  the  series  expansion  is  appropriate 
only  if  qA  Only  very  small  values  of  q  satisfy  this 

criterion  when  n— »1  (/i,— »0).  Therefore,  the  contribution  to 
the  second-order  potential  from  the  region  around  n  =  1  is 
not  well  described  asymptotically. 

There  are  other  shortcomings  in  the  form  of  Eq.  (37b). 
First,  the  (7  (f  I  term  is  the  same  as  that  obtained  from  linear 
theory,  i.e.,  it  is  the  conventional  King  integral.  Thus  the 
relation  does  not  indicate  that  there  is  depletion  of  the  funda¬ 
mental  harmonic  as  energy  is  transferred  to  higher  harmon¬ 
ics.’’  Another  important  limitation  is  that  Eq.  (37b)  de¬ 
scribes  only  the  second  harmonic,  but  higher  harmonics  are 
known  to  be  significant  to  the  distortion  process. 


The  aforementioned  items  lead  to  concern  regarding 
the  validity  of  any  prediction  of  pressure.  This  would  cer¬ 
tainly  be  the  case  if  Eqs.  (37)  were  to  be  used  directly.  The 
analysis  in  the  next  part  of  this  investigation’’  overcomes 
these  difficulties.  It  treats  the  response  obtained  from  Eqs. 
(37)  as  the  asymptotic  representation  for  small  ez'  '  and 
n^  \  o(  functional  forms  that  are  uniformly  accurate.  Such 
an  analysis  is  not  applied  directly  to  the  velocity  potential 
because  there  are  situations  where  a  portion  of  the  potential 
may  exhibit  nonuniform  growth  while  the  pressure  and  oth¬ 
er  state  variable  do  not.'*’ 
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Nonlinear  King  integral  for  arbitrary  axisymmetric  sound  beams  at 
finite  amplitude.  II.  Derivation  of  uniformly  accurate  expressions 

Jerry  H.  Ginsberg 

School  o/ Mechanical  Engineering,  Georgia  Institute  of  Technology,  Atlanta.  Georgia  30332 

(Received  IS  August  1983;  accepted  for  publication  23  May  1984) 

The  first  part  of  this  investigation  [J.  Acoust.  Soc.  Am.  76, 1201-1207  (1984)]  derived  a 
perturbation  representation  of  the  velocity  potential  for  an  axisymmetric,  monochromatic  sound 
beam  in  the  form  of  a  complex  function.  That  result,  which  used  Hankel  transforms  to  describe 
the  dependence  on  transverse  position,  lacked  uniform  validity  in  the  corresponding 
wavenumber,  as  well  as  the  distance  from  the  transducer.  Expressions  for  the  state  variables 
derived  from  the  potential  have  the  same  behavior.  The  example  of  a  planar  wave  is  used  to  adapt 
the  singular  penurbation  method  of  renormalization  to  the  case  where  the  potential  is  in  complex 
form.  The  resulting  technique  is  used  to  obtain  a  coordinate  straining  transformation  that  makes 
the  state  variables  uniformly  accurate.  The  expression  for  the  pressure  is  similar  to  the  King 
integral  in  linear  theory,  except  that  the  integrand  is  a  function  of  the  strained  coordinates. 

Comparison  of  the  predicted  waveform  properties  with  experimental  data  [Gallego-Juarez  and 
Gaete-Garreton,  J.  Acoust.  Soc.  Am.  73,  761-765  (1983)]  shows  good  correlation.  Further 
evaluations  disclose  some  new  features  of  the  distortion  phenomena  in  both  the  time  and 
frequency  domains. 

PACS  numbers:  43.25. Cb 


INTRODUCTION 

general  analysis  of  axisymmetric  sound  beams  was 
initiated  in  Part  I  of  this  study.  ‘  The  overall  goal  of  the  inves¬ 
tigation  IS  to  develop  a  theory  that  can  be  used  to  evaluate 
how  vanous  features  of  harmonic  transducer  motion  are 
manifested  within  the  entire  acoustic  field.  This  intent  obvi¬ 
ates  the  use  of  theories  founded  on  effects  that  arise  only  in 
the  near-  or  farfield. 

It  was  found  in  Part  I  that  combining  a  consistent  non¬ 
linear  wave  equation  with  a  Hankel  integral  transform  and 
an  asymptotic  integration  technique  led  to  a  representation 
of  the  velocity  potential  in  the  form  of  a  perturbation  series. 
The  second  part  of  the  analysis,  presented  here,  will  utilize 
that  potential  function  to  derive  an  integral  expression  for 
the  acoustic  pressure.  This  integral  reduces  to  the  King  inte¬ 
gral  of  linear  theory-  in  the  limit  of  an  infinitesimal  source 
pressure  level. 

Numerical  evaluation  of  this  integral  permits  compari¬ 
son  of  the  theoretical  predictions  to  those  in  one  series  of 
experiments.  In  addition  to  providing  validation,  the  exam¬ 
ple  will  yield  some  insight  into  the  unusual  distortion  phe¬ 
nomena  that  have  been  observed  in  sound  beams.^ 

I.  NONUNIFORMLY  ACCURATE  EXPRESSIONS 

In  the  first  part  of  this  investigation  z  and  R  were  nondi- 
mensional  cylindrical  coordinates  for  a  transducer  centered 
at  r  =  R  =  0,  and  t  was  nondimensional  time.  The  corre¬ 
sponding  dimensional  quantities  are  obtained  by  using  I /to 
to  define  the  time  scale  and  l/k  =  C(/a)  to  define  the  length 
scale,  where  co  is  the  frequency  of  the  monochromatic  mo¬ 
tion  of  the  transducer  and  c„  is  the  linearized  speed  of  sound. 
The  transducer  motion  was  written  in  complex  form  as 

f,  =  (f/2/)c„/(R  )  exp(;/ )  +  c.c.,  ( 1) 


where /(R )  is  an  arbitrary  complex  function  whose  magni¬ 
tude  goes  to  zero  as  The  acoustic  Mach  number  e  is 

small  compared  to  unity,  so  it  is  a  convenient  parameter  for  a 
perturbation  series.  In  general,  c.c.  will  denote  the  complex 
conjugate  of  the  preceding  terms. 

A  Hankel  transform  was  used  to  describe  the  first-order 
effects  transverse  to  the  axis  of  the  sound  beam.  For  a  speci¬ 
fied  transverse  wavenumber  n,  which  is  the  parameter  for 
the  transform,  there  is  a  corresponding  Hankel  transform  F„ 
of  the  shape  function /(R  ), 

F,  rR/(R)yo(nR)dR  (2) 

2/  Jo 

as  well  as  a  (nondimensional)  axial  wavenumber  ,  where 
the  radiation  condition  is  satisfied  by 

/i,  =/(l -n^)''^  0<n<l. 

(3) 

/I,  =(»^- 1)*'^  tt>l. 

As  a  result  of  employing  a  Hankel  transform,  the  first- 
order  velocity  potential  was  expressed  in  terms  of  the  King 
integral.  This  first-order  result  led  to  an  inhomogeneous 
wave  equation  for  the  second-order  potential.  An  asympto¬ 
tic  integration  by  Laplace’s  method  was  crucial  to  the  eva¬ 
luation  of  second-order  effects.  The  expression  for  the  veloc¬ 
ity  potential  that  was  derived  by  this  method  retained  only 
the  part  of  the  second-order  terms  that  displayed  cumulative 
growth  of  higher  harmonics.  Thus,  it  did  not  evaluate  any 
term  whose  magnitude  is  0  (e^)  at  all  values  ofz.  The  derived 
velocity  potential  was  written  as 

(fi  =  j  0dn,  (4) 

where,  for  l/r,  |’'SO(e), 
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nF,  _  ,,, 

0  =  -  € - explir  -fi^z\Jr,(nR  |  -i-  i€-0q — - 

X  exp(2/r  -  ln„z\\Ja[nR  )'  -  J^(nR  )-]  +  c.c.  +  SDT, 

(5) 

where  SDT  represents  subdominant  terms  that  do  not  grow 
with  increasing  z.  and  an  overbar  denotes  the  complex  conju¬ 
gate.  Note  that/?,,  is  the  coefficient  of  nonlinearity,  equal  to 
My  -I-  1 1  for  an  ideal  gas,  where  y  is  the  ratio  of  specific  heats. 

The  limitation  on  the  value  of fi„  for  Eq.  (5|  means  that 
the  expression  is  not  applicable  in  the  vicinity  of  n  =  1 .  This 
is  a  consequence  of  the  asymptotic  expansion  that  led  to  Eq. 
|5|,  for  which  it  was  assumed  that/x,  =  (?(1|.  Equation  (5)  is 
said  to  lack  uniform  accuracy  in  both  n  and  z,  because  the 
magnitude  of  the  second-order  term  grows  relative  to  the 
first-order  term  as  n—*  1  and  as  z— ►  ao .  The  primary  task  here 
is  to  obtain  from  Eq.  |5)  expressions  that  are  uniformly  accu¬ 
rate. 

Differentiating  Eqs.  (4)  and  (51  yields  expressions  for  the 
particle  velocity  components  and  pressure.  As  suggested  by 
the  form  of  Eq.  |4),  let  K. ,  and  P  denote  the  contribution 
of  a  specific  wavenumber  n  to  the  axial  velocity  v,,  trans¬ 
verse  velocity  Vn ,  and  pressure  p,  according  to 

V.  =  Co  r  y,  dn,  Ur  =  Co  I  Fr  dn,  (6) 

Jo  Jo 

p  =  P  dn.  (7) 

In  the  evaluation  of  the  expression  for  V. ,  the  factor  z''''  in 
the  O(e^)  term  of  0  may  be  considered  to  be  constant,  be¬ 
cause  its  derivative  leads  to  an  term,  which  de¬ 

creases  in  importance  as  z  increases.  The  physical  variables 
associated  with  Eqs.  (5)  and  (6)  are  then  found  to  be,  for  n  not 
close  to  1 , 

y 

’  dz 

=  enF^  e\^it  - p„z]\Jo[nR  1  -  2/e/7„(/jF,//2J 
X(;r/x„z|'  '-exp(/f-^,z)[yo(n/{  )•  -J^nR  t 


Fr  =  —  =  -£\n-F„/p„)t%p(it -p„z]\J'a(nR] 
oR 

-  li€po[nF„/Ji„][iTp„z)''-  exp(it-p„z} 

X  ; l»^ )]  I 

H-  c.c.  -I-  SDT,  (9) 

P=(i/MJK.  (10) 

It  is  convenient  for  later  operations  to  denote  the  derivatives 
of  Bessel  functions  by  primes,  rather  than  making  use  of 
identities  for  derivatives. 


II.  CORRECTION  OF  THE  DEPENDENCE  ON  THE 
TRANSVERSE  WAVENUMBER 

Equations  |8)-(  10)  suffer  from  the  same  lack  of  uniform 
accuracy  as  that  associated  with  Eq.  |S).  Such  behavior  is  not 
acceptable  for  state  variables  such  as  particle  velocity  and 


pressure.'*  The  recognition  that  the  integration  regarded 
'ji„  1  to  be  substantially  larger  than  zero  is  crucial  to  correct¬ 
ing  the  dependence  on  n.  It  leads  to  the  conclusion  that  Eqs. 
(8)-(  10)  are  the  asymptotic  representation  for  nonsmall  p  „  of 
alternative  functional  forms  that  behave  properly  as  — y). 
Note  that  each  0(rj  term  contains  a  factor  /xy  '  ’  higher 
than  the  corresponding  O  (e )  term.  It  follows  that  the  alterna¬ 
tive  forms  must  feature  a  function  whose  expansion  for  large 
p„z  is  proportional  to p~ 

A  variety  of  functions,  such  as  Bessel  functions,  are  pos¬ 
sible  candidates  in  this  regard.  However,  recall  that  p„  is 
either  real  or  imaginary,  depending  on  the  value  of  n.  Most 
functions  whose  asymptotic  behavior  is  appropriate  for  real 
p„  introduce  new  singularities  for  imaginary  p„,  or  vice 
versa.  The  only  function  that  was  found  to  be  acceptable  for 
all  values  of  n  was  the  complementary  error  function.  Spe¬ 
cifically,  it  is  known’  that  for  large /x,z, 

erfc[(^„z)''-] 

=  (ir/x„z)“''-exp(  -Mn4[l  ^0{\/p„z\].  Ill) 

Using  Eq.  (1 1)  to  recast  the  O(e^)  terms  in  Eqs.  (8)-l  10) 
leads  to 

V,  -  £nF„  exp(/r -p„z\\Jo{nR  )  -  liirePa\p„/p„) 

XnF„ztdc[\j:i„z)''-]  expl/rjl/olrt/f )'  -Jt[nR  )’]  | 

-I- c.c. -f  SDT,  (12) 

Fr  =  -  e  ^—2.  exp(/r  - p„z]\J o(nR)-  ) 

Pn 

X  nF,  z  erfc  [  (/X  „  z) '  ^  ]  exp(// )  [  yo(nF  )y  i  («/? ) 

-  y,(/i/J)y :(;«/?)]  I -h  c.c. -H  SDT,  (13) 

F=(///xJF,.  (14) 

Note  that  the  coefficient  ;x,//x„  is  merely  -(-  1  or  —  1 
depending  on  whether /x„  is  real  (n>  1 )  or  imaginary  (n<  1). 
Thus  Eqs.  ( 1 2)-(  14)  have  the  same  degree  of  singularity  in  ,ix , 
for  the  first-  and  second-order  terms:  they  are  descriptive  of 
the  response  for  all  n. 

III.  DERIVATION  OF  RENORMALIZATION  USING 
COMPLEX  VARIABLES 

Although  Eqs.  (12)-(14)  are  valid  for  all  n,  they  still  are 
not  uniformly  accurate  for  all  z.  One  method  for  correcting 
this  situation  is  the  renormalization  version  of  the  method  of 
strained  coordinates.'’  The  genera)  basis  for  this  method  is 
the  argument  that  nonuniform  behavior  relative  to  one  of  the 
independent  variables  (space  or  time)  is  the  result  of  improp¬ 
er  truncation  of  a  Taylor  series  expansion  in  powers  of  e.  A 
much  simpler  example  is  to  expand  sin(z  -i-  ez)  in  powers  of  e, 
and  to  truncate  such  an  expansion.  The  original  function  is 
periodic,  but  the  truncated  representation  shows  cumulative 
growth  in  ez. 

As  a  consequence  of  the  foregoing  argument,  it  may  be 
anticipated  that  there  is  a  transformation  of  the  space-time 
variables  for  which  the  response  does  not  exhibit  cumulative 
growth.  The  difference  between  the  magnitude  of  the  phys¬ 
ical  variables  and  the  transformed  ones  will  grow  cumula¬ 
tively  in  the  same  manner  as  the  O  (e^)  terms  in  Eqs.  1 1 2)-|  14). 
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The  transformed  variables  therefore  represent  a  straining  of 
the  physical  space-time  grid. 

All  prior  applications  of  the  renormalization  technique 
seem  to  have  employed  real  functional  forms.  In  order  to 
adapt  the  method  to  the  complex  functions  in  Eqs.  (12H14), 
it  is  useful  to  consider  the  analogous  steps  for  a  finite  ampli¬ 
tude  planar  wave.  Consider  such  a  wave  propagating  in  thez 
direction,  due  to  a  harmonic  particle  velocity  imposed  at 
z  =  0.  The  potential  function  in  this  case  is  governed  by 


dr 


d'6  _  y  I 

—  -  -1/^0  - 


0{<t\ 


—  =  (^— )  exp(/r  I  +  C.C., 

dz  .>0  \2iJ 

where  'le  scales  for  nondimensionalizing  position  and  time 
are  the  same  as  those  for  the  sound  beam. 

Carrying  out  a  perturbation  analysis  based  on 

d  =  fd,  +  rdi  +  •••  (16) 

yields  a  homogeneous  wave  equation  for  which  is  solved 
for  a  wave  propagating  in  the  positive  z  direction.  This  solu¬ 
tion  IS  used  to  form  the  source  terms  in  an  inhomogeneous 
wave  equation  for  dy  The  particular  solution  of  the  latter  is 
readily  obtained  by  elementary  methods,  with  the  result  that 

6  =  If  exp[/(f  -  z|]  -I-  exp[2i(/  -  z)] 

c.c  +  O  (e^),  (17) 

where  the  O  le^)  terms  not  appearing  explicitly  are  bounded 
for  all  z. 

The  expressions  for  particle  velocity  and  pressure  cor¬ 
responding  to  Eq.  ( 1 7)  are 

V,  =  Co^  =  exp[i(f  -  z)l 
dz  2i 

X  1 1  +  le^„r  exp[/(r  -  z)]  |  +  c.c.  +  0(€^).  (18) 


P=  -P<y^o—=PiPoP,- 
dt 


These  are  not  uniformly  accurate  because  the  O(e^)  term 
grows  with  increasing  z  relative  to  the  O  (f )  term.  It  is  postu¬ 
lated  that  the  dependence  on  the  strained  coordinate  varia¬ 
bles  IS  uniformly  accurate.  Because  t  only  occurs  in  conjunc¬ 
tion  with  z,  it  is  adequate  to  strain  only  the  space  variable. 
The  strained  coordinate  f  approaches  z  when  e— so  the 
transformation  is  considered  to  have  the  form 


z  =  C  +  e[5(i,V)+5(^.r)I  +....  (19) 

The  task  now  is  to  evaluate  the  function  S.  The  first  new 
feature  of  the  complex  variable  formulation  of  the  planar 
wave  is  that  S  is  considered  to  be  complex,  in  agreement  with 
the  form  ofEq.  |18).  It  is  necessary  to  introduce  the  complex 
conjugate  of  5  in  Eq.  (19)  because  the  transformation  must  be 
real. 

Presumably,  expressing  u,  and  p  in  terms  of  ^  rather 
than  z  leads  to  uniform  accuracy.  This  should  be  the  result  of 
substituting  Eq.  (19)  into  Eq.  (18).  Because  the  earlier  analy¬ 
sis  considered  only  the  first  two  orders  of  e,  enforcing  unifor¬ 
mity  will  involve  considering  only  those  powers  in  a  series 
expansion.  First  note  that  for  e<  1, 


exp(  —  iz)  =  exp[  -  /(f  +  €S  +  eS )] 

=  exp(  -  i; )[  1  -  If (5  -I-  S )]  -I-  O  (f^j  (20) 

Substitution  of  Eq.  (19)  into  Eq.  (18)  in  conjunction  with  Eq. 
(20)  leads  to  the  following  e  and  f^  terms: 

p/pA  =  Vz/Co 

=  ( 1  /2/>  exp[i(i  -  b)]  1 1  -  i€{S  S ) 

-f  exp(/(r  -  f )]  1  +  c.c.  -i-  O  {r],  (21) 

where  O  (e^)  means  terms  having  that  order  of  magnitude  at 
all  z. 

Equation  (21)  will  be  nonuniform  in  ^  unless  the  f ^  term 
is  canceled.  This  is  the  criterion  for  the  straining  function  S. 
Thus,  set 

S  =  l\/2i)/3aScxp[i{t-n]-  (22) 

The  terms  remaining  in  Eq.  (21)  are 

-^  =  ;^exp[/(/-f  )](1  -(e5)-hc.c.  -i-O(r) 

P<fZ  2/ 

=-:^exp(i(r-^)]  H- -t- c.c.  +0(f^).  (23) 

2;  4/ 

Note  that  the  f^  term  appearing  above  is  imaginary.  Ac¬ 
counting  for  the  complex  conjugate  of  each  term  yields 

-^  =  ^  =  esin(r-^-)  +  0(r)  (24) 

PePo  ‘•0 

and  the  transformation  obtained  by  substituting  S  and  5  into 
Eq.  (19)  is 

z  =  f  +  e/9afsin(f-f).  (25) 

In  order  to  demonstrate  the  correctness  of  the  wave 
described  by  Eqs.  (24)  and  (25),  Eq.  (24)  is  used  to  eliminate 
the  sine  term  in  Eq.  (25); 

z  =  ${l+0o>’z/Co}.  (26) 

Solving  this  relation  for  ^  and  substituting  the  result  back 
into  Eq.  (24)  then  yields 

-^  =  ii  =  esinfi - £ - Y  (27) 

P(PI  Co  \  l^+^oCz/CojJ 

Except  for  the  fact  that  z  and  t  are  nondimensional  here,  Eq. 
(27)  is  identical  to  Eamshaw's  closed  form  solution  for  finite 
amplitude  planar  waves’  in  the  case  of  harmonic  excitation 
at  a  boundary.  The  same  result  using  real  functions  was  ob¬ 
tained  previously.* 

The  perturbation  analysis  of  the  velocity  potential  pro¬ 
vided  an  indication  of  how  higher  harmonics  tend  to  be 
formed.  Although  the  expressions  that  resulted  were  not 
uniformly  accurate,  the  information  the  analysis  provided 
was  sufficient  to  permit  the  coordinate  straining  procedure 
to  identify  the  more  general  signal.  The  exactness  of  the  re¬ 
sult  in  this  case  was  fortuitous.  In  general,  all  that  can  be 
expected  from  the  procedure  is  that  the  error  in  the  uniform¬ 
ly  accurate  prediction  will  be  no  larger  than  C2(e^)  at  any 
location. 
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IV.  APPLICATION  OF  RENORMALIZATION  TO  THE 
SOUND  BEAM 

Only  a  few  modifications  are  required  to  apply  the 
method  in  the  preceding  section  to  Eqs.  (12H14).  As  was 
true  for  the  planar  wave,  t  only  occurs  in  combination  with  z, 
so  no  transformation  of  the  time  scale  shall  be  introduced. 
However,  R  is  an  independent  coordinate,  so  a  straining 
transformation  for  that  spatial  variable  is  also  needed.  Thus 
the  coordinate  straining  lo  be  tned  is 

2  =  sn  ^  +  C.C.], 

(28) 

R=a.  -  c.c.]. 

Note  that  a  subscript  n  has  been  associated  with  the  strained 
coordinates  in  order  to  emphasize  that  the  strained  coordi¬ 
nate  describes  only  one  transverse  wavenumber  in  a  contin¬ 
uous  spiectrum. 

In  order  to  focus  on  the  conceptual  aspects  of  renormal¬ 
ization,  Eqs.  1 12|  and  ( 13|  are  written  in  the  standard  form, 

=  €fi\z.R.t.n\  +  e^zg,\z,R,t,n\  +  c.c.  +  0{r\; 
y=1.2.  (29) 

Substituting  Eqs.  |281  into  Eq.  (29)  and  retaining  the  e  and  e' 
terms  in  a  Taylor  series  expansion  leads  to 

V  oC, 

_  ^ 

-t- -r5;,)-- — 
da, 

+  c.c. -I- 0(r);  j=  1,2. 

(30) 

The  foregoing  is  analogous  to  Eq.  (21)  for  the  case  of  planar 
waves.  The  straining  functions  S.  and  must  annihilate 
the  nonuniform  g,  term  in  Eq.  (30).  Hence,  set 

~/ 1;,  -a,  .t.n I  ^  (i:,  ,a,  .t.n ) 

aa„ 

-  -i„g,\i„,a„t,n\\  j=\,l.  (31) 

It  is  necessary  to  find  functions  Sj  and  that  satisfy 
Eq.  |31I  for  both  values  of  j.  This  is  achieved  by  using  trial 
forms  that  are  suggested  by  comparing  the  functions  g,  to 
the  denvatives  of / .  The  actual  transformation  obtained  by 
satisfying  Eq.  (31)  and  then  forming  Eqs.  (28)  is 

;  =  C,  -  j  )  exp((/ ) 

xerfc{tu,i,‘J'  -]  +  c.c.\Jg(na,], 

(32) 

^  +  2tre^i*„  j  i(F„fi„/^„]np(it ) 

xerfcfl/z,^,)''-]  -)- c.c.iy,(naj. 

In  general,  the  terms  remaining  in  Eq.  (30)  when  Eq. 
|31I  holds  are 

y,  =  r^I,-^;^(f,,a„,r,n) 

+  +  c.c.  -I-  £7(f^).  (33) 


This  expression  is  comparable  to  Eq.  (23 1  for  a  planar  wave. 
The  term  in  that  case  was  imaginary,  so  that  no  such  terms 
appeared  explicitly  when  the  real  form  of  the  solution  was 
written.  The  expressions  for  the  sound  beam  are  more  com¬ 
plicated,  primarily  due  to  the  presence  of  Bessel  functions 
describing  the  transverse  variation  of  the  sound  beam.  Here, 
renormalization  removes  any  O  (e^)  terms  in  that  depend 
explicitly  on  time. 

In  other  words,  the  coordinate  straining  transformation 
derived  from  Eq.  (31)  is  based  on  matching  the  tendency  for 
generation  of  higher  harmonics.  Such  a  transformation  in¬ 
troduces  a  mean  value  over  one  period.  The  terms  appear¬ 
ing  in  Eq.  (33)  cancel  that  mean  value  at  all  locations. 

All  of  the  state  variables  may  be  obtained  from  Eq.  1 33 1. 
The  quantity  of  primary  interest  is  the  pressure,  for  which 
the  proportionality  in  Eq.  (14)  is  used  to  determine  P.  The 
contribution  of  all  transverse  wavenumbers  n  is  obtained  by 
integrating  according  to  Eq.  (6),  with  the  result  that 

/>=PoCofJ^  ^^^jexp(  -;u„^-jjexp(;r  )y„(naj 

-  2i7reP^/iF, erfc [ (^ , c  J ' ' - ] 

X[-^o(naJ*  -  \dn  +  c.c.  (34) 


Equations  (32)  and  (34)  jointly  describe  the  pressure  at 
specified  values  of  z,  R,  and  t.  Because  of  the  complicated 
nature  of  these  relations,  quantitative  evaluations  require 
numerical  algorithms  for  solving  the  coordinate  transforma¬ 
tion,  and  for  integrating  over  all  values  of  n.  When  n  <  1,  in 
which  case  is  imaginary,  evaluating  the  pressure  accord¬ 
ing  to  Eqs.  (32)  and  (34)  requires  computing  the  complement 
of  the  error  function  for  complex  arguments.  Useful  identi¬ 
ties’  for  this  task  are 
erfc[(/<r)''’] 

=2’''  exp(  -  iV/4){  [J  -  5,(cr)]  -h  /[J  -  C-lcr)]  | , 

(35) 

erfc[  — 

=2''-  exp(  -  /7r/4)(  [J  -  C,(<7(I  -|-  ([(  -  5;(cr(]  | , 
where  S-(<z)  and  C;(<7)  denote  Fresnel  integrals,  defined  as 


5;(0-)  = 


sin  X 


dx. 


(36) 


C,(cr)  • 


cos  X 


dx. 


The  occurrence  of  Fresnel  integrals  here  is  intuitively  ap¬ 
pealing,  because  such  functions  are  known  to  enter  into  the 
evaluation  of  diffraction  effects  in  the  linear  infinite  baffle 
problem.* 


V.  RESULTS  FOR  A  PISTON 

Expressions  describing  the  response  in  terms  of  real 
functions  may  be  obtained  from  Eqs.  (32)  and  (34).  In  view  of 
Eqs.  (35),  different  expressions  are  required  for  the  propagat¬ 
ing  modes  (n  <  1)  and  the  evanescent  modes  (n  >  1).  The  sin¬ 
gularity  in  Eq.  (34)  at  n  =  1  is  regular.  Numerical  evaluation 
of  the  pressure  integral  therefore  presents  no  unforeseen  dif- 


(21 1 


J  Acoust  Soc  Am  ,  Voi  76,  No.  4,  October  1984 


Jerry  H  Ginsberg:  Beams  at  finite  amplitude  II.  Uniform 


1211 


Acuities.  As  is  usual  for  diffraction  integrals,  the  pressure 
integrand  in  Eq.  (34)  oscillates  rapidly  as  a  function  of  n  at 
large  distances  from  the  transducer  (large  r  or  R  |.  The  tran¬ 
scendental  nature  of  the  coordinate  transformation,  Eqs. 
(32l,  at  a  specified  n  makes  the  usual  n’  '.nerical  algorithms 
for  efficiently  evaluating  such  integrals  unsuitable.  How¬ 
ever,  a  Gauss-Chebychev  integration  formula’  is  particular¬ 
ly  well-suited  to  the  \/n„  singularity.  This  matter,  as  well  as 
the  algorithm  by  which  the  coordinate  transformation  is 
evaluated  at  discrete  values  of  n,  will  be  described  elsewhere. 

The  net  result  is  an  algorithm  that  is  relatively  costly  for 
extensive  computation.  Nevertheless,  it  yields  a  prediction 
of  the  instantaneous  acoustic  pressure  that  can  be  utilized  to 
generate  waveforms  or  spatial  profiles  of  the  signal. 

Such  predictions  are  limited  to  locations  where  shocks 
do  not  form  in  the  individual  duct  modes.  Shocks  are  mani¬ 
fested  by  a  vanishing  value  of  the  Jacobian  of  the  coordinate 
transformation,  corresponding  to  a  multivalued  solution. 
Selecting  the  appropriate  solution  in  the  presence  of  a  shock 
requires  considerations  not  addressed  in  the  present  study. 

The  experiments  recently  reported  by  Gallego-Juarez 
and  Gaete-Garreton'"  for  propagation  in  air  provide  useful 
data  for  validating  the  analysis.  The  transducer  for  those 
experiments  was  a  circular  plate  whose  spatial  vibration  pat¬ 
tern  was  combined  with  steps  on  the  plate  in  order  to  simu¬ 
late  a  piston.  For  a  piston  of  diameter  2a,  the  shape  function 
/\R  )  in  Eq.  (1)  is 


T. 

0, 


R  <  ka, 
R  >  ka. 


(37) 


w  hich  leads  to  the  following  Hankel  transform  amplitude; 


f,  =  [kanin)J,\n  ka).  (38) 

The  transducer  diameter  in  the  aforementioned  study 
was  200  mm  and  the  frequency  was  20.4  kHz.  This  corre¬ 
sponds  to  ka  =  37. 1  when  Cg  =  345  m/s.  (The  ambient  con¬ 
ditions  were  not  specified.)  Only  measurements  beyond  a 
distance  2/fc  =  900  mm  were  reported,  whereas  the  last  axial 
maximum  on  axis  occurred  at  2/k  =  218  mm.  Comparing 
the  experimental  and  computed  results  will  therefore  indi¬ 
cate  how  well  farfield  propagation  properties  are  predicted. 

Figures  1-3  reproduce  the  amplitude  response  curves  in 
Ref.  10  for  the  fundamental  frequency  and  the  first  three 


SCLjRCE  SP_  (dB) 


FIG.  2.  Amplitude  response  at  2  ro  on-axis.  Harmonic  number  as  indicated. 

*  =  371  m“',  0  =  6.1  m,  /3„=1.2.  . :  Theory  for  source 

SPL  =  132.5  dB. - :  Ref.  10. . :  Linearized  response. 


higher  harmonics  obtained  at  three  locations  on  the  axis  of 
symmetry.  The  sloping  dotted  line  is  the  closed  form  linear 
solution  of  the  Rayleigh  integral  for  the  pressure  amplitude 
on-axis.'’"  Comparing  the  measured  fundamental  at  low 
source  pressure  levels  with  the  Rayleigh  prediction  for  Fig.  2 
causes  concern.  Mechanisms  such  as  dissipation  which  af¬ 
fect  the  fundamental  at  very  low  source  levels,  where  linear 
theory  is  valid,  do  not  explain  why  the  measured  amplitude 
should  be  higher  than  the  theoretical  one.  Also,  increasing 
the  source  level  in  the  linear  domain  (<115  dB)  did  not  ex¬ 
actly  increase  the  received  fundamental  SPL  in  the  experi¬ 
ments  by  the  same  amount. 

Because  of  this  uncertainty  regarding  the  measure¬ 
ments,  comparisons  with  the  present  analysis  based  on 
source  SPL  might  be  erroneous.  The  method  used  for  com¬ 
parison  was  to  evaluate  the  analytical  waveform  at  a  nominal 
source  level.  The  waveform  was  then  Fourier  analyzed  to 
determine  the  corresponding  amplitude  levels  (and  phase  an¬ 
gles)  for  all  harmonics.  The  comparable  experimental  data 
point  was  selected  by  matching  the  predicted  and  measured 
amplitudes  of  the  fundamental.  This  matching  is  indicated 
by  the  dotted  horizontal  lines  marked  by  the  number  1. 

Figure  I,  which  is  for  the  most  remote  location,  shows 
that  the  second  and  third  harmonic  levels  agree  to  within  2 
dB  between  theory  and  experiment,  while  the  fourth  har- 


'  ^35  155 

SCARCE  SPL  CdB) 


FIG  I.  Amplitude  response  at  3.2  mon-axis.  Harmonic  number  as  indicat¬ 
ed.  *  =  371  m"'.  a  =  0,1  m,  . .  Theory  for  source 

SPL  =  132.5  dB. - :  Ref.  10. . :  Lineanzed  response. 


SOURCE  SPL  (dB) 


FIG.  3.  Amplitude  response  at  1 .2  m  on-axis.  Harmonic  number  as  indicat¬ 
ed  *  =  371  m"'.  0a  -  . ;  Theory  for  source  SPL  =  142.5  dB 

- -  Theory  for  source  SPL  =  132.5  dB. - .  Ref  10. . 

Linearized  response. 
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monies  are  even  closer.  The  discrepancies  between  theory 
and  experiment  are  slightly  larger  in  Fig.  2,  but  they  are  well 
within  the  uncertainty  associated  with  the  difference 
between  the  measured  fundamental  and  the  prediction  of 
linear  theory  at  low  levels.  It  is  possible  that  the  disagree¬ 
ment  results  from  the  experimental  configuration,  which 
employed  a  stepped  transducer  face.  Even  though  the  funda¬ 
mental  was  believed  to  match  well  with  an  ideal  piston," 
minor  discrepancies  might  have  a  substantial  effect  on  the 
higher  harmonics.  Also,  it  is  conceivable  that  the  motion  of 
the  transducer  face  was  not  exactly  axisymmetric.  This 
could  substantially  alter  the  on-axis  diffraction  effects. 

Figure  3  presents  the  response  curves  for  the  closest 
location  in  the  experiments.  As  indicated  by  the  amount  of 
reduction  in  the  fundamental  amplitude,  most  measure¬ 
ments  of  the  third  and  fourth  harmonics  were  taken  at 
source  levels  for  which  the  effects  of  shock  formation  are 
significant.  The  coordinate  transformation  for  a  range  of  val¬ 
ues  of  the  transverse  wavenumber  n  has  multiple  solutions 
when  shocks  are  present.  The  theory  is  not  formally  valid  in 
this  case,  but  the  numerical  algorithm  was  implemented  to 
select  the  value  of  the  strained  coordinate  closest  to  the 
value  of  2  when  a  shocks  occurs. 

Two  source  levels  were  evaluated  for  Fig.  3.  The  dashed 
lines  for  a  theoretical  142.S  dB  SPL  correspond  to  significant 
shocking  effects.  The  agreement  between  theory  and  experi¬ 
ment  for  the  lowest  three  harmonics  is  remarkable,  especial¬ 
ly  in  view  of  the  uncertainty  about  how  shocks  should  be 
treated  theoretically.  The  broken  line  for  a  theoretical  SPL 
of  132.5  dB  show  the  same  degree  of  agreement  for  the  sec¬ 
ond  harmonic.  Extrapolating  the  curve  for  the  third  har¬ 
monic  back  to  this  level  would  show  that  the  theory  closely 
predicts  this  harmonic  also. 

Figures  1-3  do  not  provide  a  complete  picture.  First, 
they  describe  locations  which  are  in  the  farfield.  Further¬ 
more.  the  response  curves  do  not  display  the  phase  angles  for 
the  various  harmonics.  The  higher  harmonics  for  planar 
waves  are  in-phase  with  the  fundamental,  as  are  cylindrical 
and  spherical  waves  in  the  farfield.  Asymmetry  between  the 
compression  and  rarefaction  phases,  which  was  observed  by 
Mellen  and  Browning,^  corresponds  in  the  frequency  do¬ 
main  to  out-of-phase  conditions. 

Figures  4  and  5  display  waveforms  at  axial  locations  in 
the  nearfield  and  in  the  farfield,  respectively.  Reference  10 


’  !  ME  Cr'Of^d  1  “n.  > 


FIG  4.  Waveform  at  0.184  m  on-axis  for  142.5  dB  source  SPL  —  371 
m  a  =i  0. 1  m.  =  1.2. - :  Theory. .  :  Linearized  response. 
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FIG.  5.  Waveform  at  2.1  m  on-axisfor  133.5  dB  source  SPL  =  571  m  " 

a  =  0.1  m,  ^0  =  1-2- -  Theory. . :  Measured  in  Ref.  10  for 

130.5  dB  source  SPL. . ;  Linearized  response. 


measured  a  waveform  at  the  same  location  as  that  described 
by  Fig.  5,  but  for  a  source  level  of  1 30.5  dB.  The  theoretical 
waveform  obtained  from  the  present  analysis  for  the  same 
source  level  would  show  less  distortion  than  the  observed 
one.  This  discrepancy  is  consistent  with  those  discussed  in 
regard  to  Fig.  2,  which  describes  a  nearby  location.  For  this 
reason,  the  theoretical  source  SPL  was  increased  by  3  dB  for 
the  comparison.  The  source  level  for  Fig.  4  was  selected  to 
give  a  comparable  amount  of  distortion  at  the  closer  loca¬ 
tion. 

Several  features  of  the  measured  waveform  in  Fig.  5 
must  be  noted.  Although  an  interval  of  two  periods  is  shown, 
the  shape  for  the  second  period  does  not  duplicate  that  for 
the  first.  Also  the  null  pressure  level  was  not  indicated  in 
Ref.  10;  the  experimental  waveform  reproduced  in  Fig.  5  has 
been  plotted  to  give  a  zero  mean  value.  Finally,  the  original 
waveform  was  drawn  to  a  very  small  scale,  so  its  enlarge¬ 
ment  to  obtain  Fig.  5  may  have  introduced  additional  inac¬ 
curacies.  In  view  of  these  uncertainties,  the  agreement 
between  the  experimental  and  theoretical  waveforms  is  quite 
good,  particularly  for  the  second  period. 

The  distortion  of  the  waveform  in  Fig.  4  is  the  type 
reported  by  Browning  and  Mellen.^  One  effect  of  nonlinear¬ 
ity  is  to  shift  the  extrema  in  the  manner  that  a  plane  wave 
distorts.  Nonlinearity  also  enhances  and  narrows  the  com¬ 
pression  phase,  while  it  has  the  opposite  effect  on  the  rarefac¬ 
tion  phase. 

The  waveform  in  Fig.  5  is  for  a  lower  source  level  than 
the  one  for  Fig.  4,  so  less  of  the  distortion  is  associated  with 
the  neaifield.  This  gives  rise  to  a  phenomenon  that  is  some¬ 
what  different  from  the  one  appearing  in  Fig.  4.  It  is  instruc¬ 
tive  to  compare  the  waveforms  in  both  figures  to  the  corre¬ 
sponding  linearized  signals.  It  seems  that  the  maximum 
compression,  which  occurs  sooner  than  the  linear  signal  in 
the  nearfield,  is  retarded  as  it  propagates  until  it  matches  the 
linear  maximum.  Also,  the  maximum  rarefaction  is  retarded 
more  in  the  farfield.  Recall  that  in  the  Rayleigh  formulation, 
the  individual  harmonics  appear  to  be  nonuniform  spherical 
waves  in  the  farfield.  The  aforementioned  retardation  effect 
might  be  a  consequence  of  the  90°  phase  shift  that  spherical 
waves  undergo  in  the  transition  from  the  neatfield  to  the 
farfield.  In  addition,  the  fact  that  Fig.  4  describes  a  nearfield 
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TABLE  I.  Founcr  senes  data  for  the  waveforms  m  Figs.  4  and  5 


Location 

Harmonic 

number 

10'  p 

p<Fi 

Amplitude 

received 

SPL  IdBI 

Phase  lag 

re\  fundamental 

:/k  =  0.184  m 

linear 

2.646 

142  5 

-4  1* 

1 

2.479 

142  0 

■> 

0.426 

126.7 

6.9‘ 

3 

0.208 

120.4 

54.8’ 

4 

0  083 

112.4 

47.2* 

5 

0  055 

108  8 

66.6' 

z/k  =  2  1  m 

linear 

0.4010 

126.1 

-IXf 

I 

0.3384 

124.7 

0.0623 

110.0 

64  2* 

3 

00247 

101  9 

122.3‘ 

4 

00103 

94.3 

-  177  1* 

5 

0.0044 

86.9 

-  129  6‘ 

location  suggests  that  the  asymmetry  is  produced  by  diffrac¬ 
tion  associated  with  cancellations  of  the  wavelets  emanating 
from  the  vanous  points  on  the  piston.  The  generality  of  these 
observation  needs  to  be  explored  further. 

The  waveforms  were  Fourier  analyzed  using  a  retarded 
time  in  which  the  fundamental  is  a  pure  sine.  The  form  of 
this  representation  is 

p=  Vp,  sin[;it  -  to) 

-  180*<t^^<180‘. 

The  amplitudes  p  and  the  phase  lags  of  the  higher  har¬ 
monics  relative  to  the  fundamental  arc  presented  in  Table  I. 
There  is  no  obvious  pattern  to  the  phase  lags  in  the  nearfield 
case,  but  the  farfield  location  shows  that  the  increment  in  the 
phase  lag  from  a  harmonic  to  next  higher  one  is  nearly  uni¬ 
form.  This  same  pattern  was  observed  in  the  computed 
waveforms  at  other  farfield  locations,  with  the  increments 
ranging  between  40°  and  70”.  The  underlying  mechanism  for 
this  effect  is  open  to  conjecture  at  this  time. 


VI.  CONCLUSION 

The  theory  for  finite  amplitude  sound  beams  developed 
in  Part  1  and  here  provides  a  versatile  algorithm  for  evaluat¬ 
ing  the  effects  of  nonlinearity.  The  theory  has  been  shown  to 
agree  well  with  experimental  data  for  on-axis  responses. 


both  qualitatively  and  quantitatively.  The  derived  expres¬ 
sions  are  applicable  to  off-axis  responses,  but  no  such  results 
have  been  computed  thus  far. 

Comparisons  with  other  available  experiments  are  now 
underway,  but  past  experiments  have  been  limited  in  scope. 
Measurement  of  amplitude  levels  in  the  higher  harmonics 
has  received  more  attention  than  measurement  of  the  corre¬ 
sponding  phase  angles  and  overall  waveforms  (for  under¬ 
standable  reasons).  Also,  measurements  in  the  true  nearfield 
have  been  sparse.  Hence  the  theory  has  not  yet  been  fully 
confirmed,  but  insight  into  basic  phenomena  has  already 
been  gained. 
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The  propagation  of  finite  amplitude  waves  radiating  from  a  baffled 
piston  has  been  described  in  terms  of  directional  spherical  waves  [J.  C. 
Lockwood,  T.  G.  Muir,  and  D.  T.  Blackstock.  J.  Acoust.  Soc.  Am.  53. 
1 148-1  153  (1973)].  That  analysis  predicts  waveforms  at  large  distances, 
provided  that  comparable  information  is  known  at  a  reference  location  in 
[he  farficld.  Lockwood  ct  al.  used  this  approach  based  on  assuming  that 
linear  theory  is  accurate  at  the  reference  location.  Such  an  assumption  is 
inaccurate  when  the  source  pressure  level  is  sutlicicnt  to  generate  signifi* 
cant  nonlinear  effects  (growth  of  higher  harmonics  and  depletion  of  the 
fundamental)  within  the  near  field.  The  present  work  desenbes  the  inter¬ 
lacing  of  the  spherical  propagation  theory  and  the  nonlinear  King  integral 
[J.  If.  Ginsberg,  J.  Acoust.  Soc.  Am.  Suppl.  I  71.  S30(I982)J.  The  latter 
theory  is  used  in  this  approach  to  evaluate  the  nonlinear  waveform  at  the 
reference  location  Companng  the  results  of  interfacing,  and  of  direct 
propagation  according  the  nonlinear  integral  formulation,  with  expen- 
mental  data  provides  a  strong  validation  for  both  theones.  The  advantage 
of  using  directive  sphcncal  wave  theory  lies  m  its  superior  computational 
efficiency  and  its  ability  to  desenbe  shock  formation.  (Work  supported  by 
ONR,  Code  425-UA.  and  NSF,  Grant  MEA-8I01 106.J 
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VV4,  Analysis  of  nonlinear  harmonic  generation  for  arbitrary  dual 
frequency  transducer  excitation.  Mosaad  A.  Foda  and  Jerry 
H.  Ginsberg  (School  of  Mechanical  Engineering,  Georgia  Institute  of 
Technology.  Atlanta,  GA  30332) 

An  earlier  study  of  finite  amplitude  axisymmetric  sound  beams  [J.  H. 
Ginsberg,  J.  Acoust.  Soc.  Am.  Suppl.  1  73,  Sg3  (1983)]  considered  the  case 
of  monochromatic  excitation.  That  work  featured  a  singular  perturbation 
analysis  combining  asymptotic  integration  and  the  renormalization  ver¬ 
sion  of  the  method  of  strained  coordinates.  The  present  paper  initiates  an 
extension  of  those  techniques  to  the  case  of  a  dual  frequency  source.  The 
parametric  array,  in  which  the  pnmary  beams  are  at  closely  spaced  fre¬ 
quencies,  has  already  received  much  attention.  The  system  discussed  here 
permits  disparate  frequencies.  Aside  from  a  restriction  to  axisymmetry, 
the  excitation  at  each  frequency  is  arbitrary.  The  analysis  thus  far  has 
obtained  the  first  two  orders  of  approximation  for  the  velocity  potential. 
This  expression  describes  the  manner  in  which  nonlinear  effects  accumu¬ 
late  for  the  various  sum  and  difference  frequencies.  It  is  the  foundation  for 
a  future  derivation  of  an  expression  for  the  pressure  that  is  descriptive  of 
the  entire  field.  In  addition,  the  trend  for  harmonic  generation  indicated 
by  the  analysis  suggests  that  conversion  efficiency  in  the  parametric  array 
might  be  improved  by  altering  the  transverse  vibration  pattern  of  the 
individual  primary  beams.  [Work  supported  by  ONR.  code  425-UA.l 
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K7.  Evaluation  of  the  overall  sound  field  properties  for  a  finite  amplitude 
sound  beam.  J.  H.  Ginsberg  (School  of  Mechanical  Engineering,  Georgia 
Institute  of  Technology,  Atlanta,  GA  30332) 

The  nonlinear  King  integral  [J.  H.  Ginsberg,  J.  Acoust.  Soc.  Am.  (to 
be  published)]  provides  a  general  algorithm  for  finite  amplitude  axisym- 
metric  waves  radiating  from  a  harmonically  vibrating  transducer.  The 
derivation  of  that  result  was  based  on  asymptotic  analyses  of  the  trans¬ 
verse  wavenumber  spectrum  near  the  axis  for  almost  planar  modes  and  far 
off  axis.  The  validity  of  the  analysis  is  confirmed  here  by  a  change  of 
variables  that  yields  an  overall  measure  of  the  associated  error.  Previous 
evaluations  using  the  nonlinear  King  integral  provided  temporal  and  fre¬ 
quency  spectrum  predictions  at  selected  locations,  primarily  on  axis.  The 
present  paper  reports  on  an  extensive  mapping  of  the  field  for  a  moderate¬ 
ly  high  frequency  in  terms  of  amplitude  and  relative  phase  lags  for  the 
fundamental  and  several  higher  harmonics.  This  mapping  is  cross-refer¬ 
enced  to  waveform  displays  that  show  the  changing  nature  of  the  asym¬ 
metrical  distortion  process  associated  with  transition  from  the  nearfield  to 
the  farfield.  [Work  supported  by  ONR,  Code  425-UA.] 
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3;3< 

V9.  Finite  amplitude  distortion  and  dispersion  in  a  bard-walled 
wareguide.  J.  H.  Ginsberg  and  H.  C.  Miao  (School  of  Mechanical 
Engineering,  Georgia  Institute  of  Technology,  Atlanta,  GA  30332) 

The  fundamental  symmetric  two-dimensional  mode  in  a  hard-walled 
rectangular  waveguide  is  decomposed  into  a  pair  of  obliquely  propagating 
planar  waves,  in  order  to  treat  the  effect  of  nonlinearity.  A  perturbation 
analysis  of  the  nonlinear  wave  equation  for  the  velocity  potential  identifies 
amplitude  dispenion  as  one  source  of  singularity.  The  only  interaction 
between  the  oblique  waves  attributable  to  this  effect  is  a  change  in  the 
distance  parameter  affecting  the  magnitude  of  the  higher  harmonics.  An¬ 
other  singularity  arises  when  the  frequency  u  or  width  L  is  large.  The 
oblique  waves  are  then  closely  aligned  with  the  axis,  resulting  in  resonant 
interaction  with  the  true  planar  mode.  Harmonic  generation  m  this  case 
has  the  appearance  of  a  spatial  beating  pattern.  A  set  of  coordinate  trans¬ 
formations  make  the  representation  uniformly  valid.  Analyses  of  limiting 
forms  are  confirmed  by  quantitative  examples.  Small  values  of  luL  are  well 
described  by  an  earlier  general  solution  in  terms  of  groups  of  nondisper- 
sive  modes  [J.  H.  Ginsberg,  J.  Acoust.  Soc.  Am.  4S,  1127-1133  (1979)], 
while  large  toL  yields  a  quasiplanar  signal.  The  transition  at  moderate  uL 
is  characterized  by  frequency,  as  well  as  amplitude,  dispersion.  The  dis- 
tonion  of  waveforms  then  is  very  close  to  that  obtained  in  the  nearfleld  of 
sound  beams.  [Work  supported  by  ONR.) 


Joe  C.  Thompson  Conference  Center 
The  University  of  Texas  at  Austin 
Austin,  Texas 
8-12  April  1985 
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Finite  Amplitude  Distortion  and  Dispersion 
of  a  Symmetric  Mode  in  a  Waveguide 


Jerry  H.  Ginsberg 
Hsu-Chiang  Miao 

School  of  Mechanical  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 


Abstract 

The  perturbation  method  of  renormalization  is  used  to  study  the  effect  of 
nonlinearity  in  a  hard-walled  waveguide.  The  excitation  would  induce  only  the 
fundamental  symmetric  mode  if  the  system  was  linear.  The  analysis  develops  a 
solution  that  satisfies  a  nonlinear  wave  equation  for  the  velocity  potential, 
as  well  as  all  boundary  conditions.  The  response  consists  of  a  pair  of 
oblique  planar  waves  that  interact  through  second  order  excitation  of  the  true 
planar  mode. 

The  investigation  discloses  that  v^en  the  transverse  width  is  much  larger 
than  the  axial  wavelength,  the  signal  has  a  quasi-planar  behavior.  In 
contrast,  when  the  axial  wavelength  is  large,  the  oblique  waves  are 
essentially  independent.  The  distortion  is  then  a  result  of  self-refraction, 
in  which  the  particle  motion  shifts  the  wavefronts  and  rays.  The  transition 
between  the  long  and  the  short  wavelength  approximations  is  marked  by  the 
appearance  of  nonlinear  frequency  dispersion  which  produces  asymmetrical 
distortion  of  the  waveform. 


1.  INTRODUCTION 


Finite  amplitude  effects  in  a  waveguide  feature  multidimensional 
phenomena  involving  interacting  waves.  In  linear  theory  a  mode  in  a 
hardwalled  waveguide  may  be  constructed  from  pairs  of  oblique  planar  waves 
that  are  reflected  from  the  walls.  The  present  study  will  employ  the  same 
type  of  decomposition  to  show  that  distortion  resulting  from  nonlinearity 
displays  a  phenomenological  change  as  the  excitation  frequency  is  increased. 
This  transition  is  associated  with  an  anomaly  contained  in  previous  studies, 
which  only  considered  the  low  frequency  case. 

Initial  explorations  of  finite  amplitude  nonplanar  modes  in  waveguides 
employed  the  perturbation  method  of  multiple  scales  in  a  rudimentary  fashion 
that  considered  selected  aspects  of  wave  interaction  [1-3].  A  different 
method  of  vestigation  was  developed  to  study  waves  radiating  from  a  flat 
plate  [4-9].  To  a  certain  extent  the  latter  studies  were  academic  in 
nature.  The  system  they  treated  featured  a  periodically  supported  plate  of 
infinite  extent.  They  assumed  periodicity  of  the  signal  parallel  to  the 
plate,  which  meant  that  energy  was  propagating  inward  from  infinite 
boundaries.  This  apparent  violation  of  the  uniqueness  condition  nevertheless 
proved  to  be  instructive,  because  the  system  could  be  studied  by  a  variety  of 
analytical  techniques.  The  perturbation  methods  of  multiple  scales  and  of 
renormalization,  and  the  method  of  characteristics  mutually  agreed  for  the 
case  of  a  spatially  sinusoidal  excitation.  One  significant  aspect  of  their 
result  was  the  prediction  of  self-refraction,  in  which  the  wavefronts  and  rays 
of  constant  phase  are  distorted  by  the  particle  velocity. 

Although  the  plate  problem  did  not  treat  a  physically  realizable  system, 
the  relevance  of  these  investigations  to  waveguides  was  recognized  in  a 
subsequent  investigation  [10].  The  basis  of  that  work  was  that  there  are 


nodal  lines  in  the  plate  system  along  vvhich  the  velocity  component  parallel  to 
the  surface  of  the  plate  vanishes.  Such  lines  are  perpendicular  to  the  plate, 
as  they  are  in  linear  theory.  This  observation  led  to  the  conclusion  that  the 
infinite  plate  analyses  had  actually  derived  a  single  mode  in  a  waveguide. 

The  treatment  of  general  excitation  in  a  waveguide  performed  in  Reference 
[10],  which  was  a  straightforward  extension  of  the  method  of 
renormalization,  disclosed  a  type  of  superposition  principle.  Modes  having 
identical  phase  speed  were  found  to  form  distinct  groups  whose  distortion  in 
self-refraction  was  a  consequence  of  only  the  particle  velocity  arising  from 
that  group.  The  overall  response  consisted  of  a  linear  combination  of  the 
response  in  each  group. 

A  similar  analysis  had  been  used  to  study  waves  radiating  from  cylinders 
[11-14].  One  of  those  studies  [12]  identified  a  paradox  associated  with  very 
long  axial  wavelengths.  One  would  expect  that  if  the  wavelength  along  the 
axis  of  a  cylinder  is  large,  so  that  the  rate  of  variation  in  that  direction 
is  very  gradual,  then  the  response  would  approach  that  for  the  case  of  a  two- 
dimensional  system,  in  which  the  axial  wavelength  is  actually  infinite.  This 
was  found  to  be  the  case,  except  that  the  distortion  phenomena  in  the  limit 
were  found  to  be  too  weak  by  a  factor  of  one  half.  This  dilemma  was  resolved 
by  noting  that  distinct  modes  in  the  case  of  axial  variation  coalesce  only 
when  the  wavelength  is  actually  infinite. 

These  observations  also  apply  to  the  investigation  of  waveguides  [10]. 
For  example,  as  the  width  of  a  waveguide  is  increased,  the  earlier  analysis 
predicts  that  the  distortion  of  the  planar  mode  will  be  twice  as  strong  as 
that  of  the  fundamental  symmetric  (2,0)  mode.  Although  the  explanation  of 
coalescing  effects  for  infinite  transverse  wavelength  (i.e.,  the  planar  mode) 
is  plausible,  it  nevertheless  is  unsettling  from  a  physical  viewpoint. 
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Distortion  arises  from  higher  harmonic  sources  that  are  generated  by 
nonlinearity  in  the  entire  acoustic  field.  Could  it  be  that  minor 

discrepancies  between  the  long  and  infinite  wavelength  cases  accumulate  to 
create  the  discrepancy?  Lack  of  experimental  data  prevented  an  earlier 
response  to  this  question,  but  discussions  with  researchers  currently  involved 
in  such  activity  [15]  sparked  the  present  authors'  interest  in  exploring  these 
concerns. 

The  analysis  presented  herein  treats  an  excitation  of  only  the  (2,0)  mode 
in  a  hard-walled  waveguide.  It  will  be  shown  that  this  mode  excites  the 
planar  mode  in  an  insignificant  fashion,  unless  u  L/c^  »  2it,  where  L  is  the 
transverse  width,  u  is  the  (circular)  frequency,  and  Cq  is  the  linear  speed  of 
sound.  The  phase  speed  of  the  (2,0)  mode  then  differs  slightly  from  that  of 
the  planar  mode.  This  sets  up  a  spatial  beating  phenomenon  that  leads  to  a 
smooth  transition  to  the  planar  mode  response  in  the  manner  one  would 
expect.  The  analysis  will  confirm  the  earlier  theory  for  waveguides 

when  a)L/c„  is  not  large.  It  will  also  show  that  the  transition  from  the 
earlier  theory  to  the  short  wavelength  case  is  marked  by  frequency  dispersion, 
in  which  the  waveforms  are  remarkably  similar  to  those  observed  in  the 
nearfield  of  intense  beams  of  sound  [16]. 

2.  FORMULATION 

A  pressure  excitation  of  the  fundamental,  symmetric,  two-dimensional  mode 
in  a  hard -walled  waveguide  may  be  written  as 

P  I2  =  0  “  ^*’0*^0  sin((i)t)cos(k^x),  e  «  1,  -  L/2  <  x  <  L/2  (1) 


where  is  the  ambient  pressure,  Cq  is  the  speed  of  sound  at  ambient 


conditions,  and  the  transverse  wave  number  is  related  to  the  duct  width  L 

by 


=  2WL  (2) 

The  question  to  be  addressed  here  is  the  effect  of  nonlinearity  associated 
with  the  finiteness  of  e  on  the  waves  that  propagate  in  the  positive  z 
direction  as  a  result  of  this  excitation. 

The  equations  of  continuity,  momentum,  and  state  may  be  combined  to  form 
a  single  nonlinear  wave  equation  governing  the  velocity  potential  [17]  under 
isentropic  conditions. 

2  2 

c/  '7^  -  ^  =  1^  [-^2  (3^  <  1)  (  If  )  +  •  7*]  +  0(*^)  (3) 

3t  Cq 

where  the  nonlinearity  coefficient  is  the  constant  associated  with  the 

second  order  term  in  a  polynomial  expansion  of  the  pressure  perturbation  p  as 
a  function  of  the  density  perturbation  p  at  fixed  entropy. 

P/(Po  ^0^  =  p^Pq  +  •••  (4) 

The  pressure  is  related  to  the  potential  by 


/ 


iL 


If  +  I  7^  •  7^  =  0 


(5) 


From  equations  (4)  and  (5)  p,  p  and  ^  have  the  same  order  of  magnitude,  so 
elimination  of  p  from  these  relations  yields 


The  boundary  conditions  for  ^  are  obtained  by  making  the  particle 


velocity  normal  to  the  walls  vanish. 

If  =  0  at  X  =  (7) 

as  well  as  by  matching  Eq.  (6)  at  z  =  0  to  Eq.  (1).  Also,  for  uniqueness,  it 
is  required  that  the  signal  consist  of  a  wave  propagating  in  the  positive  z 
direction. 

The  initial  stage  of  the  solution  technique  employs  a  regular 
perturbation  expansion  of  the  potential  in  terms  of  the  small  parameter  e. 


(j»  =  E 


(8) 


Matching  like  powers  of  e  in  the  differential  equation  and  boundary 
conditions  leads  to  a  sequence  of  equations  in  the  usual  manner.  The  order  e 


terms  are 


=  i  -T-  -  k  x)] 

'  z  =  0 

+  exp  [i (tut  +  ^ 


where  CC  in  general  shall  denote  the  complex  conjugate  of  all  preceding 

2 

terms.  The  order  e  perturbation  equations  are 


2 

o  O  ^  ‘f’?  a  1  ^ 

s  ^  ^2- -r-k 

at  c^ 


(10a) 


X  =  +  L/2 


(10b) 
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^  Z=0  "  ^  ^  •  V4.i]| 


z  =  0 


(10c) 


3.  EVALUATION  OF  THE  POTENTIAL 

It  is  a  straightforward  matter  to  solve  Eqs.  (9)  by  separation  of  variables, 
with  the  result  that 


4)1  =  -^  (exp  [i(u)t  -  k^  X  -  k^  z)] 


+  exp  C  i  (u)t  +  k^  X  -  k^  z)]}  +  CC 


where 


k  =  u/c^  ,  9  *  sin"^  (k„/k),  k,  =  k  cos  9  =  (k^  -  k  (12) 
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Only  the  case  of  propagating,  rather  than  evanescent,  waves  is  of  interest, 
which  means  that  <  k.  This  condition  is  obtained  whenever  u  exceeds  the 
cutoff  frequency  for  the  fundamental  mode,  w  >  2ti  • 

Equation  (11)  represents  the  first  order  solution  as  two  trains  of  planar 
waves  propagating  symmetrically  relative  to  the  centerline  x  =  0.  These  waves 
are  depicted  in  Figure  1,  where  e^  and  are  the  individual  directions.  The 
angle  9  measures  the  direction  in  which  these  waves  propagate  relative  to  the 
centerline.  Each  wave  represents  the  reflection  of  the  other  from  the  rigid 
walls.  Increasing  either  the  frequency  u  or  width  L  decreases  9.  In  the 
limit  9+0,  the  two  trains  of  waves  merge  into  the  planar  mode. 

The  first  step  in  deriving  is  to  use  Eq.  (11)  to  form  the 
inhomogeneous  terms  in  Eq.  (10a).  This  yields 


{exp  [2i(u)t  -  k^x  -  k^z)] 


+  exp  [2i  (ut  +  k^x  -  k^z)]} 

-  J  -  2  ^  )exp  [2i(»t  -  k^z)]  +  CC  (13) 

The  first  two  exponentials  in  Eq.  (13)  excite  second  harmonics.  Such  signals 
propagate  parallel  to  the  two  waves  forming  iji^,  which  are  homogeneous 
solutions  of  the  linearized  wave  equation.  The  corresponding  particular 


solution  may  be  obtained  by  the  method  of  variation  of  parameters,  in  which 
the  amplitude  of  the  homogeneous  solution  is  considered  to  be  an  unknown 


8 


function.  The  last  inhomogeneous  term  is  a  planar  second  harmonic.  Such  an 
excitation  matches  the  planar  mode  for  the  waveguide  when  =  k.  Hence, 
decreasing  k^  brings  the  planar  part  of  the  excitation  into  close  coincidence 
with  the  planar  mode  for  that  frequency,  which  means  that  this  excitation  is 
nearly  resonant  at  small  k^.  The  method  of  variation  of  parameters  will  also 
yield  the  solution  associated  with  this  term.  Thus,  let 

((>2  =  u(x,z)  exp  (2iu)t)  +  CC 

u  =  C(z)  {exp  [-  2i  (k^z  +  k^x)]  +  exp  [-2i  (k^z  -  k^x)]} 

+  D(z)  exp  (-2ik^z)  (14) 

It  should  be  noted  that  the  unknown  functions  C  and  D  depend  on  the  axial 
distance  only.  The  periodic  nature  of  the  excitation  eliminates  dependence  of 
these  functions  on  t.  Similarly,  the  rigid  wall  conditions,  Eq.  (10b), 
imposed  along  x  =  ^  ir/k^,  could  not  be  satisfied  if  C  or  D  were  functions  of 

X. 

The  result  of  requiring  that  Eq.  (14)  satisfy  Eq.  (13)  is  a  set  of 
uncoupled  differential  equations  for  the  amplitude  functions.  After  Eq.  (12) 
for  k^  is  applied,  these  equations  are  found  to  be 

C"  -  4  ik^  C  =  -  i  iSg  (u 

D"  -  4i  k^  O'  +  4k^^  0  =  -  i  ^ 


(15) 
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where  a  prime  denotes  differentiation  with  respect  to  z. 

The  particular  solutions  of  Eqs.  (15)  are  readily  found  to  be 


$0  0, 

32  ^ 


(16) 


It  is  convenient  to  let  the  constant  coefficients  of  Cp  and  Dp  appear 
explicitly  in  the  corresponding  complementary  solutions,  which  are  therefore 
written  as 


''c  = 


6  (1) 


32 


-  [Cj  +  C2  exp  (4ik^  z)] 


°c  = 


16  ^,^2  ■  ^2^  ^°1 


+  O2  exp  (X2  z)] 


(17) 


where  and  X2  ere  roots  of  the  characteristic  equation 


'1.2  -  '"z  h.2  ♦  “ 


(18a) 


The  roots  are  found  with  the  aid  of  Eqs.  (12)  to  be 


X,  =  2i(k  -  k) ,  X,  =  2i(k,  +  k) 


(18b) 


The  expressions  for  obtained  by  substituting  Eqs.  (16)  and  (17)  into 
Eqs.  (14)  must  satisfy  the  radiation  condition.  In  order  for  (j)^  to 
represent  an  out-going  wave  in  the  z-direction,  it  must  only  contain  negative 
imaginary  exponentials  in  the  z  variable.  Satisfaction  of  this  condition 
requires  that  ^2=02=  0.  The  remaining  terms  yield 

3  ti) 

u  =  32  I  (z  +  C^)  (  exp  C-  2i(k^  ^  '^x 

+  exp  [-2i(k^  ^  ■  ‘‘x 

■"lf“  (-2ik^2)  +  Dj  exp(-2ikz)]  (19) 

*^x  ^ 

Note  that  describes  complementary  solutions  of  the  wave  equation  associated 
with  second  harmonics  of  the  oblique  waves,  whereas  Dj  is  the  planar  eigenmode 
at  the  second  harmonic  frequency. 

The  case  =  0  corresponds  to  a  true  planar  mode,  which  is  governed 
by  the  Earnshaw  solution  for  a  nonlinear  planar  wave.  However, 
letting  k^  +  0  in  Eq.  (19)  results  in  a  singularity  in  the  coefficient  of  the 
last  terms. 

Such  behavior  resembles  the  case  of  resonance  in  a  one-degree-of-freedom 
oscillator  whose  equation  of  motion  is 

2 

X  +  (i)  X  =  F  sin  £jt 


(20) 


seems  to 


When  A -►  (j,  the  amplitude  of  the  particular  solution  for  ft  =  « 
become  infinite.  This  ignores  the  presence  of  the  complementary  solution, 
which  forms  a  beating  response  when  it  is  added  to  the  particular  solution. 
In  the  limit  n  =  u,  the  resonant  response  reduces  to  a  harmonic  at 
frequency  u  whose  amplitude  grows  in  proportion  to  t. 

In  the  same  manner  the  singularity  of  Eq.  (19)  at  ♦  0  be  removed 

by  an  appropriate  selection  of  the  coefficient  of  the  homogeneous  solutions. 
The  coefficient  Cj  is  not  used  for  this  purpose  because  the  singularity  is 
associated  with  the  planar  mode. 

In  order  to  study  0,  the  troublesome  terms  in  Eq.  (19)  are  expanded 
in  a  Taylor  series  about  k^/k. 


k,  H  (k2  .  k^2)l/2  ^  k 


1 

1  X, 


exp(-ik^2)  =  exp  C-i(k  -  j  +  •••] 


ik^ 

=  (1  +  z  +  •••)  exp(-ikz) 


(21) 


The  corresponding  asymptotic  form  of  the  planar  terms  in  Eq.  (19)  is 


To) 
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(— ^  -  [exp(-2ik^z)  +  Dj  exp(-2ikz)] 


•1.2 

i  8.  o  Ik  z 

(  _i-  -  I-  )  [1  +  —  +  D  ]  exp  (-2ikz) 
16  \  2  ^2  2k  1  1 


(22) 


The  singularity  for  k^  +  0  is  cancelled  if  the  leading  term  in  =  -1, 


where  the  coefficient  D*  may  depend  on  in  any  manner  that  satisfies  the 
condition 


lim  D*  , 

~2  ~  ° 
k  +  0  k^ 

X  X 


(23b) 


where  5  is  a  bounded  number.  Similarly  the  coefficient  Cj^  is  restricted  to 
depend  on  k  in  any  manner  that  is  not  singular  as  k  0. 

X  X 

The  second  order  potential  is  now  found  from  Eqs.  (14)  and  (19)  to  be 


6^01 


4>2  =  Cj)exp(2iajt)  [exp(-2i4ij)  +  exp(-2it|)2)] 


io)  /”p  2 

16 


(—2  *  (2i(ot)  {expC-i(4-j  +  *2)] 


+  (-1  +  0  )  exp  [-i(il»j  +  <|»2)k/k^]}  +  CC 


(24) 


where 


s  k,z  +  k  X 
1  z  X 

'*’2  *  “^z^  ■  '^x’^ 

The  foregoing  expression  for  satisfies  the  wall  conditions,  Eqs.  (10b). 
At  this  juncture,  412  does  not  satisfy  the  boundary  condition,  Eq.  (10c), 
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which  specifies  that  there  should  be  no  second  order  contribution  to  the 

pressure  at  z  =  0.  This  condition  could  be  satisfied  by  appropriate  selection 

of  the  coefficients  and  D*.  However,  both  of  these  describe  homogeneous 

solutions  for  and  they  are  not  singular  as  k  +  0.  Thus,  they  represent 

2 

effects  that  are  0(e)  at  all  locations.  In  contrast,  observable  distortion 

phenomena  are  associated  with  second  order  terms  that  grow  with  increasing 

distance.  Therefore,  setting 

Cj  =  0*  =  0  (26) 

leads  to  insignificant  errors.  The  corresponding  potential  function  then 
obtained  from  Eqs.  (8),  (11),  and  (24)  is 

2 

c 

I  =  e  ^  exp(ia)t)Cexp(-i.|»j)  +  exp(-ii|)2)] 

+  e  z  exp(2ia)t)  [exp(-2itjij)  +  exp(-2iij)2)] 

6 

-  ||  (  -|  -  ^)  exp(2ia)t)  {exp[-i(^j  +  n,^)] 

k  k 

X 

-  exp[-i(<|,j  +  t2)kA2]}  +  CC  +  O(e^)  (27) 

2 

where  0(  e  )  refers  to  terms  having  that  order  of  magnitude  at  all  locations. 

4.  EVALUATION  OF  THE  PRESSURE 

Prior  formulations  of  nonlinear  propagation  using  the  velocity  potential 
have  generated  the  potential  in  the  form  of  a  separation  of  variables 
solution.  Specifically,  the  expression  was  a  product  of  functions  of  each 


space  variable  and  time.  In  that  situation,  it  was  necessary  to  consider 
individually  the  state  variables  of  particle  velocity  and  pressure. 


The  present  case  is  different  because  the  potential  is  now  represented  as 

two  planar  waves,  each  of  which  is  described  by  a  single  propagation  distance 

parameter.  In  general,  proper  behavior  of  the  expression  for  pressure  in  a 

simple  planar  wave  ensures  comparable  results  for  the  other  state  variables. 

The  pressure  is  related  to  the  potential  function  by  Eq.  (6).  Omission  of  the 

2 

quadratic  products  in  that  relation  ignores  terms  that  are  uniformly  0(e  )  , 
which  is  comparable  to  the  error  in  Eq.  (27)  for  4,  .  Thus, 


c 

0 

-  ~  ei  exp(iait)[exp(-i4,j)  +  exp(-i,j)2)] 

2 

-  jg  ^  z  exp(2i(ut)Cexp(-2’4»j) 

1  2 

+  exp(-2i.l,2)]  -  8  e  ‘  2)  exp(2iu,t) 

*^x 

X  {expC-i(4.^  +  -  exp 


+  CC  +  0(£^) 


(28) 


2 

The  first  set  of  0(  e  )  terms  grows  with  increasing  z  in  all  cases,  and 
the  second  set  grows  when  k^/k  is  very  small.  Such  functional  behavior  is  a 
result  of  using  z  and  x  as  position  variables,  neither  of  which  consistently 
match  the  spatial  scaling  of  the  nonlinear  processes.  In  order  to  ascertain 
the  correct  spatial  dependence,  a  near-identity  transformation  in  the  form  of 


a  coordinate  straining  is  employed.  A  different  transformation  is  introduced 


for  each  wave  variable  i|)^  and 

2 

The  presence  of  0(e  )  terms  in  Eqs.  (27)  that  depend  on  4;^  + 
indicates  that  the  waves  interact.  Further  examination  of  the  form  of  Eq. 

(27)  suggests  the  trial  transformations 


il"-  =  a-  +  e  CP.(a,,a-)  exp(ia)t) 

w  V  w  ^  ^ 

"^j  (oj  ,02 )  6xp(-i(ijt)]  +  ...;  j  *  1,2 


(29) 


where  the  complex  conjugate  term,  denoted  by  an  overbar,  is  introduced  in 
order  to  ensure  that  the  transformation  is  real.  Substitution  for  and  [ii^ 
in  Eq.  (28),  followed  by  expansion  in  Taylor  series  in  powers  of  e,  yields 

— ^  =  -  {  ei  exp(i(ut)Cexp(-iaj)  +  exp(-ia2)] 

‘^o^'o 

-  j  e^[F^  exp(2iu)t  -  ia^)  +  exp(-iaj)  (30) 

+  F2  exp(2iut  -  ia2)  +  ^3 
2 

-  jg  iBo  ^  z  exp(2iu)t)  [exp(-2iaj) 

2 

+  exp(-2ia2)]  -  ^®o  T  ■  exp(2iu»t) 

*^x 

X  {expC-i(a^  +  02)]  -  exp[-i (a^  +  02)k/k^]} 

+  CC  +  O(e^) 
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The  task  now  is  to  determine  the  functions  and  ^2  cancel  all 

O 

0(e)  second  harmonic  terms  which  grow  with  increasing  z.  For  this,  the  terms 
that  depend  on  +  02  apportioned  equally  between  F^  and  F2.  The 

appropriate  choice  is  found  to  be 


ie.k  ,  .2 

Fj  =  -  j-j —  z  exp{-iaj)  -  j  “2  ■  {exp(-ia2) 

^  *^x 

-  exp  ■^«2  ^ 

iS-k^  1 

F2  =  -  —  z  exp{-ia2)  -  4  “2  -  2)  {exp(-iaj) 

^  ‘‘x 

-  exp  C-ia,  (t - 1)  -  io,  ^  ]} 


These  straining  functions  do  not  cancel  all  0(e  )  terms  in  the 
pressure.  The  remaining  terms,  which  are  created  by  the  complex  conjugates  of 
F^  and  F2,  contain  combinations  of  the  and  02  variables.  Their  presence  is 
not  a  problem,  because  they  are  independent  of  t.  Their  role  is  to  cancel  a 
mean  value  of  the  pressure  that  is  created  by  the  cordinate  transformation. 

It  is  convenient  at  this  juncture  to  write  the  coordinate  transformations 
and  pressure  resulting  from  Eqs.  (28)  -  (30)  in  real  functional  form.  The 
pressure  is  governed  bv 
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V* 

■s.* 


i 


E 

where 
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The  foregoing  relations  fully  define  the  pressure.  The  value  of  p  at 
specified  x,  z,  and  t  may  be  determined  by  solving  Eqs.  (33)  simultaneously 
for  the  values  of  and  t  and  then  using  those  values  to  compute  p.  It 
will  be  noted  that  the  terms  in  which  Oj  and  02  couple  do  not  explicitly  grow 
with  z.  However,  their  magnitude  increases  as  k^/k  +  0,  so  the  spatial 
beating  phenomena  created  by  this  interaction  takes  on  the  appearance  of 
growth  in  the  limit.  This  matter  is  treated  in  detail  in  the  next  section. 


5.  ASYMPTOTIC  TRENDS 

Equations  (32)  and  (33)  are  generally  valid,  but  examination  of  the 
behavior  at  limiting  values  of  k^^/k  provides  important  insights. 

For  k^/k  <<  1  (ujL/c^  >>  2tt),  the  coordinate  transformation  may  be  expanded  in 
a  power  series  in  k^/k.  First,  apply  the  identity  for  the  cosine  of  a  sum  to 
the  last  term  in  Eq.  (33a). 

1  k^ 

4'i  =  <Xi  +  ^  cSq  Z  sin(ujt  -  Oj) 

2 

+  e(8Q  -  2)  sin[a)t  ~  ^2  ~  \  - 02)] 

*^x  ^ 

X  sin  [j  (^  -  l)(a^  +  02)]  (34) 

2  2 

Since  k/k  ->•  1  +  k  /2k  +  ...,  the  leading  terms  in  a  Taylor  series 

Z  X 

expansion  of  Eqs.  (33)  are 

'^1  ~  “1  ^  i  e0o'^^  sin(a)t  -  a  )  +  i  gB  (a,  +  02)sin(ujt  -  02) 


(35a) 
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m 

k-‘ 


When  the  same  operations  are  performed  on  Eq.  (33b),  the  result  is 


:  ii 

I 

■ 

'  ■  < 

-  S-1 


: 


<ii2  ~  ^2  ^  Y  ^®o^^  sin((jDt  -  (X2) 

+  +  02)  sin(ut  -  ttj) 


(35b) 


According  to  these  relations  the  values  of  and  02  niay  be  estimated 


as  +  0(e  kz). 


Hence,  the  factor  e(a^  +  82)  may  be  replaced 


hy  e(>J)^  +  <1^2^  =  2  e  k^z  ,  which  is  approximately  2ekz  because  of  the  smallness 
of  kjj/k.  Thus,  the  coordinate  transformations  have  the  common  limiting  form 


tjii  ~  a^.  +  j  [sin(ait  -  a^^)  +  sin((ot  -  02)] 


^  “2v 


+  eSgkz  sin  (ut  -  - )  cos  { - ^ - ) 


(36) 


from  which  it  follows  that 


-  <i)2  =  2.  k^x  ~  <*1  ”  <*2 


“1  ®2,  /l  ■  =^2, 


+  ii2  =  2  k^z  ~  a^  +  02  +  2£8Qkz  sin  (ut - 2 - )  ( - 2 - ) 


The  same  analysis  is  now  applied  to  Eq.  (32).  Series  expansion  in  powers  of 
k^/k  yields 


4 

•  •  •  v> 


I..'.  V-J 

4 


t  i 


I  i 


•  < 


t  i 


i  i 


2  [sin(u)t  -  +  sin((ut  -  a^)] 


1  ?  k 

+  -g  e  (^0  ~2  ~  COS(aj  -  a2)  -  COS[(aj  -  a2) 

*^x 

u  2  .2 

k  k 


‘1  “2'  2k^  ^  ^ 


01-1  ~  cxp  cx  ^  ^ 

E  sin  (ait- - — -)  cos{-~ — -) 


The  next  step  is  to  substitute  the  first  of  Eqs.  (37)  into  the  foregoing, 
and  to  use  the  resulting  expression  for  p  to  eliminate  ■*■  32  between  the 
second  of  Eqs.  (37)  and  Eq.  (38).  The  pressure  expression  that  is  derived  in 
this  manner  is 


■j  ~  e  sin  (uit  -  k^2  +  B^kzp)  cos  (k^x) 


If  k^  =  0,  this  expression  reduces  to  the  well  known  solution  for  a 
planar  finite  amplitude  wave  at  moderate  amplitudes  [18].  For  very 

small  k^/k,  the  signal  described  by  equation  (39)  is  a  quasi-planar  wave. 
The  distortion  is  measured  by  the  value  of  S^kzp,  the  change  in  the  axial 
phase  variable  from  its  value  ait  -  k^z  in  linear  theory.  The  wave  is  not 
truly  planar  because  the  amplitude  varies  with  transverse  position 
as  cos(k^x).  Comparable  phenomena  are  encountered  in  the  far  field  of 
cylindrical  and  spherical  waves  whose  amplitude  is  not  uniform  in  the 
transverse  direction  [11,19]. 

Suppose  that  the  limits  of  Eqs.  (32)  and  (33)  for  small  k^/k  had  been 
derived  without  considering  the  interaction  terms  (those  containing 
both  and  32).  The  result  would  have  been  the  same,  except  that  in  such 
an  expression  would  have  been  replaced  by  ^  g  .  In  other  words,  half  the 
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nonlinear  effect  when  k  <<  k  is  due  to  interaction  between  the  oblique  waves. 

X 


The  situation  for  comparatively  low  frequencies  (exceeding  cutoff)  can 


also  be  examined  asymptotically.  Suppose  that  k^/k  =0(1)  (recall 


that  <  k  for  propagating  modes).  In  that  case  the  interactive  terms  in 


Eqs.  (32)  and  (33)  are  not  associated  with  beating  interactions,  so  they 


remain  0(e)  at  all  locations.  Such  effects  may  be  ignored.  The  remaining 


terms  may  be  written  as 


P  =  Pi  +  P2. 


^  £  sin(ajt  “  J  “ 


(40a) 


where 


♦j  ■  “j  *  “o  kj  ^  Pj 


(40b) 


The  coordinate  straining  for  each  wave  p.  is  reminiscent  of  that  for  a 

J 


planar  wave,  with  an  important  exception.  The  nonlinear  effect  is  measured  by 


the  difference  between  the  nonlinear  and  linear  spatial  phases,  a.  -  .  In 

<3  3 


an  isolated  planar  wave,  this  difference  is  proportional  to  the  propagation 


distance,  which  would  be  (k^z  _+  k^x)/k  for  waves  propagating  in  the  direction 


of  either  oblique  wave.  Instead,  the  distance  parameter  for  each  wave  in  Eq. 


(40b)  is  z  k/k  .  It  follows  that  although  Eqs.  (39)  specify  a  superposition 


of  the  oblique  waves,  the  presence  of  one  affects  the  other  by  altering  the 


spatial  dependence  for  the  distortion  phenomena. 


Another  viewpoint  for  the  low  frequency  (long  axial  wavelength)  case  may 


be  obtained  from  a  different  resolution.  Define  new  strained 


coordinates  n,  C  such  that 


ai=C''-nt  a2  =  5-  n 


•  'A' 


.I 

A  V-M 
1 


ll 


Return  now  to  Eqs.  (32)  and  delete  the  second  0(e)  term  in  each,  because  those 
terms  are  not  growth  effects  when  k^/k  =  0(1).  The  variables  and  a.^  are 
removed  from  the  functional  dependence  by  forming  the  sum  and  difference  of 
those  equations  after  substitution  of  Eqs.  (41).  This  yields 


1  Ic^ 

^  I  '^^o  IT  ^  sin(ut  -  c)  cos  (n) 


1 

kxX  =  n  -  7  —  z  cos  (o»t  -  c)  sin  (n) 


The  corresponding  expression  for  pressure  obtained  from  Eq.  (32)  is 


=  G  sin(tot  -  5)  cos(n)  +  0(e  ) 


The  significance  of  this  representation  of  the  signal  becomes  apparent 
when  the  particle  velocity  is  evaluated.  For  this,  the  oblique  planar  wave 
decomposition  is  useful.  The  approximation  v  =  p/pc^  is  applicable  to  weakly 
nonlinear,  as  well  as  linear,  planar  waves.  The  propagation 
directions  e’j  and  in  Figure  1  may  be  used  in  conjunction  with  Eqs.  (32) 
and  (41)  to  represent  the  individual  contributions.  Thus 


1  2 
2  =  -2  CoG  [  e j  sin(a)t  -  5  -  n)  +  sin(a»t  -  c  +  n)]  +  0(e  ) 


(44a) 


The  components  of  particle  velocity  are  therefore 


v_  =  V  •  e_  =  c  G  sin(wt  -  c)cos(n) 

2  U 


(44b) 


1  =  -  V  Ic  ■  5)  sin  (n) 
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These  expressions  may  be  substituted  into  Eqs.  (41),  with  the  result  that  the 
new  strained  coordinates  are  found  to  be  governed  by 


z 


1  k''  '^x 


This  form  was  derived  in  the  earlier  analysis  that  assumed  noninteracting 
modes  [10].  Constant  values  of  5  and  n  are  wavefronts  and  rays,  respectively, 
for  the  phase  of  the  wave  in  £q.  (43).  The  velocity  components  transverse  to 
these  lines  are  v^  and  v^,  respectively.  Hence,  the  dependence  of  the 
wavefronts  on  v^,  and  of  the  rays  on  v^^,  was  ascribed  to  self-refraction  in 
the  earlier  work. 


6.  EXAMPLE 

The  trends  identified  in  the  previous  section  indicate  that  at  low 
frequencies  (k^  =  0(k))  the  distortion  process  involves  only  the  harmonics  of 
the  fundamental  mode  for  the  waveguide.  In  contrast,  at  high 
frequencies  (k^  «  k)  the  tendency  is  to  form  a  quasiplanar  wave  that 
propagates  like  the  true  planar  mode.  Identification  of  these  trends  leaves 
the  questions  of  when  the  transitions  to  each  situation  occur,  and  what 
happens  in  the  intermediate  regime? 

These  matters  may  be  addressed  by  numerical  examples.  Quantitative 
results  in  general  are  obtained  by  solving  the  coupled  transcendental  Eqs. 
(33)  for  the  strained  coordinates  and  02  corresponding  to  specified  values 
of  x,  z,  and  t.  These  values  then  yield  the  pressure  according  to  Eq.  (32). 
If  desired,  a  waveform  may  be  generated  by  incrementing  ut  through  an 
interval  2  ir  ,  and  that  result  may  be  Fourier  analyzed  to  determine  the 


*  V  * 


*  •  •  A 


.  t,... 


r  .V**  .* 


fre-i'jency  response.  One  simplification  in  performing  a  numerical  evaluation 

2 

is  that,  for  specified  properties  of  the  fluid,  the  value  of  P/pqC^  obtained 
from  Eqs.  (32)  and  (33)  depends  only  on  the  independent  variables  kx,  kz, 
and  a)t  and  on  the  value  of  kL,  (because  k^/k  =  2ir/kL)  .  For  the  discussion 

3 

that  follows,  the  fluid  is  air  (p^  =  1.2  kg/m  ,  c^  =  343  m/s,  =  1.2) 
and  0)  =  10  kHz. 

A  case  of  comparatively  low  frequency  is  illustrated  in  Figure  2,  for 
which  L  =  0.20  meter  and  e  =  0.0014166  corresponding  to  an  excitation  of  140 
dB  re  20  pPa  at  the  origin.  For  comparison,  the  noninteractive  theory,  Eqs. 
(40),  and  the  quasi-planar  limit,  Eq.  (39)  are  also  shown  in  the  figure.  The 
unimportance  of  the  mixing  between  the  oblique  waves  is  apparent,  as  is  the 
fact  that  the  distortion  associated  with  the  planar  theory  is  stronger. 

Altering  the  frequency  for  the  next  example  would  change  the  overall 
degree  of  nonlinearity.  For  example,  the  distance  for  shock  formation  in  the 
planar  wave  is 

a  =  (46) 

Since  the  degree  to  which  wave  interaction  is  significant  depends 
(nondimensional ly)  only  on  the  value  of  kL,  the  various  phenomena  shall  be 
explored  by  changing  L.  Thus,  the  next  case,  illustrated  in  Figure  3,  is  for 
L  =  2  meters,  with  the  other  parameters  unchanged.  The  quasi-planar 
approximation  is  now  very  close  to  the  new  theory. 

The  situation  for  a  transitional  case  is  shown  in  Figure  4,  which 
corresponds  to  L  =  0.5  meter.  Neither  approximation  is  accurate  here.  The 
difference  between  the  axial  phase  speeds  of  the  planar  harmonic  created  by 
nonlinearity  and  the  true  planar  mode  is  relatively  small.  This  leads  to 
frequency  dispersion  in  combination  with  the  usual  amplitude  dispersion  that 
is  associated  with  a  sawtooth  waveform.  The  effect  is  asymmetrical  between 


compression  and  rarefaction;  it  is  remarkably  similar  to  the  near  field 
distortion  observed  for  baffled  transducers  [16]. 

The  relatively  drastic  transition  from  one  approximate  theory  to  another 
resulting  from  increasing  kL  by  a  factor  of  10  has  a  direct  explanation.  The 
frequency  dispersion  phenomenon  is  attributable  to  spatial  beating  described 
by  the  last  terms  in  the  coordinate  transformations,  Eqs.  (33).  The 
trigonometric  identity  for  the  difference  of  cosines  applied  to  these  terms 
shows  that 

COS(a)t  -  a^)  -  COS  [ojt  “  Oj  -  ("jj - 1)  cij)3 

=  -2  Sin  [(a .  +  aj  -  1)]  sin  [ut 
+  (oj  -  a^)  -  j  (aj  +  )  |-];  i , j  =  1,2,  i*j  (47) 

The  first  sinusoidal  factor  is  independent  of  time;  it  governs  the 
wavelength  of  the  beats.  When  the  argument  of  that  sine  term  is  very  small 
compared  to  w,  the  factor  is  well  approximately  by  (oj  +  a^)(k/k^  -  1). 
Since  and  02  iMy  be  approximated  by  k^z,  small  values  of  the  aforementioned 
argument  correspond  to  cumulative  growth  of  the  frequency  dispersion  effect. 

It  follows  that  the  prominence  of  frequency  dispersion  is  indicated 
by  :i/[2k^z(k/k^  -  1)].  In  contrast,  the  significance  of  the  sawtooth 
distortion  effect  is  measured  by  the  ratio  of  the  axial  distance  z  to  the 
planar  shock  distance  0.  A  comparison  of  the  two  nondimensional  factors 
indicates  whether  frequency  dispersion  will  be  noticeable  in  the  presence  of 
sawtooth  distortion.  Thus,  define  a  beating  parameter  8  according  to 
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-  ee„(kz)^  [ - o-  .  VTO  ] 


(48) 


This  parameter  is  5.08,  0.05,  and  0.798  for  Figures  2-4,  respectively.  Cases 
where  B  is  substantially  greater  than  unity  can  be  anticipated  to  be  well 
described  by  the  earlier  noninteractive  theory  for  duct  modes,  whereas  values 
that  are  much  less  than  unity  will  closely  fit  the  planar  wave  approximation. 

Another  aspect  of  the  distortion  process  is  displayed  in  Figures  5  and  6, 
which  are  waveforms  at  off-axis  locations.  The  line  x/L  =  1/4  is  a  node 
according  to  linear  theory,  as  well  as  the  quasi-planar  nonlinear 
approximation.  However,  Figure  5,  which  corresponds  to  such  a  location,  shows 
that  only  the  odd  harmonics  are  nulled  in  the  oblique  wave  theories.  Hence 
the  fundamental  frequency  of  the  signal  at  the  “nodes"  is  twice  the  excitation 
frequency.  Note  that  both  oblique  wave  theories  indicate  that  the  tendency  to 
form  a  sawtooth  profile  is  still  present. 

The  nulling  of  the  odd  harmonics  was  explained  in  the  earlier  analysis  of 
the  plate  problem  as  being  a  result  of  self-refraction  [4,6].  The  rays  in  the 
non-interactive  theory  were  shown  to  be  distorted  in  the  direction  of  the 
transverse  velocity  component.  This  caused  the  nodal  ray  to  cross  the  axial 
line  of  zero  linearlized  pressure  twice  per  axial  wavelength,  thereby  setting 
up  the  second  harmonic  signal.  It  is  apparent  from  Figure  5  that  this  effect 
also  occurs  in  the  presence  of  frequency  dispersion  resulting  from  interaction 
of  the  oblique  waves. 

A  waveform  for  a  general  location  appears  in  Figure  6.  The  even 
harmonics  are  more  prominent  than  they  were  in  Figure  4  because  the  odd 
harmonics  are  lessened  by  the  proximity  to  the  nodal  line.  This  effect  is 
accompanied  by  amplitude  dispersion,  as  evidenced  by  the  tendency  to  a 
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sawtooth  profile,  and  by  frequency  dispersion,  as  indicated  by  the  asynmetry 
between  compression  and  rarefaction. 

A  different  perspective  is  offered  by  the  amplitude  and  phase 
distribution  curves  in  Figures  7-9.  These  curves  were  obtained  by  Fourier 
series  decomposition  of  the  computer  waveforms  into 

—2—  =  J  sin[n  u(t  -  t^)  -  Xp];  Xj  =  0  (49) 

where  t^  is  the  arrival  time  of  the  fundamental  in  the  interacting  oblique 
wave  theory.  The  amplitudes  p^,  are  displayed  for  the  three  nonlinear 
theories.  However  the  phase  lags  displayed  only  for  the  latest  theory 

--  they  vanish  in  the  other  descriptions  in  which  the  waveform  distorts 
symmetrical ly. 

Although  only  three  harmonics  are  displayed  in  Figures  7-9,  their  trends 
are  also  indicative  of  higher  harmonics.  The  earlier  observation  of  the 
increased  relative  contribution  of  the  even  harmonics  in  the  vicinity  of  the 
''nodal"  line  x  =  L/4  is  evident  in  Figures  8  and  9.  In  addtion.  Figure  7 
shows  that  the  phase  of  each  harmonic  tends  to  lag  behind  that  of  its 
predecessor  by  a  uniform  amount  that  increases  as  the  signal  propagates.  This 
effect  was  also  predicted  for  sound  beams  [20],  whose  waveform  in  the  near 
field  is  much  like  Figure  4. 

7.  CONCLUSION 

The  excitation  of  the  true  planar  mode,  which  provides  a  mechanism  for 
the  interaction  of  the  oblique  waves  forming  the  fundamental  symmetric  mode, 
has  been  shown  to  be  significant  for  large  values  of  kL.  In  the  limit, 
multidimensionality  is  only  manifested  as  sinusoidal  variation  in  the 
transverse  direction,  much  like  the  directivity  factor  for  nonuniform 


spherical  waves  in  the  far  field  [19], 

In  the  earlier  (small  kL)  theory  the  modes  are  formed  from  obliquely 
propagating  waves  whose  interaction  is  only  manifested  by  a  change  in  the 
distance  parameter  governing  the  distortion.  If  each  wave  were  truly 
independent,  that  parameter  would  have  bee.i  the  distance  over  which  the  wave 
had  propagated.  Instead  the  distortion  of  the  oblique  waves  depends  on  the 
axial  distance.  That  theory  has  been  shown  here  to  be  valid  when  the 
underlying  assumption  of  distinct  phase  speeds  is  valid.  In  that  case,  kL  is 
moderately  larger  than  2ir  ,  so  that  the  scales  with  which  the  signal  varies  in 
the  transverse  and  axial  directions  are  comparable.  The  transition  from  small 
to  large  kL  is  predicted  by  the  present  theory  to  exhibit  frequency  dispersion 
that  is  responsible  for  distortion  of  the  waveform  that  is  not  symmetrical 
between  compression  and  rarefaction. 

The  same  mechanism  can  be  expected  to  enter  into  other  situations  in  a 
waveguide.  For  example,  suppose  two  modes  are  excited.  If  they  have 
different  phase  speeds,  they  superpose  according  to  the  noninteractive  theory, 
[10].  If  the  two  modes  have  identical  phase  speeds,  the  modes  combine  to  form 
a  nondispersive  group,  for  which  the  earlier  theory  is  also  valid.  In  the 
transitional  situation,  the  two  modes  have  phase  speeds  that  are  nearly 
identical.  The  interaction  of  such  inodes  may  be  anticipated  to  lead  to 
frequency  dispersion  phenomena  of  the  type  identified  here. 
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FIG.  1  Geometry  of  the  oblique  waves. 

FIG.  2  Waveform  on-axis  at  z  =  3.05  m  for  UO  dB  at  the  origin,  L  =  0.2 

m,  0)  =  10  kHz. - :  Interacting  waves; - :  Noninteractive 

theory;-  —  — :  Quasi-planar  wave. 

FIG.  3  Waveform  on-axis  at  z  =  3.05  m  for  140  dB  at  the  origin,  L  =  2.0 

m,  0)  =  10  kHz. - ;  Interacting  waves; - :  Noninteractive 

theory; -  Quasi-planar  wave. 

FIG.  4  Waveform  on-axis  at  z  =  3.05  m  for  140  dB  at  the  origin,  L  =  0.5 

m,  0)  =  10  kHz.  -  :  Interacting  waves; - :  Noninteractive 

theory;  _  _  Quasi-planar  wave. 

FIG.  5  Waveform  at  x  =  0.125  m,  z  =  3.05  m  for  140  dB  at  the  origin,  L  =  0.5 

m,  u  =  10  kHz. - :  Interacting  waves; - :  Noninteractive 

theory;— - :  Quasi-planar  wave. 

FIG.  6  Waveform  at  x  =  0.1  m,  z  =  3.05  m  for  140  dB  at  the  origin,  L  =  0.5 

m,  0)  =  10  kHz.  -  :  Interacting  waves; - :  Noninteractive 

theory;— - :  Quasi-planar  wave. 

FIG.  7  Axial  dependence  of  frequency  response  along  x  =  0  for  140  dB  at  the 

origin,  L  =  0.5  m,  u  =  10  kHz. - :  Interacting  waves; - : 

Noninteractive  theory;  :  Quasi-planar  wave. 

FIG.  8  Axial  dependence  of  frequency  response  along  x  =  0.1  m  for  140  dB  at 

the  origin,  L  =  0.5  m,  w  =  10  kHz.  -  :  Interacting 

waves; - :  Noninteractive  theory;— - :  Quasi-planar  wave. 

FIG.  9  Transverse  dependence  of  frequency  response  along  z  =  3.05  m  for  140 

dB  at  the  origin,  L  =  0.5  m,  u  =10  kHz.  - :  Interacting 

waves; - :  Noninteractive  theory;—  —  — :  Quasi-planar  wave. 
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Summary 

This  paper  analyses  the  interaction  between  a  vibrating  circular  membrane 
contained  in  an  infinite  baffle  and  the  resulting  sound  field  radiated 
into  a  fluid  medium  contained  in  the  half-space  above  the  membrane.  The 
case  of  resonant  excitation  of  the  membrane  leads  to  nonlinear  coupling 
between  transverse  and  in-plane  displacement.  Nonlinearity  within  the 
fluid  medium  is  described  by  a  recent  general  treatment  of  finite 
amplitude  sound  beams  resulting  from  boundary  motion.  The  vibratory 
response  of  the  membrane  is  evaluated  in  a  perturbation  technique  based  on 
the  modes  of  free  vibration.  The  results  give  amplitude-frequency  rela¬ 
tions  for  the  plate  that  account  for  the  Inertial  and  damping  impedances 
of  the  fluid,  as  well  as  expressions  that  may  be  solved  for  the  pressure 
Signal  in  the  fluid. 
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Finite  Amplitude  Effects  in  Acoustic  Radiation 

Consider  a  circular  membrane  of  radius  a  that  is  fixed  at  its  edges 
to  an  infinite  baffle.  An  excitation  at  frequency  u  close  to  a  natural 
frequency  is  applied  on  one  side,  and  a  fluid  medium  occupies  the  half¬ 
space  on  the  other  side.  The  speed  of  sound  and  density  of  the  fluid  at 
ambient  conditions  are  c,  and  p,,  respectively.  Dimensional  cylindrical 
coordinates  are  R/k  and  z/k  within  the  fluid.which  occupies  z  >  0,  and 
dimensional  time  is  t/o),  where  k  »  co/c,.  A  recent  study  described  the 
effects  of  material  and  convective  nonlinearity  on  the  acoustic  radiation 
resulting  from  an  arbitrary  motion  on  the  boundary  [1].  Suppose  that  the 
normal  velocity  on  the  boundary  is 


'  -  -rr  e  Co  f(R)  exp(it)  ♦  c.c. 

^  z-0 


where  c  <<  1  is  the  acoustic  Mach  number  and  f(R)/R''‘  is  bounded  for 
large  R.  Then  the  pressure  at  any  (R,  z),  omitting  a  mean  value  correc¬ 
tion  term,  was  found  to  be  described  by 


(2) 


2  cp 0^ 0 


r® 

an  . 
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expfit  -  uO  clt"  *  c.c. 


In  general,  c.c.  shall  denote  the  complex  conjugate  of  all  preceding 
terms.  The  parameters  m  and  y  are  transverse  and  axial  wavenumbers, 
respectively,  JgC  )  denotes  the  Bessel  function  of  zero  order,  and  V  is 
the  Hankel  transform  of  the  spatial  pattern  f(R), 

u  =  id  -  m  )  if  m  <  1  &  y  -  (m  -  1 )  ^  if  m  >  1 


V  =  R  f(R)  Jo(mR)  dR  (3) 

The  parameters  (a,  5)  are  strained  coordinates  defined  in  implicit  form  by 

z  =  C  -  ireBg^  [(mV/y)  exp(it)  erfc[(yC)^^^]  c-c}  Jg(na) 

R  =  a  +  TteSo^  l(mV/y)  exp(it)  erfc[(yC)  ]  ♦  c.c}  J,(na)  {^) 

The  evaluation  of  the  response  of  the  fluid-membrane  system  requires  that 

eqs.  (3)  and  (4)  be  interfaced  with  the  equations  governing  the  membrane. 

Equations  of  Motion  for  the  Membrane 


An  elastic  membrane  undergoing  finite  deformation  due  to  a  resonant 
excitation  was  studied  by  Chobotov  and  Binder  [2].  Several  shortcomings 
of  the  earlier  work  shall  be  corrected  here.  Small  errors  associated  with 
using  an  assumed  mode  function  shall  be  addressed  by  employing  the  exact 
Bessel  function  mode.  Furthermore,  the  study  here  shall  describe  the 
situation  for  resonance  of  any  mode,  rather  than  only  the  fundamental. 
The  last  matter  is  that  the  resistive  and  reactive  portions  of  the  fluid 
impedance  will  be  derived  analytically.  In  Ref,  [2]  the  acoustic  im¬ 
pedance  was  based  on  a  low  frequency  approximation  that  did  not  account 
for  diffraction. 

Chobotov  and  Binder  began  with  a  derivation  of  displacment  equations 
of  motion  that  accouhted  for  In-plane  deformation  and  geometrical 
nonlinearity.  Examination  of  these  equations  reveals  that  the  membrane 
displacements  are  scaled  such  that 

w  2  u 

transverse  displacement  -  e  -  ;  in-plane  displacement  »  c  r  (5) 


Based  on  these  definitions  and  the  fact  that  the  acoustic  Mach  number  e  in 
eq.  (1)  is  very  small,  the  equation  for  in-plane  motion  may  be  rewritten 
from  its  original  form  in  [2]  as 
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1  ~  3^u 

2  3R''3R‘'  2R  ^3R-'  “  ^  2  ..2 

C  3t 


(6) 


The  corresponding  equation  for  transverse  motion  is 


1  1_ir  lii  +  if  tr  Ih  in 

R  3R‘  3R  eo  ^  3R  3R 


3w  1  _  f It 

“  3R  2  ^  ^3R^  ^ ' 


+  — ^  [Q(R)  cos(t) 
Eka 


(7) 
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In  the  foregoing  r  =  PoC,  a/Ooh  (h  is  the  thickness  of  the  membrane),  Q(R) 
is  the  radial  profile  of  the  excitation  applied  to  the  membrane,  and 


c^  =  (E/p[l  -  ;  c^  -  (Oo/P)’^^  ;  eo  -  (c„/c^)^  (8) 

Note  for  later  use  that  c  <<  c  .  In  addition,  the  magnitude  of  Q  is 

w  u 

required  to  be  sufficiently  small  to  induce  a  transverse  displacement 
whose  peak  velocity  actually  is  a  small  fraction  of  o^. 


Vibratory  Response 

The  eigenfunctions  for  a  membrane  whose  edges  are  fixed  are 

•  JoCn)  where  n  =  i.R/ka  and  -  0  (9) 

J  J  3 

Some  approximate  values  are  =  2.^05,  •  5.520,  A^  ■  8.65^.  Proximity 

to  resonance  is  specified  by  u  ■  AjC^/a. 

In  general,  the  transverse  displacement  may  be  expanded  in  a  series 

of  modes.  When  the  excitation  is  close  to  one  of  the  natural  frequencies, 

the  corresponding  mode  may  be  expected  to  dominate  the  response.  Since 

2 

the  nonlinear  terms  in  eq.  (7)  are  0(g  ),  the  transverse  motion  in  this 
case  should  satisfy 

w  •  « j  Wj ( t )  ♦  £  ^w ( R ,  t )  (10) 

where  w;R,t)  is  orthogonal  to  iij  over  the  interval  0  S  R  S  ka. 


Substitution  of  eq.  (10)  into  eqs.  (6)  and  (7)  shows  that  all  deriva¬ 
tives  have  comparable  magnitudes.  The  largeness  of  relative  to  c^,  in 
combination  with  the  resonance  of  w,  makes  it  permissible  to  neglect  the 
inertia  of  in-plane  motion.  An  expression  for  u  may  then  be  found 
analytically.  Toward  this  end  the  0(1)  term  for  w  in  eq.  (10)  is  sub¬ 
stituted  into  eq.  (6).  Using  the  chain  rule  to  replace  R  by  n  results  in 
an  inhomogeneous  ordinary  differential  equation  for  u. 


|-(nu)l 
3n  n  3n 


(1  -  v) 


[J>(n)2]l 


This  equation  may  be  integrated  twice,  after  which  the  constants  of  in¬ 
tegration  are  selected  so  as  to  satisfy  the  condition  that  u  =  0  at  R  *•  0 
end  R  =•  ka.  Repeated  application  of  the  recursion  relations  for  Bessel 
functions  yields  the  following  expression. 

^  j  2 

u  =  w.  U 
ka  j 

U  -  ^(1  ♦  v)  Jo(n)  J.(n)  ♦  TTv  n  [J.(A.)^  -  Jo(n)^  -  J,(n)^]  (12) 


Amplitude-Frequency  Relations 

The  equation  for  w^(t)  is  obtained  by  substituting  eqs.  (10)  and  (12) 
into  eq.  (7),  then  applying  orthogonality  of  w  with  respect  to  The 

J 

result  is 


r  i  ^ '  2=0  r 

♦  — —  ^  n  - 5  dn  =  -r—  cos(t)' 

eka  „  2  eka 

0  p,Co 


J  n 


Q  Jo(n)  dn 


where  6  is  a  coefficient  of  nonlinear  elasticity  that  is  found  to  be 
J 


(u  [(1  +  v)  J,(n)^  -  2v  Jo(n)  J,(n)]  +  -n  Ji(n''  )  dn  (1^) 


values  of  this  parameter  for  the  first  three  resonances  when  v  »  0.29  are 

6  •  0.08206,  6,  “  0.1  1996,  6,  »  0.114069  .  Chobotov  and  Binder  used  the 
lad  2 

approximate  fundamental  mode  -  1  -  (R/ka)  to  obtain  6^  =  0.084214. 

Because  of  the  smallness  of  the  nonlinear  term  in  eq.  (13),  w  will 

s) 


I 


be  harmonic  in  the  first  order  approximation.  Let  A  be  the  complex 
amplitude,  and  let  w  denote  the  dimensional  transverse  displacement  found 
from  eq.  (5). 

^  exp(it)  +  c.c.  !  exp(it)  +  c.c.  +  O(e^)  (15) 

The  Hankel  transform  V  of  the  normal  velocity  at  the  membrane  may  be 
evaluated  by  comparing  the  (dimensional)  time  derivative  of  w  to  eq.  (1). 
This  defines  the  shape  factor  f(R)  for  R  <  ka.  For  R  >  ka,  f(R)  =  0 
because  the  baffle  is  stationary.  The  transform  is  found  in  this  manner 
to  be 

'ka 

V  =  -AG;  G=  R  Jo(X.R/ka)  J<,(mR)  dm  (16) 
Jq  J 

The  coefficient  G,  which  may  be  evaluated  in  closed  form,  leads  to  an 
expression  for  p  according  to  eq.  (3)*  The  acoustic  loading  applied  to 
the  membrane  is  readily  found  from  that  expression,  because  5  “  0  snd  a  - 
R  at  2  =  0,  see  eqs.  (4). 


P a^o 


nG 
n  P 


exp(it  -  yz)  J,(mR)  dm  +  c.c. 


The  final  step  leading  to  the  amplitude-frequency  relations  is  to  use 
the  method  of  harmonic  balance  to  equate  all  terms  in  the  equation  of 
motion  that  are  proportional  to  exp(it).  (This  procedure  is  equivalent  to 
removing  secular  terms  in  singular  perturbation  schemes.)  The  equation 
obtained  in  this  manner  relates  the  amplitude  eA  to  the  frequency  w  and 
the  generalized  modal  amplitude  for  the  excitation. 


3  2  * 
35,e^A^A 


'  ^^d'  -  ^m^  -  *  ■  --  3 

j  w  2eoLi<aJ,  (Aj  )  J 


r  ^  Q,  (18) 


» 

where  gA  denotes  the  complex  conjugate.  The  coefficients  7  and  are 

m  d 

reactance  (added  mass)  and  resistance  (damping)  impedance  coefficients, 
respectively,  resulting  from  the  acoustic  radiation. 


Equivalent  expressions  were  derived  by  Bouwkamp  [3]  by  using  a  scattering 
integral . 

The  real  and  imaginary  parts  of  eq.  (18)  may  be  solved  for  the  mag¬ 
nitude  and  phase  of  eA;  the  latter  represents  the  phase  lag  of  the 
response  relative  to  the  excitation.  Then  the  acoustic  signal  may  be 
evaluated  by  substituting  the  value  of  gA  into  eqs.  (16)  to  form  the 
Hankel  transform  V. 

Results 

Figures  1  and  2  depict  the  resistance  and  reactance  coefficients  as  a 

function  of  ka.  The  peaks  in  Y  are  centered  around  ka  =  X .  for  the 

m  j 

respective  resonances,  whereas  rises  almost  linearly  beyond  those 

locations.  Suppose  that  the  fluid  is  water,  in  which  case  c  /Cj  <<  1  (oj 
2 

<<  PoCo  because  o  o  cannot  exceed  the  yield  stress).  Furthermore,  note 

that  ka  =  A .  c  /c„  when  ui  =  ui.,  and  that  ka  >>  2it  for  closely  confined 
j  w  j 

sound  beams.  It  follows  that  resonant  excitations  at  high  ka  values  for 
beam  forming  only  arise  for  very  high  order  modes. 


Fig.  1.  Resistance  Y^  for 
the  three  lowest  modes. 


Fig.  2.  Reactance  Y  for 
the  three  lowest  moSes. 


The  dependence  of  the  harmonic  amplitudes  forming  the  pressure 
waveform  along  the  axis  of  a  sound  beam  is  depicted  in  Figure  3.  Only  the 
fundamental  frequency,  marked  #1,  is  treated  by  a  linearized  analysis. 


amplitudef^lO 


The  corresponding  result  for  a  piston  (uniform  velocity)  Is  shown  In 
Figure  4.  In  both  cases  ka  =  HO  and  e  -  0.00102,  which  produces  a  maximum 
sound  pressure  level  in  the  piston  case  of  250  dB  re  1  yPa.  The  fluctua~ 
tions  in  the  piston  case  are  due  to  diffraction  effects  that  alternately 
reinforce  and  cancel  the  fundamental  frequency  in  the  near  field,  z  < 
{ka)^/2u.  Diffraction  effects  are  much  less  significant  for  the  membrane 
because  the  particle  velocity  is  continuous  across  the  edge,  so  the 
propagation  curves  are  smoother,  and  the  peak  values  are  lessened. 


axial  distance  z 


Fig.  3-  Range  dependence  of 
the  amplitude  of  the  three 
lowest  harmonics  in  a  sound 
beam  generated  by  a  membrane, 
as  a  fraction  of  PqCo** 


fO 


Fig.  H.  Range  dependence  of 
the  amplitude  of  the  three 
lowest  harmonics  in  a  sound 
beam  generated  by  a  piston, 
as  a  fraction  of  PoCj,^. 
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9J0 

S3.  Fourier  series  represenlation  of  finite  nrai-litude  sound  beams.  Hsu- 
Chiang  Miaoand  Jerry  H.  Ginsberg  (School  ofMechanicaJ  Engineering. 
Georgia  Institute  of  Technology.  Atlanu,  GA  30JJ2I 

An  earlier  analysis  of  finite  amplitude  sound  beams  derived  general 
expressions  from  the  behavior  in  the  off-axis  legion  [J.  H.  Ginsberg,  J 
Acoust.  Soc.  Am.  ^6,  I20I-I2I4  (1984)].  An  asymptotic  analysis  of  the 
region  very  close  to  the  beam  axis  confirms  the  earlier  result  for  the  veloc¬ 
ity  potential.  The  signal  is  rewritten  in  a  form  that  makes  the  contribution 
of  each  wavenumber  in  a  rjniinuous  spectrum  appear  to  be  the  sum  of  two 
waves  traveling  transversely,  as  well  as  axially.  The  coordinate  transfor¬ 
mations  required  to  renormalize  this  form  lead  to  a  temporal  Fourier 
senes  ihat  is  reminiscent  of  the  Fuoini  solution  for  finite  amplitude  planar 
waves.  The  comp'cx  amplitude  of  each  harmonic  is  obtained  from  an 
ntegration  over  the  transverse  wavenumber  The  computational  effi- 
cency  of  this  repres.  '■lation  permits  extensive  evaluation  of  propagation 
properties  An  example  c  -npares  the  signal  derived  from  a  piston  to  lhat 
obtained  .'rom  the  one-dimensional  assumption  that  p  =  pen,  on  the 
boundary,  which  has  been  employed  in  prior  investigations  using  approxi¬ 
mate  parabolic  .'quations  [Work  supported  by  ONR,  Code  425-UA.J 
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Session  S.  Physical  Acoustics  IV:  Nonlinear  Acoustics 
Joseph  E.  Blue,  Chairman 

Underu^ater  Sound  Reference  Detachment,  Naval  Research  Laboratory,  P.  0.  Box  8337,  Orlando,  Florida  32856 


S4.  Finite  amplitude  acoustic  waves  generated  by  a  baffled  dual  frequency 
transducer.  Mosaad  A.  Fuda  and  Jerry  H  Ginsberg  (School  of 
Mechanical  Engineering,  Georgia  Institute  of  Technology,  Atlanta,  GA 
J0332) 

In  a  previous  presentation  [J  Acoust.  Soc.  Am.  Suppl.  I  75,  S92 
1 19841],  a  nonuniformly  accurate  expression  for  the  velocity  potential  was 
derived  using  a  singular  perturbation  procedure.  Tlic  linearized  signal 
was  represented  by  a  dual  King  integral,  and  the  cumulative  growth  was 
evaluated  asymptotically.  In  this  paper,  the  renormalization  version  of 
the  method  of  strained  coordinates  is  employed  to  annihilate  the  secular 
terms.  The  result  is  a  uniformly  valid  expression  for  the  acoustic  pressure 
at  all  locations  preceding  the  formation  of  a  shock.  Additional  discussion 
is  devoted  to  a  simpler  model  that  interfaces  the  dual  frequency  King 
integral  in  the  nearheld  with  nonlinear  spherical  distortion  in  the  farfield. 
The  results  obtained  from  both  descriptions  in  the  case  of  a  parametric 
array  (proximate  primary  frequencies)  agree  better  than  earlier  theories 
with  previous  measurements  performed  for  a  wide  range  of  parameters 
[Work  supported  by  ONR,  Code  425-UA.] 
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FINITE  AMPLITUDE  ACOUSTIC  WAVES  GENERATED 


ARBITRARY  EXCITATION  OF  TWO  PRIMARIES 


FIRST  ORDER  SOLUTION  FOR  THE  VELOCITY  POTENTIAL 

DUAL  FREQUENCY  KING  INTEGRAL 
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SECOND  ORDER  POTENTIAL 


SECOND  ORDER  SOURCE  TERMS 
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(4)  Asymptotic  integration  by  Laplace's  method 
(5>  Procedure  as  by  Ginsberg  JASA  <1904) 


RENORMALIZATION  OF  THE  SIGNAL 


coordinate  transformations, 
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QUANTITATIVE  EVALUATION 
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<a)  No  coordinatQ  straining  required. 

Cb)  Input  frequency  response  given  directly, 
(c)  Not  a  thorough  description  of 
distortion  in  the  nearfield. 


METHOD  II  -  UNIFORM  SPHERICAL  PROGAGATION  CUSP) 
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•  ••  Measurements  in  the  First  Tank 
AAA  Measurements  in  the  Second  Tank 
I  ■  Predicted  Results 
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